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1 . 0  RESEARCH  OBJECTIVES 


1 . 1  A  Fast  PN  Acquisition  Scheme 

In  this  project  we  studied  acquisition  techniques  for 
DS  spread  spectrum  signals.  The  technique  utilized  two  thres¬ 
holds  rather  than  a  single  threshold.  As  a  result,  the  time 
to  acquire  was  reduced  by  an  order  of  magnitude  as  compared  to 
the  techniques  proposed  by  other  investigators  who  use  single 
threshold  techniques.  The  procedure  can  be  employed  with  most 
acquisition  techniques  merely  by  replacing  the  single  threshold 
with  the  two  threshold  system. 

1 . 2  Tracking  of  Frequency  Hopped  Spread  Spectrum  Signals 
in  Adverse  Environments 

The  objective  of  this  research  is  to  design  tracking 
strategies  in  adverse  environments  such  as  fading  and  jamming. 
This  led  to  optimization  of  tracking  in  the  sense  of  maximizing 
the  mean  time  to  loss  of  lock. 

1 . 3  Comparison  of  Schemes  for  Coarse  Acquisition  of  FH 
Spread  Spectrum  Signals 

The  objective  of  this  research  was  to  compare  three  of 
the  more  popular  acquisition  techniques  on  the  basis  of  prob¬ 
ability  of  miss,  time  to  acquire  and  complexity. 

1 . 4  Two-Dimentional  Optical  Filtering  of  1 -D  Signals 

The  object  of  this  study  was  to  implement  a  generalized 

space-spatial  frequency  filter  and  study  its  operation. 


1 . 5  Coherent  Optical  Production  of  the  Hough  Transform 
The  objective  of  this  research  was  to  study  the  operation 
of  the  forward  Radon  transform  (FRT)  and  compare  it  to  the  stan¬ 
dard  Fourier  transform. 


1 . 6  Estimation  of  two  Closely  Spaced  Frequencies 

The  objective  of  this  research  was  to  use  linear  program¬ 
ming  techniques  to  improve  the  resolution  of  a  spectrum  analyzer. 

1 . 7  Restoration  of  Discrete  Fourier  Spectra 

The  objective  of  this  study  is  to  examine  the  use  of  linear 
programming  for  superresolution.  A  technique  is  presented  which' 
greatly  enhances  the  resolution  of  a  spectrum  analyzer. 

1 . 8  Study  of  Random  Access  Scheme  for  Multi-Beam  Satellite 

The  objective  of  this  study  is  to  examine  the  study  of  a 

reservation  access  scheme  for  local  area  networks.  We  also 
examine  the  study  of  an  algorithm  for  scene  matching,  using  the 
sum  of  the  absolute  values  of  the  differences  of  image  intensity 
between  the  template  and  the  scene  as  the  measure  of  similarity. 
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ABSTRACT 

This  paper  describee  a  new  technique  to  achieve 
rapid  acquisition  in  a  direct  sequence  spread 
spectrum  communication  system.  It  is  based  in 
part  on  an  upper  bound  for  the  partial  correlation 
function  of  a  PN  sequence,  and  hence  it  is  a 
fairly  general  procedure,  applicable  to  a  vide 
variety  of  systems. 

INTRODUCTION 

The  problem  of  acquisition  of  a  spread  spectrum 
signal  consists  of  correctly  estimating  the  phase 
of  a  received  PN  sequence  of  arbitrary  length 
L  ■  2N  -  1 ,  where  N  refers  to  the  number  of  delay 
stages  in  the  generating  shift  register  or,  equiv¬ 
alently,  to  the  degree  of  the  generating  polynom¬ 
ial  which  describes  the  feedback  connections. 

The  acquisition  process  consists  of  repeatedly 
comparing  the  received  PN  sequence  g(t  -  JTC) , 
with  a  local  generated  PN  sequence,  g(t  -  1TC), 
until  the  phases  ITC  and  _JTC  are  found  to  be  equal. 
At  that  point  the  locally  generated  PN  sequence 
la  said  to  be  synchronized  to  the  received  PN 
sequence  and  the  acquisition  process  is  ended. 
Throughout  this  paper,  the  carrier  phase  and  chip 
timing  are  assumed  to  be  in  perfect  synchronism. 
The  in-phase/not  in-phase  decision  is  based  on  a 
partial  correlation  process.  As  is  well  known, 
the  partial  correlation  function  of  a  PN  sequence 
is  not  nearly  as  well  behaved  as  the  autocorrela¬ 
tion  function,  it  being  a  function  of  the  parti¬ 
cular  starting  phases  1  and  J,  the  specific  feed¬ 
back  connections,  and  the  length  of  the  correlation 
YTc,  where  y  equals  the  number  of  chips  (an 
integer)  and  Tc  equals  the  chip  duration.  To 
circumvent  all  of  these  difficulties,  various 
bounds  have  been  found  ([l]-[3])  to  the  partial 
autocorrelation  function,  and  an  approximation  to 
one  of  these  bounds  is  used  In  the  analysis 
presented  in  the  next  section* 

ANALYSIS 

Reference  [2]  describes  a  recursive  method  for 
generating  a  non-linear  sequence  which  serves  as 
an  upper  bound  to  the  partial  autocorrelation 
function  of  a  PN  sequence  (via  the  shlf t-and-add 
property).  This  bound,  unfortunotely ,  is  defined 
via  a  recursive  algorithm  as  opposed  to  an  analyt¬ 
ical  form.  It  can  be  shown  that  a  good  approxi¬ 
mation  to  this  bound  Is  given  by  Y(1  -  Y/L),  and 
In  what  follova,  we  will  use  this  approximation  in 
order  to  obtain  a  simple  analytical  result.  Using 
tha  approximate  bound,  it  follows  that  for  1  t1  J  , 


K 

]0  «(*  * 


iTjdt;  tcy/i  -i) 


(i) 


The  output  of  the  correlator  in  Fig.  1  is 


+  nw(t)J  g(t  -  lTjdt  .  (2) 


which  consists  of  a  signal  term  S  (yT  )  and  a 
noise  term  n0(fTc)  •  The  noise  term  can  easily 


be  shown  to  be 


'Gaussian  with  zero  mean  and 
‘  YTr 


variance  0^(yTc)  -  j  \Tr ,  where  £  1®  the  two-sided 
power  spectral  density  of  the  input  noise.  The 
signal  term,  Sq ( YTc )  is  given  by 


(YTc 

oK)  ’  1  ^  8  (t  -  JTj  g  (t  -  iTc)dt 

■Jo 


yT 

S  c 


if  i  -  J 
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The  synchroniso/no  synchronism  decision  is  seen 
only  to  involve  the  choice  of  which  signal  term, 
So(yTc),  is  present  at  the  correlator's  output. 

Due  to  the  presence  of  noise,  this  decision  can 
only  be  made  with  a  specified  accuracy.  If  sa 
shown  in  Fig.  2,  we  wait  for  s  time  T  -  YjT  befora 
making  a  decision,  then  we  can  define  two  voltage 
regions  such  that  if  the  correlator’s  output  falls 
in  region  I,  we  decide  synchronism  is  attained, 
and  if  we  fall  in  region  II,  we  decide  the  two  PN 
sequences  are  not  synchronized.  The  probability 
of  an  erroneous  decision,  given  we  are  not  looking 
at  the  correct  phase  position,  la  equal  to 

P.  •  rlVo(Vc)  >  Vt] 

r  /  \  .  W 

*  pK(Vc) >  ei] 

where  fj  Is  shown  in  Fig.  2  and  V  the  threshold, 
is  given  below  [see(7)].  4 
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Note  that  both  in  Fig.  2  end  from  Eq.  (A),  for  This  Implies  that  the  two  arguments  of  the  erfc 

simplicity  the  decision  boundary  has  been  located  function  must  be  equal.  Hence 

halfway  between  the  two  signals.  Therefore, 


(5) 


-2T  -***- 


a  j  c 

The  probability  of  error  of  Eq.  (A)  then  becomes 


P 

a 


(6) 


€ 
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/20o(VcJ 


or. 


€ 


°o[Vc) 


(io> 


■  i* 

Eq.  (10)  defines  the  parameter  <2  euch  that  at  any 
tine  T  «  Y?T  a  decision  can  be  attempted.  Fig.  2 
shows  qualitatively  the  decision  regions  at  any 
given  Y2  <  Yj  . 


where  E  -  Y4E  ,  E  i»  the  energy  per  chip  and 

*  *  1C  c 

and  equals  P  T  ,  and 
^  sc* 


If  at  time  T  -  Y2TC  the  corre^ator,#  output  volt¬ 
age  equals  V  ,  the  following  decisions  can  be  made 
with  a  given  probability  of  error  P^: 


A  2  V 

erfc(x)  ■  -  l  e  y  dy  . 

L  /T*  J  ^ 

Finally,  V^YjTJ  Is  given  by 

MVc)  '  *^VlTc  -  £1 


■  ^iTc[l  _a]  (7) 

In  what  follows,  we  will  show  how  the  acquisition 
process  can  be  Improved  by  shortening  the  acqui¬ 
sition  time. 


The  Decision  Process  Prior  to  “  Yj  . 

At  times  It  may  not  be  necessary  to  wait  the 
full  Yj  chips  to  make  a  decision.  The  only  re¬ 
quirement  is  that  at  any  time  a  decision  is  made, 
the  probability  of  error,  Pe,  be  kept  constant. 

If  any  attempt  at  making  a  decision  after,  say, 

Y2  chips  is  made,  then  the  threshold  voltage  must 
be  set  euch  that 


(8) 

Note  that  Eq.  (8)  la  identical  to  Eq.  (A)  with  the 
exception  of  the  pertinent  parameters  (i.c.  the 
decision  mechanism  Is  Identical).  But  If  the 
probability  of  error,  P  ,  is  to  be  kept  constant, 
then  e 


e|at  t  -  YjTc  e|at  t  -  Y,T 


2  c 


(9) 


y*rfe[^W]’I*r£c[/53^7]  ‘ 


1)  Synchronization  has  been  attained  (1  ■  J)  1 

V  >  V  -  /5F  y,T  (l  - 
o  s  a  2  c \  LI  2 

2)  Synchronization  will  not  be  attained  (1  i  J. 
if 


V  <  V 
o  n 


/2P  y,T 
8  2  c 


e 
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3)  No  decision  can  be  made;  continue  the 
correlation  If 


V  <  V  <  V 

nos 

Thus,  V^Y2)  «nd  V^Yjj  become 


W  '  '  *  2 


V  2  YlTc[Z  Yj  L  /  Yjj 


and 


(U) 


Va(Y2) 


♦  c 
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From  Fig.  2  the  probability  of  no  decision  can  be 
seen  to  equal  P(no  decision) 


-P.-PfVCs+ncVl  ( 

nd  I  n  o  o  s  J 

With  V  and  V  given  in  equations  (11)  and  (12), 
respectively,  the  probability  of  no  decision  at 
t  line  Y2Tc  becomes 
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If  Q,  is  defined  as  the  probability  of  making 
a2 

a  decision  at  T  *  y~T  ,  then 

^  c 


The  realization  that  a  decision  can  be  made 
prior  to  t  ■  witbout  compromising  the  system's 

error  performance  can  be  shown  to  decrease  the  PN 
acquisition  time.  A  detailed  analysis  of  the 
actual  Improvement  in  the  acquisition  performance 
obtained  by  the  application  of  this  technique  is 
rather  lengthly  and  is  given  in  [4]. 

NUMERICAL  RESULTS 

Figures  (3a)  to  (3f)  show  curves  of  probability 
of  decision  versus  the  chip  number  (i.e.  the  value 
of  >2)  f°r  various  sets  of  system  parameters.  The 
ratio  of  energy-per-chip-to-noise  spectral  density 
E  /n  has  a  value  of  -lOdB,  -5dB  or  OdB,  and  the 
probability  of  making  an  error  is  taken  to  be 
either  0. I  or  0.01.  The  period  of  the  PN  sequence 
is  1023  (i.e.,  N  *  10)  in  all  the  figures. 

On  each  figure  there  are  two  curves,  one 
showing  the  probability  of  making  a  decision  at  or 
before  any  time  T,  and  the  other  showing  the  con* 
ditional  probability  of  making  the  decision  at  any 
specific  time  T  (i.e.  ^ -  It  can  be  seen  from 

the  figures  that  there  15  a  high  probability  that 
a  decision  can  be  made  before  time  Yj  with  no  loss 
of  system  performance. 
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A  Comparison  of  Schemes  for  Coarse  Acquisition  of 
Frequency-Hopped  Spread-Spectrum  Signals 

CLIVE  A.  PUTMAN.  STEPHEN  S.  RAPPAPORT,  senior  member,  hu  .  ano  DONALD  L.  SCHILLING,  nu.ow,  ieee 


Abstract — Three  schemes  for  course  acquisition  of  spread-spectrum 
signals  are  compared:  stepped  serial  search,  matched  filler,  and  two- 
level.  Analytical  models  and  formulas  are  developed  which  charac¬ 
terize  performance  in  adverse  environments.  Comparisons  are  on  the 
basis  of  miss  prohahilit) .  mean  time  to  search  uncertainty  region,  and 
relative  complexity. 

1.  INTRODUCTION 

r|1HE  synchronization  process  required  for  spread-spectrum 
^systems  consists  of  two  parts- acquisition  and  tracking. 
This  paper  describes  and  compares  several  methods  for  acquisi¬ 
tion.  Three  schemes  arc  considered:  stepped  serial  search, 
matched  filter,  and  two-level  acquisition.  These  schemes  have 
been  well  described  in  the  technical  literature  [l]-[4j.  The 
intention  here  is  to  develop  analytical  models  and  mathemati¬ 
cal  formulas  which  characterize  the  performance  of  the  schemes 
in  adverse  environments,  and  to  compare  the  schemes  in  terms 
of  miss  probability,  uncertainty  region  search  time,  and 
relative  complexity. 

The  investigation  is  primarily  motivated  by  consideration 
of  ground  radio  frequency-hopping  schemes.  It  is  assumed 
that  communicators  operate  in  a  push-to-talk  mode  and  that  a 
synchronization  process  is  initiated  with  each  transmission. 

A  bnef  description  of  the  schemes  under  consideration  is 
followed  by  a  general  derivation  of  performance  measures. 
Models  for  Rician  fading  and  user  or  jammer  interference 
environments  are  characterized  and  used  to  predict  perform¬ 
ance.  Methods  used  in  obtaining  numerical  results  are  out¬ 
lined  and  the  resulting  performance  curves  are  discussed  with 
reference  to  particular  spread-spectrum  applications. 

II.  GENERAL  DESCRIPTION 

In  the  stepped  serial  acquisition  scheme  shown  diagramat- 
ically  in  Fig.  1,  the  timing  epoch  of  the  local  PN  code  is  set 
and  the  locally  generated  signal  correlated  with  the  incoming 
signal.  If,  at  the  end  of  an  examination  interval,  the  threshold 
is  not  exceeded,  the  search  control  inhibits  a  clock  pulse  to  the 
PN  code  generator  so  that  the  local  code  phase  slips  to  the 
next  cell,  n  cells  per  chip,  and  the  process  is  then  repeated. 
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l  ig.  1.  Stepped  serial  search  scheme. 


Fig.  2.  Matched  filter  scheme. 

The  system  thus  performs  an  active  correlation  with  a  received 
code  sequence. 

The  matched  filter,  a  passive  correlator,  searches  the  in¬ 
coming  code  in  real  time  The  essential  structure  is  shown  in 
Fig.  2.  A  sequence  of  M  consecutive  frequencies  is  chosen  by 
the  receiver  to  establish  a  code  start  epoch  and  the  matched 
filter  performs  a  near  optimal  noncoherent  detection  as  the 
sequence  is  received. 

The  two-level  scheme  combines  passive  and  active  correla¬ 
tion  and,  as  such,  combines  the  capability  of  searching  the 
code  in  real  time  with  integration  over  a  large  number  of 
chips.  As  shown  in  Fig.  3,  a  matched  filter  detects  a  relatively 
short  M  hop  synch  prefix  and  applies  code  start  signals  to  a 
bank  of  c  active  correlators.  Each  signal  causes  the  next  idle 
correlator  to  cycle  through  K  hops  at  the  end  of  which  any 
correlator  output  exceeding  the  second  threshold  causes  a 
synch  indication;  otherwise,  the  correlator  is  again  made 
available  to  the  common  bank. 

III.  PERFORMANCE  MEASURES 
Detection  Probabilities 

For  ground  radio  systems  operating  in  a  push-to-talk  mode 
the  prime  concern  is  to  establish  the  desired  signal  code  phase 
during  some  brief  initial  acquisition  period.  We  thus  charac¬ 
terize  the  performance  of  an  acquisition  scheme  in  terms  of 
the  probability  that,  at  the  end  of  the  examination  interval, 
the  desired  signal  code  sequence  is  not  detected.  This  is  a  func- 
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Fig.  3.  Two-leveJ  scheme. 


tion  of  the  threshold  level  setting,  and  for  the  purposes  of 
comparing  schemes  we  set  this  parameter  by  specifying  the 
probability  that  the  threshold  is  exceeded  when  the  desired 
signal  code  sequence  ts  not  present.  The  performance  meas¬ 
ure  is  then  the  miss  probability  for  a  given  false  alarm  proba 
bility. 

The  input  waveform  in  the  ith  hop  interval  can  be  charac¬ 
terized  as  a  narrow  band  waveform 

Suit)  cos [*><(/) +  0<]  T*f<0  (l) 

where  wt  is  determined  by  a  known  pseudorandom  sequence, 
and  £,{/)  is  a  Gaussian  process  with  variance  a2.  This  signal 
can  be  expressed  as 

sr{0  *  r<(t)  cos  (w,<r)  +  6t  +  0!<f)|  (2) 

where  r,{t)  is  a  Rician  distributed  process. 

In  the  case  of  an  active  correlator  the  local  waveform  is 
hopped  at  the  same  rate  as  the  input  signal  and  the  detector 
output  is  integrated  over  M  hop  periods.  In  tire  case  of  a  pas¬ 
sive  correlator  the  outputs  from  M  similar  mixers  and  de¬ 
tectors  for  successive  frequencies  are  stored  and  added.  The 
system  thus  forms  the  sum  of  squares  of  independent  Rician 
variates  and  the  output  Z  can  be  shown  1 5 J  to  be  noncentral 
chi-squared  distributed  with  2M  degrees  of  freedom.  For  con¬ 
venience  we  define  a  normalized  variable  Z0,  such  that  Z0  = 
Zjo1 .  Its  first-order  probability  density  is  then  given  by 

where  the  noncentral  parameter 

<5  =  2  A? /a1 

i=i 

and  IM  .  |  is  the  modified  Bessel  function  of  the  first  kind, 
order  M-  I. 

If  the  output  exceeds  a  threshold  level  Vr  a  synch  detec¬ 
tion  is  declared.  The  detection  probability  is  then 

PD=P{Z>VT}  =  j  P20(z0)dz0  (4) 


where 

Qs,(a.b)=  (  x(x/af-le-<tUxl^2lM.l{ax)dx. 

Its  complement  is  the  miss  probability  expressed  as 

I'm  =  l  -Pi>=Qmc(s/S,VI).  (6) 

The  false  alarm  probability  is  the  probability  that  the 
threshold  is  exceeded  when  the  correct  code  sequence  is  not 
present  In  the  case  of  the  benign  environment  (noise  only- 
no  fading  or  mterferers)  considered  thus  far,  the  false  alarm 
probability  is 

Pva=Qm(  O.VT).  (7) 

Search  Time 

it  is  to  be  emphasized  that  the  acquisition  scheme  is  re¬ 
quired  to  reliably  detect  the  first  occurrence  of  the  correct 
code  epoch  during  the  search  process.  We  do  not  consider 
systems  which  can  tolerate  one  or  more  misses  of  the  code 
epoch  and  continue  to  seek  acquisition.  Let  the  uncertainty 
region  length  ;Vr  be  the  maximum  expected  delay  between  the 
locally  generated  and  desired  signal  code  phases  measured  in 
hop  intervals.  Thus,  a  second  performance  measure  which 
must  be  considered  is  the  time  taken  to  search  the  uncertainty 
region,  Ts. 

IV.  ENVIRONMENT  MODELING 
Fading  Channel  Model 

We  consider  a  Rician  fading  channel  so  that  the  received 
signal  with  which  synchronization  is  desired  is  given  during  the 
ith  hop  interval  by 

SR(t)  =  \/7a~Ps  cos  (w,l  -  8,) 

+  2  Cj\/2PS  cos  (wtt  —  ()•)  +  nw{t)  (9) 

/ 

consisting  of  specular,  scatter,  and  noise  components,  respec¬ 
tively,  The  scatter  and  noise  components  are  Gaussian  and 
independent  and  the  scatter  component  is  signal  dependent. 
Equation  (9)  defines  the  fading  model,  and  for  at  =  0,  is 
identical  to  (hat  given  in  (7| . 

For  comparative  analysis  we  equate  the  sum  of  powers  in 
the  specular  and  scatter  components  to  a  constant  Ps ,  repre¬ 
senting  the  average  received  signal  power  over  many  hops.  In 
particular  we  take 

a  +  b=  I  (10) 

where 


where  L  =  Ly/o2.  This  can  be  expressed  in  terms  of  the  gen 
eralized  Marcum  0-function  [6]  as 


o~-  2  <*t  and  6=2  ci1- 
M  |  i 


(51  Here  a  is  summed  over  many  hops  and  b  is  summed  over  many 


Pd=Qm(VS'VI) 
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multipaths.  The  quantities  a  and  b  are  indicative  of  the  relative 
strength  of  the  specular  and  scatter  components,  respectively. 
It  is  assumed  that  these  parameters  are  fixed  over  tire  acquisi¬ 
tion  period,  although  they  may  have  long  term  slow  variations. 

Following  a  derivation  similar  to  that  in  [7],  the  above 
signal  can  be  expressed  in  the  form  of  (2).  The  Rician  dis¬ 
tributed  envelope  at  the  input  to  the  detector  now  has  param¬ 
eters 

A{^y/2aP^ 

and 

a1  =  N0Bip  +  bP,.  (11) 

where  JV0  is  the  one-sided  noise  power  spectral  density  and 
B[f  is  the  detector  bandpass  filter  bandwidth. 

Multiple  Access  Model 

We  now  consider  the  interference  generated  by  other 
users  who  are  hopping  over  the  set  of  F  frequency  slots.  The 
probability  that  no  other  user  causes  interference  during  a 
given  hop  when  N  other  users  are  active  is  ( 1  —  \/F)N .  The 
probability  that  at  least  one  other  user  causes  interference  is 

p=  1  -(1  -  UFf.  (12) 

Since  we  are  integrating  over  M  hops,  the  probability  that  j 
slots  have  at  least  one  user  present  is 

=  (13) 

Jamming  Model 

In  application,  the  receiver  could  provide  for  limiting  of 
received  signal  strength  to  slightly  above  that  of  the  desired 
signal.  This  mitigates  the  effects  of  CW  jamming.  An  optimized 
jammer  would  be  forced  to  jam  a  random  selection  of  fre¬ 
quencies.  When  the  number  of  frequencies  available  F  is 
large,  the  above  model  is  reasonable  for  random  hop  or  comb 
jamming  if  iV  is  taken  as  the  number  of  interferes  or  inter¬ 
fering  tones,  respectively. 

The  Combined  Model 

To  quantify  performance  for  a  fading  environment  with 
interferes  present  it  is  analytically  convenient  to  assume  that 
the  power  of  interfering  signals  lies  in  the  specular  compo¬ 
nents,  and  that  hopped  interfering  signals  are  (hop)  synchro¬ 
nous  with  the  desired  signal.  In  addition,  all  interferes  are 
assumed  to  be  of  equal  strength  as  “heard"  at  the  receiver 
Interferes  are  assumed  to  be  statistically  independent  RF 
sources  with  random  phase.  Our  approximate  performance 
analysis  proceeds  by  noting  that  the  sum  of  interfering  and 
desired  sinusoids  produces  a  resultant  sinusoidal  signal  at  the 
receiver.  Because  of  the  assumptions  above,  the  average  power 
in  the  resultant  is  the  sum  of  the  powers  in  the  components, 
and  we  proceed  by  using  this  effective  sinusoid  as  the  resultant 
specular  component  at  the  receiver. 


The  false  alarm  probability  given/  slots  are  jammed  is  then 
P  (false  alarm)/)  =  V^)  04) 

with 

•$/  =  /2JF =f'2JFc/N0 

where  Ec  is  the  received  energy  per  hop,  J  is  the  jammer 
power  (or  limited  signal  power)  to  desired  signal  power  ratio, 
and  we  assume  that  BlF  -  1  !Th,  Th  being  the  hop  interval. 
The  overall  false  alarm  probability  is  then  determined  by  multi- 
pl>  ing  by  the  probability  of  /  slots  being  jammed  and  summing 
over  /, 

M 

^FA  =2  bU'lUpYiMis/Sj.sfL).  (IS) 

/- 0 

Similarly,  the  miss  probability  is 
M 

I'm  =  £  b(i,  M,p)Qmc(s/S^.  y/C)  (16) 

/-0 

with 

[j(J  +  a)  +  {M-j)a\2P  [jJ  +  Mi]2Ec/N0 
‘  J  ~  N0 B,F  +  bP,  ~  1  +  bEJN0 

and 

L'  =  L/(l  +  bFc/N0y 

V.  ACQUISITION  SCHEME  PERFORMANCE 
Stepped  Serial  Search 

There  are  many  ways  to  control  the  search  process  (8J.  In 
the  scheme  considered  here  it  is  assumed  that  three  successive 
threshold  exceedences,  or  hits,  are  required  to  initiate  the 
tiacking  phase.  Any  one  miss  causes  the  cell  to  be  rejected, 
lienee,  the  overall  probability  of  false  alarm  is 

/’fa  =  Pfh3  ov 

where  Pyih  the  probability  of  a  false  hit,  is  given  by  (15),  and 
the  overall  miss  probability  is 

/’u  =  l-V  (18) 

where  Ph  is  the  probability  of  a  hit. 

The  magnitude  of  the  detected  energy  is  dependent  on  the 
relative  delay  T  between  the  locally  generated  and  desired 
signal  code  sequences.  Tins  can  be  seen  from  the  correlation 
diagram  in  Fig.  4.  Due  to  the  stepping  nature  of  this  scheme, 
the  worst  case  correlation  will  occur  for  a  delay  of  ±Thl2n. 
Assuming  equally  likely  delays,  the  average  loss  from  full 
correlation  can  be  found.  In  the  calculation  of  the  parameter 
.Yj,  the  desired  signal  energy  must  be  multiplied  by  the  factor 
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It  should  be  noted  that  this  scaling  is  used  only  for  the  stepped 
serial  search  scheme.  This  is  because,  in  practical  implementa¬ 
tions  of  both  the  matched  filter  and  two  level  schemes,  the 
largest  output  of  the  passive  correlator  (in  the  interval  where 
that  output  exceeds  the  threshold)  would  be  used  to  deter¬ 
mine  the  coarse  code  epoch,  so  there  will  be  essentially  no 
correlation  loss. 

Since  each  cell  is  examined  for  a  period  MTh  the  search 
time  for  no  false  hits  would  be 

Ts  =  {nNc  +  2)MTh.  (20) 

False  hits  tend  to  increase  this  time,  but  since  PF H  is  small 
the  effect  will  not  be  significant. 

Matched  Filter  Acquisition 

The  analysis  is  straightforward.  False  alarm  and  miss  prob¬ 
abilities  are  given  by  (15)  and  (16).  Since  the  search  time  is 
just  the  time  taken  to  search  the  uncertainty  region  once, 

T,  =  {Nc+M)Th‘  (21) 

Two-Level  Acquisition 

This  scheme  is  analyzed  in  [4] .  A  false  alarm  occurs  when 
the  matched  filter  gives  a  false  start  signal  which  finds  an  idle 
correlator,  which  in  ti'  i  generates  the  false  synch  indication. 
The  probability  of  this  occurring  is 

^FA  -/Vai^FAiO  *.(22) 

where  /Vai  a°d  ^FA2  are  the  fake  alarms  conditions  for  the 
passive  and  active  correlators,  respectively,  and  B(c,g)  is  the 
blocking  probability  for  the  queue  ofc  active  correlators  given 
by  the  Erlang  loss  formula  with#  =  KPyA  , . 

A  correct  code  acquisition  epoch  will  be  missed  either  when 
it  is  missed  by  the  matched  filter,  or  when  it  is  correctly  de¬ 
tected  by  the  matched  filter  but  is  either  missed  by  an  active 
correlator  or  finds  all  active  correlators  engaged.  The  proba¬ 
bility  of  this  occurring  is 

Pm  =  Pm i  +  V  ~ PMMe.g)  +  V  -B(c,g))PMJ)  (23) 

where  Pm  t  and  Fvf  2  are  the  miss  probabilities  for  the  passive 
and  active  correlators,  respectively. 


The  detection  and  false  alarm  probabilities  for  the  sub¬ 
systems  are  provided  by  (15)  and  (16)  with  K  substituted  for 
M  in  the  case  of  the  active  correlators 

The  search  time  is  given  by 

Ts  =  (Nc  +  M  +  K)Th.  (24) 

VI.  RESULTS 

Starting  with  the  serial  search  scheme,  the  basic  perform¬ 
ance  curve  is  presented  in  Fig.  5  as  the  miss  probability  versus 
energy  per  chip  to  noise  density  ratio  for  various  integration 
periods  M  (hops)  in  a  benign  environment.  It  is  seen,  for  ex¬ 
ample,  that  a  miss  probability  of  about  10“  3  is  obtainable  for, 
say,  L'd A0  =  5  dB  with  M  of  20.  To  take  full  advantage  of 
this  scheme,  performance  can  be  improved  by  increasing  M 
without  increasing  hardware  requirements  but  increasing 
search  time.  However,  the  accompanying  set  of  results  for 
an  environment  where  half  the  received  energy  is  in  the  scat¬ 
ter  component  and  there  are  five  interferes  present  shows  that 
longer  integration  times  are  in  fact  necessary  to  produce  miss 
probabilities  of  the  same  order. 

The  performance  curves  for  a  matched  filter  scheme  are 
provided  in  Fig.  6.  If  the  matched  filter  complexity  M  is 
equated  to  the  stepped  serial  search  integration  time  in  hops, 
we  find  that  the  serial  search  scheme  is  slightly  superior  in 
performance.  This  is  due  to  the  fact  that  the  serial  search 
scheme  considered  requires  three  successive  hits  for  a  synch 
indication,  whereas  the  matched  filter  relies  on  a  single  thres¬ 
hold  exceedence. 

With  the  two-level  scheme  it  is  meaningful  to  seek  those 
thresholds  L\  and  L2  which  minimize  the  miss  probability, 
while  the  false  alarm  probability  is  held  fixed.  For  each  given 
parameter  set  consisting  of  A/,  K ,  c,  £C/Ar0,  relative  strength 
of  scatter  component  £>,  and  number  of  users  N,  the  speci¬ 
fied  false  alarm  probability  establishes  L2  for  a  fixed  L\. 
Thus,  each  computation  of  miss  probability  is  minimized  over 
L\  subject  to  the  constraint  of  a  given  false  alarm  probability. 
Curves  for  the  two-level  scheme  are  presented  in  Fig.  7.  Here 
schemes  with  one  and  three  active  correlators  are  compared. 
The  active  correlation  period  is  10  hops.  It  can  be  seen  that 
an  improvement  in  performance  of  1-2  dB  in  the  region  of 
interest  can  be  gained  by  using  a  bank  of  three  active  correla¬ 
tors  as  opposed  to  one.  Comparison  with  ihe  matched  filter 
results  show  that,  for  similar  performance,  the  two-level 
scheme  with  K  -  10  requires  a  passive  filter  complexity  only 
half  of  the  matched  filter  on  its  own. 

VII.  CONCLUSIONS 

The  serial  search  scheme  has  been  shown  to  provide  good 
detection  capabilities  even  in  adverse  environments  and  is 
relatively  simple  to  implement.  The  cost  is  a  long  uncertainty 
region  search  time  compared  to  that  for  the  matched  filter  or 
two-level  scheme.  If  this  can  be  made  small,  such  as  in  systems 
with  reasonably  fast  hop  rates  and  short  codes,  or  if  some 
rough  code  epoch  can  be  maintained  with  clocks,  then  the  ac¬ 
quisition  times  obtained  can  be  acceptable,  and  the  system  Is 
preferable  in  terms  of  detection  reliability. 

Matched  filter  detection  requires  a  complex  hardware 
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structure  for  similar  performance.  However,  the  matched 
filter  searches  the  uncertainty  region  in  real  time,  offering 
a  clear  advantage  where  it  performs  adequately.  It  is 
therefore  preferable  in  systems  operating  in  less  adverse 
environments  where  simpler  filters  will  provide  the  desired 
performance. 

The  two-level  scheme  combines  the  advantages  of  the 
rapid  search  time  of  the  matched  filter  with  the  detection 
reliability  of  an  active  correlator.  It  would  therefore  find 
application  in  systems  which  do  not  rely  on  a  time  reference 
and  which  utilize  long  codes  but  require  rapid  and  reliable 
acquisition  in  adverse  environments.  In  situations  with  im¬ 
paired  signal  detectability  due  to  interference  and/or  scat¬ 
ter,  the  parallel  bank  of  active  correlators  in  this  scheme 
provides  additional  versatility.  Roughly,  even  while  some  cor¬ 
relators  which  may  have  been  falsely  engaged  are  operating, 
successful  signal  acquisition  is  still  possible  with  an  idle  cor¬ 
relator.  This  parallelism  is  not  available  in  the  other  schemes 
considered  here.  These  advantages  are  at  the  expense  of  some¬ 


what  higher  complexity,  which,  however,  is  not  necessarily 
formidable. 
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ABSTRACT 

Tracking  requirements  for  a  frequency  hopped 
spread  spectrum  system  with  a  view  to  application 
in  mobile  radio  are  investigated.  Models  charac¬ 
terizing  Rician  fading  and  multiple  access  or  jam¬ 
ming  environments  are  used  to  formulate  signed, 
detection  probabilities.  These  determine  the  mean 
ti*e  to  loss  of  lock.  Optimized  tracking  in  these 
environments  is  then  considered. 


.1.  INTRODUCTION 

The  frequency  hopping  (FH)  waveform  as  a  pseudo¬ 
random  modulation  technique  has  the  attraction  of 
ease  of  implementation  and  inherent  frequency  div¬ 
ersity  C X ] .  A  central  feature  of  a  FH  system  is 
the  pseudorandom  code  generators  at  both  the  trans¬ 
mitter  and  receiver,  capable  of  producing  identical 
codes  with  proper  synchronization.  The  pseudoran¬ 
dom  code  sequence  is  used  to  switch  the  carrier 
frequency  via  a  frequency  synthesizer  and  wideband 
mixer.  When  the  synthesizer  in  the  receiver  is 
switched  with  the  synchronized  sequence  the  fre¬ 
quency  hops  on  the  received  signal  will  be  removed, 
leaving  the  original  modulated  signal. 

The  synchronizat ion  of  FH  spread  spectrum  signals 
can  be  broadly  divided  Into  two  phases:  acquisition 
and  tracking.  A  typical  system  would  first  obtain 
coarse  acquisition  of  the  received  pseudo-noise 
encoded  signal  [2],  followed  by  fine  acquisition 
by  the  tracking  loop.  At  the  end  of  the  acquisi¬ 
tion  period  loop  parameters  would  be  adjusted  from 
acquisition  to  tracking  settings  and  the  system  is 
said  to  be  in  lock.  There  are  two  possibilities 
at  this  point.  Depending  on  the  scheme  used, 
either  one  or  a  series  of  false  alarms  may  have 
causea  the  system  to  incorrectly  enter  the  lock 
■tate,  or  the  correct  code  epoch  may  have  been 
detected  causing  the  system  to  correctly  enter  the 
lock  state.  The  requirements  for  the  tracking 
phase  are  to  maintain  tracking  of  the  correct  code 
phase  In  the  latter  case,  yet  to  quickly  reject 
the  false  lock  condition  to  allow  resumed  search 
for  the  correct  code  phase  in  the  latter  case. 
Although  this  paper  focuses  primarily  on  tracking 
.  of  FH  signals,  much  of  the  analysis  applies  direct¬ 
ly  or  with  little  modification  to  other  spread 
spectrum  schemes  such  as  direct  sequence  or  hybrid 
systems.  In-lock  detection  reliability  la  charac¬ 
terized  by  mean  time  to  loss  of  lock  when  mean 
time  to  reject  false  lock  Is  given. 

2.  THE  TRACKING  LOOP 

The  coarse  acquisition  schemes  described  in  (?I 


will  provide  acquisition  of  the  received  pseudo¬ 
random  modulating  sequence  to  within  about  half  a 
code  bit,  or  chip.  Gince  finer  acquisition  is 
needed,  and  since  reliance  cn  clock  stability  to 
maintain  accurate  synchronization  is  generally  not 
sufficient,  even  for  short  transmission  push  to 
talk  operation,  a  tracking  loop  Is  usually  employed. 

The  early-late  gate  tracking  loop  as  used  for 
range-gate  tracking  in  ra iar  systecs  [?),  lends  it¬ 
self  veil  to  the  frequency  hepped  signal  applica- 
*ion.  This  loop  ar.i  it r.  associated  waveferrs  are 
rhown  in  figs.  la,t.  7h*»  rating  wav*  f '-nr.  g't)  from 
the  VCO  is  alternately  negative  and  positive  with 
the  minus  to  plus  transitions  coinciding,  with  the 
chip  to  chip  frequency  hopping  instants  of  the 
locally  generated  waveform  G  (t).  The  error  signal 
is  generated  by  integrating  ?he  multiplied  gating 
signal  and  envelope  detector  output  v(t)  and  is 
proportional  to  the  delay  T  between  the  local  and 
Incoming  code  phases.  The  resulting  discriminator 
character! st i c  is  shsvr  in  Fig.  lc.  The  error 
voltage  is  used  to  advance  or  rttard  the  local 
waveform  into  alignment . 

Inherent  in  the  loop  operation  will  be  seme  phase 
Jitter  which  will  increase  with  decreasing  signal 
to  noise  ratio.  The  variance  of  this  phase  error 
is  given  in  (3l  for  a  single  received  pulse  as 


T  T 

*  JSLX 

8E/N 

o 


where  is  the  gate  width,  T  is  the  pulse  width 
and  E/NJ  is  the  pulse  energy  to  single-sided  noise 
power  spectral  density  ratio.  For  a  frequency  hop¬ 
ped  system  T„  is  equal  to  the  chip  duration  Tc ,  and 
assuming  that  the  channel  correlation  bandwidth  is 
greater  than  the  hop  rate,  T  is  also  equal  to  Tr. 
Hence 


8mE  /N 
c  o 


U) 


where  m  In  the  number  of  h«’pn  in  the  integration 
period  and  F.r  is  the  energy  per  chip.  This  result 
also  agrees  wit:,  that  in  { b  ]  with  the  loop  band¬ 
width  written  as  =  l/mTc. 

3.  I N- LOCK  MONITORING 

Once  the  system  enters  the  leek  state,  it  becomes 
necessary'  to  monitor  the  tracking  of  the  received 
code  phase.  This  is  achieved  using  the  detector 
of  Fig.  Ta.  The  input  waveform  in  the  i^h  hop 
interval  can  be  characterized  as  a  narrowband  wave¬ 
form  . 

r.R(t)  *  .'TFcos(w1(t-iTc)*81]  ♦  Rjlt-IT.)  (2) 
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vhere  la  determined  by  a  known  pseudorandom  se¬ 
quence.  and  g^{t)  Is  a  gaussian  process  with  vari¬ 
ance  o2.  The  envelope  of  this  waveform  is  then  n 
Rician  distributed  process,  and  since  integration 
is  over  M  hops,  the  output  Z  is  the  sum  of  squares 
of  M  independent  Rician  variates.  For  convenience 
we  define  a  normalized  variable  Z0,  such  that 
ZQ  =  Z/a^.  Its  first  order  probability  density  is 
then  given  by  [5} 


(M-l)/2  -(*  *S)/2 

PZ  (io>  =  i(VS)  e  Vl 

o  2 


where  1^  is  the  Modified  Bessel  function  of  the 
first  kind,  order  M-l ,  and  the  noncentrality  para¬ 
meter  is 


S  =  2MP  jo2 . 
3 


(3) 


If  after  MTc  seconds  the  output  exceeds  a  thresh¬ 
old  level  V-ji  a  synch  detection  is  declared.  The 
detection  probability  is  then 

pd  *  p(z>vT)  =  Tpz  U0)d«  (M 

L  o 

vhere  L  =  V^/o^. 

This  can  be  expressed  in  terms  of  the  generalized 
Marcum  Q-function  as 

PD  =  (y.S.v'L)  (5) 

where  Q^(a,b)  =  /  x(x/a)^_*e”^a "+x  ^^1^  ^(ax).ix. 
b 

Its  complement  is  the  miss  probability  expressed 
as 

PM  =  1  '  PD  =  i6) 

The  false  alarm  probability  for  this  system  is 
the  probability  that  the  threshold  is  exceeded 
/hen  the  correct  code  sequence  is  not  present.  In 
the  case  of  the  benign  environment  considered  thus 
far,  the  false  alarm  probability  is 

PFA  *  V°.^  •  (7) 


The  tracking  loop  phase  error  will  reduce  the 
correlation  between  received  and  locally  generated 
waveforms,  thereby  degrading  in-lock  detector  per¬ 
formance,  and  must  be  included  in  the  detection 
probability  formulations.  This  Is  done  by  assum¬ 
ing  that  the  phase  error  is  approximately  Gaussian 
distributed  (3]  so  that  90  percent  of  th^.tiroe  the 
phase  error  is  less  than  3 oT.  It  can  then  be  seen 
from  the  diagram  of  Fig.  2b  that  the  correlation 
1b  reduced  by  a  factor  less  than  3o<j./Tc  and  honce 
the  signal  power  at  the  detector  output  is  scaled 
by  a  factor  greater  than 

d  -  (1  -  3aT/Tc)2 


Hence  the  miss  probability  is  given  by  (6)  with, 
from  ( 3) , 


S 


2MdP 
_ 8_ 

Vo 


(0) 


where  B^p  is  the  system  noise  bandwidth. 


Combining  the  tracking  loop  with  a  stepped  serial 
search  acquisition  loop  results  in  the  circuit  of 
Fig.  2c,  with  each  loop  providing  fine  and  coarse 
VCO  control  respectively.  Separate  envelope  detec¬ 
tors  are  used  for  tracking  and  signal  detection 
since  it  is  likely  that  their  characteristics  will 
differ.  It  can  be  noted  that  the  in-lock  detector 
is  inherent  in  the  same  circuit  with  probable  para¬ 
meter  changes.  However,  if  matched  filter  acquisi¬ 
tion  were  utilized,  such  circuitry  would  be  addi¬ 
tional  . 


5.  CONTROL  STRATEGY 

Information  received  from  the  threshold  detector, 
that  is  whether  or  not  the  threshold  is  being  ex¬ 
ceeded,  is  interpreted  by  a  control  system  which 
decides  whether  or  not  to  continue  tracking  or  re¬ 
turn  to  the  search  phase.  This  decision  will  be 
based  on  some  form  of  control  strategy  [6J,  In  its 
simplest  form,  a  single  occurrence  of  the  threshold 
not  being  exceeded,  a  miss,  causes  the  return  to 
the  search  phase.  Alternatively  a  series  of  such 
occurrences  may  be  required  before  this  action 
takes  place  as  shown  in  Figs.  3a, b.  Here  n  is  the 
number  of  additional  in-lock  states.  The  differ¬ 
ence  between  the  two  strategies  shown  is  that  a 
threshold  exceedence,  or  hit,  either  causes  control 
to  return  to  the  first  in-lock  state  or  the  previ¬ 
ous  one . 

The  control  strategy  can  be  represented  by  a 
finite  Markov  chain  [7).  It  is  assumed  that  the 
probabilities  of  a  hit  or  a  miss  are  constants. 

The  Markov  chains  for  the  strategies  of  Figs.  3a, b 
are  3hovn  in  Figs.  ba,b  respectively,  where  p  is 
the  probability  of  a  hit  and  q  is  the  miss  proba¬ 
bility.  The  first  state  has  a  reflecting  wadi,  and 
the  last  state  is  an  absorbing  state  representing 
rejection  of  the  current  code  phase.  The  Markov 
chain  can  be  described  in  terms  of  it9  transition 
matrix  P  whose  element  is  the  probability  of 
transition  from  state  1  to  state  J.  For  the"  Markov 
chain  of  Fig.  Ua,  the  (n+2)x{n +2)  transition  matrix 
is 


p  q  0 . 0 

P  0  q 

0  p  0 

P  = 
a 

0  q  0 

:  P  0  q 

0 . 0  0  1 

and  for  that  of  Fig.  Ub, 

p  q  0 . 0 

p  0  q  : 

p  0  0  : 


p 


b 


:  *0  q  0 

t>  0  0  q 

0 . 0  0  1 


respectively.  In  each  case  the  matrix  can  be  parti¬ 
tioned  as 


P 


Q 

0 


R 

1 


so  that  the  (n+1 )x<n+l )  matrix  Q  describes  the 
Markov  chain  without  the  absorbing  state,  R  Is  a 
column  vector  of  n  zeros  and  a  q,  and  0  is  a  vector 
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of  n*l  zeros. 

The  mean  transit icn  time  from  any  state  i  to  ab¬ 
sorbing  state  r.*l  is  given  ty  the  element  of 
the  vector  [?  ] 

t  *  li  -  Qf  :t  *  :;_1t 

where  I  is  an  identity  matrix  and  t  Is  a  vector 
whose  element  a  t^  are  the  times  spent  in  state  i. 
The  mean  time  to  loss  of  lock.  Tp,  is  thus  given  by 
the  first  el  omen  t  cf  T,  the  mean  t  ran  a  i  t  i  on  t  i  me 
from  state  0  to  state  r.*l .  Using  elementary'  ma¬ 
trix  algebra, 


u  # 

T.  ■  I  i ' V  i 

L  1=0  det  N 


where  is  the  cofactcr  of  j .  If  t( ,  the 
threshold  detector  Integration  times,  are  a  con¬ 
stant  KTC,  then 


n 

T  =  MT  *:  Li0  . 
c  1=0  det  K 


(9) 


This  formula  upq  lies  to  any  control  strategy  ar.d 
is  easily  calculated  when  matrix  algebra  routines 
are  available. 

A  slightly  different  approach  yields  closed  form 
expressions  for  a  giver,  control  strategy  of  n 
stages.  The  state  transition  diagrams  can  be 
transformed  into  signal  flow  graphs  IS].  The 
strategies  of  Figs.  3a, b  can  then  also  be  repre¬ 
sented  by  the  signal  flew  graphs  of  Figs.  5a,t 
where  the  variable  t  ru presents  the  unit  path  leng¬ 
th,  in  this  case  the  fixed  integration  time.  The 
neon  path  length  from  rode  C  to  node  n*l  is  then 
the  mean  time  to  loss  of  lock.  Since  the  mean  path 
length  is  Just  the  sum  of  all  path  lengths  times 
their  probability  It  can  te  obtained  from  the 
graph  transfer  function  Gn(t)  by  di f ferentlating 
with  respect  to  t  and  setting  t=l  (3]: 

T *  =  MT  djn  U)l  (10) 

dt  | t=l 

The  graph  transfer  function  for  Fig.  5b  is  ob¬ 
tained  by  inspection  using  Mason's  formula  (9l  a.3 


Gn(t) 


n+T  n+i 
S>  t _ 

i  -  pt  :  qh1 
1=0 


Applying  equation  (10)  and  u3ing  p  =  1-q  it  is 
found  that 


MT  1- 


n  +  1 


i ,  \  n+l 

(l-q)q 

Hate  that  as  the  miss  probability  approaches  unity 
ve  find 


llmTL  =  M(nH)  (n) 

q-*l 

This  reoult  will  be  used  as  a  basis  for  comparison 
of  various  control  strategies .  A  general  closed 
fora  expression  for  7^  for  the  graph  of  Fig.  5a 
has  not  been  found,  but  the  same  technique  yields 
An  expression  for  each  value  of  n.  In  each  case 


equ&ti.n  (11)  applies. 

,  ^  ir*  c'r-t‘r  -  determine  the  roar.  time  to  lo;s 
cl  It:-  given  that  the  correct  cote  epcch  has  been 
detect *•:,  cr  the  mean  hold-in  time  TH,  ve  suteti- 
tete  1m  for  q  in  the  exrresuior.  ci.taVr.ed  for  It 
for  tr.e  control  strategy  need.  Similarly,  insti¬ 
tuting  1  - 1  for  q  yields  the  mean  time  to  loos  cf 

lock  given  incorrect  acquisition,  cr  mean  time  for 
false  1 :ck  Tv* . 


Fad  in?  Channel  M.  ■  •  •  1 

-’e  cc  :is  i  do  r  a  Kiel  an  fa  I  i  r.g 
received  signal  with  which  ayn 
sired,  is  gii;en  during  the  ith 

Sa( t  ,  =  r2a.  P  cos (w, t-0  .  )+Ic  , 
n  is  i  j  J 


channel  sc  that  the 
hror.izuticn  i3  ce- 
■op  interval  by 

^coS(w.t-^)+rlu(0 

(11) 


consisting  of  specular,  scatter  ar.d  r.cise  compo¬ 
nents  respectively.  The  scatter  and  noise  compo¬ 
nents  are  Gaussian  ar.d  ir.dej  enier.t  and  the  scatter 
component  is  signal  dependent,  equation  (11;  de¬ 
fines  *  he  fading  model  and  for  a, =0  is  identical  to 
that  given  in  1 10;.  1 

Per  comparative  analysis  wo  equate  the  sum  cf 
powers  in  the  op ecular  and  scatter  components  to  a 
constant,  ?s ,  representing  the  average  received 
signal  power  over  many  hops.  In  particular  we 
take 


a  +  b  =  1 

where  a  =  1  I  a.  and  b 
M  i 


Here  a  is  summed  over  many  h  ps  and  b  i3  summed 
over  many  nultipatha.  The  quur.titles  a  a:.  ;  b  are 
indicative  cf  the  relative  strength  of  the  specular 
and  scatter  components  respectively.  It  is  assumed 
t.nat  these  parameters  are  fixed  over  the  acquisi¬ 
tion  period,  although  they  may  have  lcr.g  term  slow 
variations . 

tell  owing  a  derivation  uiml I  nr  to  that  in  [ll], 
th'Q  at^ve  3iFr*al  ~a:-  be  expressed  in  the  form  of 
(2).  The  noneent rul i ♦ y  parameter  C  in  the  equation 
for  tl.e  tiss  prot&MUty  lo)  ther.  beccr.es 


s  =  zylipjk  (13) 

2 

where  e*-  =  N  +  b? 

o  IF  9 

Multiple  ftcrers  rc  ?■•! 

now  consider  the  interference  generated  by 
other  users  who  are  hcppir.g  over  the  set  of  F  freq¬ 
uency  slots.  The  probability  that  no  other  user 
causes  interference  during  a  given  hop  when  S  other 
users  ure  active  in  (l-l/K)x.  The  probability  that 
at  least  one  other  user  causes  interference  is 


p  =  i  -  u-i/k:s. 


(1U) 


Since  we  are  integrating  over  M  hops,  the  probabil¬ 
ity  that  J  slot 3  have  at  least  one  user  present,  is 


b(j,M,p) 


PJU-p)M-J 


(15) 


Jamming  Model 

In  application  the  receiver  could  provide  for 
1  inviting  of  received  signal  strength  above  that  of 
the  desired  signal.  This  mitigates  the  effects  of 
CV  Jamming.  An  optimized  Jammer  would  be  forced 
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1 


to  Jam  a  random  selection  of  frequencies.  When  the 
number  of  frequencies  available  F  h  large,  the 
above  model  is  reasonable  for  random  hop  or  comb 
Jamming,  if  N  Is  token  as  the  number  of  interfer- 
ers  or  interfering  tones  respectively. 

The  Combined  Model 

To  quantify  performance  for  a  fading  environment 
with  interferers  present  it  Is  analytically  conven¬ 
ient  to  assume  that  the  power  of  interfering  sig¬ 
nals  lies  in  the  specular  components  and  hopped 
interfering  signals  axe  assumed  to  be  synchronous 
with  the  desired  signal.  In  addition,  all  inter¬ 
ferers  are  assumed  to  be  of  equal  strength,  and 
where  interferera  and  desired  signal  are  present 
in  the  same  slot  their  powers  are  assumed  to  be 
additive. 

The  false  alarm  probability  given  J  slots  are 
Jammed  is  then 

P(false  alarm  |  J)  -  (l6) 


with 

Sj  “  J2JVNoBIF  =  j2JEc/No 

where  J  is  the  Jammer  power  (or  limited  signal 
power)  to  desired  signal  power  ratio  and  we  assume 
that  Bjp^l/T  .  The  overall  false  alarm  probability 
is  then  determined  by  multiplying  by  the  probabil¬ 
ity  of  J  slots  being  Jammed  and  summing  over  J , 


P 


FA 


M 

Z  b(J,M,p)a,(*'s7,  /L). 
J=0  J 


Similarly  the  miss  probability  is 


(17) 


(18) 


with 

S, 


[j(J+a)+(M-j)a]2Ps 

N  Bt„  +  bP 
o  IF  d 


ljJ+Mal2Ec/No 

1+bVNo 


and 

L’  =  L/(l  +  bE  /N  ). 

c  o 


6.  RESULTS 

The  performance  parameter  of  interest  is  the 
mean  hold-in  time  T„  given  the  mean  time  for  false 
lock  Tp. .  For  graphical  convenience  we  generate 
the  logarithm  of  these  times  normalized  to  the  hop 
period  T  .  Thus  given  T^*  M,  m,  and  n,  Pp.is 
specified  so  that  for  any  set  of  parameters^Br/No, 
relative  strength  of  scatter  component  b,  and  num¬ 
ber  of  interferers  N,  the  detector  threshold  level 
and  hence  P  and  finally  TH  can  be  calculated. 

Consider  a  tracking  system  employing  the  strat¬ 
egy  of  Fig.  5a.  The  basic  performance  curve  is 
presented  in  Fig.  6  as  the  log  of  normalized  mean 
hold-in  time  vs  energy  per  chip  to  noise  density- 
ratio  for  varying  number  of  control  strategy 
states  n.  For  the  loop  parameters  M  and  a  chosen, 
the  need  for  a  control  strategy  with  one  or  more 
additional  states  ia  evident.  For  example,  if  the 
mean  hold  in  time  is  required  to  be  greater  than 
an  hour  at  an  Ec/N  of  6dB  with  a  hop  rate  of 
lOKHz  and  a  mean  time  to  reject  false  lock  of  lOmS, 

then  a  value  of  zero  for  n  (l.e.,  a  single  miss 
causes  immediate  return  to  the  search  phase/  is 


sufficient  in  a  benign  environment.  However,  when 
the  environment  consists  of  a  fading  channel  where 
half  the  received  energy  is  in  the  scatter  compo¬ 
nent  ani  there  arc  50  interferers  present,  cne  or 
two  additional  in-lock  states  are  needed. 

Minimum  tracking  system  requirements  for  a  speci¬ 
fied  mean  hold- in  time  in  a  given  environment  are 
determined  by  the  specified  mean  time  to  reject 
false  lock.  To  see  this  consider  the  curves  of 
Fig.  7»  log  of  normalized  mean  hold-in  time  vs  log 
of  normalized  false  lock  time,  which  show  the 
trade-eff  in  performance  for  reduced  mean  false 
lock  time.  With  reference  to  the  results  for  an 
environment  where  half  the  received  energy  is  in 
the  scaler  components  and  there  are  50  interferers 
present,  we  see  that  for  the  parameters  chosen,  if 
TfL  *s  re1u^re(^  to  be  less  than  10ms  then  two  addi¬ 
tional  in-lock  states  should  be  used. 

The  curves  of  Figs.  8  and  9  show  the  effect  of 
the  relative  strength  of  scatter  component  and  num¬ 
ber  of  interferers  respectively  on  the  mean  hold-in 
time.  The  effect  of  the  fading  environment  is  seen 
to  be  more  severe  than  that  of  the  multiple-user  or 
Jamming  environment,  and  control  strategy  require¬ 
ments  can  be  determined  from  these  results. 

Similar  results  were  obtained  for  a  system  using 
the  control  strategy  of  Fig.  3b,  and  were  found  to 
have  little  significant  difference  from  the  above 
results  for  the  values  of  n  considered.  Either 
strategy  would  thus  suffice. 

Finally,  there  were  also  found  to  be  little  signi¬ 
ficant  difference  in  the  performance  obtained  from 
a  system  with  n=Q  and  M=20  compared  to  a  system 
with  n=2  and  M=10  (application  of  equation  (11) 
yields  a  normalized  mean  time  to  loss  of  lock  of 
30  in  each  case).  This  implies  that  adequate  per¬ 
formance  can  always  be  obtained  simply  by  increas¬ 
ing  the  in-lock  detection  integration  time,  provi¬ 
ded  the  constraint  Tp^/Tc  >>  M  is  met.  However, 
the  analytical  results  predict  average  performance 
and  do  not  reflect  the  disastrous  effect  of  events 
such  as  deep  fades  which  are  often  prevalent  in  a 
mobile  environment.  The  in-lock  states  and  shorter 
integration  times  would  be  preferable  since  a 
single  miss  does  not  indicate  loss  of  lock. 

7.  CONCLUSIONS 

A  system  for  the  tracking  of  frequency  hopped 
signals  has  been  evaluated  in  terms  of  mean  hold 
in  time  for  given  mean  time  to  reject  false  lock. 

It  was  found  that  tracking  system  requirements  for 
a  benign  environment  are  greatly  extended  by  the 
Rician  fading  environment,  and  only  to  a  limited 
degree  by  a  multiple  access  or  Jamming  environment. 
System  design  should  therefore  utilize  the  predict¬ 
ed  results  for  an  expected  environment. 

In  the  case  of  a  benign  environment,  performance 
of  the  tracking  system  In  terms  of  mean  hold-in 
time  cart  always  be  improved  by  increasing  the  in¬ 
lock  detector  integration  time,  providing  that  the 
mean  time  to  reject  false  lock  Is  acceptable.  How¬ 
ever,  a  search  lock  strategy  with  one  or  more  addi¬ 
tional  in-lock  states  gives  similar  performance 
for  the  same  mean  time  to  reject  false  lock  and  is 
preferable  in  the  case  of  the  fading  environment. 
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la.  farly-late  gate  tracking  loop. 


2c.  Serial  search  acquisition  and  tracking  system 
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Coherent  optical  production  of  the  Hough  transform 


George  Eichmann  and  B.  Z.  Dong 


Tht*  tbuigh  tr.uislitrni  I II  Tl  IS  .til  Htu  it- lit  -.h.tpr  »htr*  t.-r  lh.it  mu|*k  Mf<w:ht  into  a  iw..  p.ir.trm  l»  r  Ira 
til  re  >[>act*.  KiTt'iitly  it  lias  hren  poinh'tl  tail  that  thr  hrwttrii  Ration  tranUorm  i  FR  I  i.  ui.i  know  n  tnan  lh<- 
theory  of  eomputeii  tomography .  ami  (hi*  HT  are  t  ijtnvaif  nt  h< r  hmarv  miagr-  in  fin**  paper,  anal  *^  <  (»her 
out  opt u  al  implementation  of  the  h  K  I  i^  iIim  u**Mti  I  he  t  K  1  will  not  i>mx  he  ■  >t  use  in  implement t  he 
HT  shape  descriptors  hut  also  act  as  «  coherent  opt.mJ  prepnaessor  for  the  impleme -ruinm  of  multidi 
mensional  convolution,  correlation,  and  spei  tral  ana!\sis  inii^  1  i )  mousi(M*ptit  al  signal  pr*K  essin^  devn  es 
Several  differ(*nt  cidierent  oplnal  F  It  i  '  ,irt  lotei  fur*  >  are  |>n*a  n!#*<i  Fxprrni. m!aJ  results  using  (..n\m 
tional  coherent  Fourier  transform  configuration  are  civen  The  relationship  between  the  mlu  rent  optual 
implementation  of  the  FRT  and  the  inverse  Ration  transform,  an  important  tool  in  computed  toniographv. 
is  also  detailed. 


I.  Introduction 

The  Hough  transform1  for  detecting  straight  lines  has 
been  discussed  in  the  context  of  image  pattern  analy¬ 
sis.-  4  The  Hough  transformation  maps  a  point  in  the 
a  y  plane  into  a  sine  curve  in  the  {pM)  transform  plane 
using  the  formula  /)  =  t  cosf)  4-  y  sinfT  A  straight  line 
in  the,t-y  plane  is  mapped  into  a  point  in  the  transform 
plane.  Kverv  picture  edge  clement,  which  may  he 
considered  as  a  unit  length  line  segment,  is  mapped  into 
a  point  in  the  transform  plane.  All  collinear  edge  ele¬ 
ments  are  mapped  into  the  same  point  in  the  transform 
plane.  The  total  number  of  unit  length  collinear  line 
segments  form  a  histogram  at  a  point  in  the  transform 
plane.  A  single  long  or  many  short  collinear  line  seg¬ 
ments  can  yield  the  same  histogram  value.  Images  that 
contain  edge  elements  only,  the  so-called  shape  func¬ 
tion,  can  always  he  decomposed  as  a  set  of  approximate 
line  segments.  The  Hough  transform  maps  these  line 
segments  into  a  smaller  feature  space  to  form  the  1  lough 
transform  shape  descriptors. 

Recently  it  has  been  pointed  out  by  l)eansr’  that  the 
Hough  transform  for  a  binary  image  is  equivalent  to  the 
forward  Radon  transform1'  (FRT)  of  this  image,  a 
transform  that  has  been  known  for  a  half-century.  In 
fact  the  generalized  versions  of  the  Hough  transform7 
that  detect  other  shapes  such  as  rectangles,  ellipses,  or 
parabolas  can  he  shown  to  he  equivalent  to  a  generalized 
version  of  the  FRT.  In  this  case,  the  generalized  FRT 
maps  the  2-1)  shape  function  into  a  multidimensional 
feature  space.  The  Radon  transform  theory  Inis  not 
only  applications  in  pattern  recognition  but  also  in 
nuclear  medicine,  imaging  by  nuclear  magnetic  reso¬ 
nance,  determination  of  the  electron-momentum  dis¬ 
tribution  in  solids,  scattering  theory,  and  plasma  diag¬ 
nostics.8  10  While  digital  methods  for  both  const  ruc- 
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tion  and  reconstruction  of  the  Radon  transform  are 
available,11  14  optical  analog  reconstruction  techniques 
have  also  been  discussed.1  ‘  19  A  recent  survey  of  both 
coherent  and  incoherent  optical  reconstruction  tech 
niques  is  available.*0  More  recently,  Barrett-1  has 
discussed  use  of  incoherent  optical  FRT  processors  for 
performing  multidimensional  convolution,  correlation, 
and  spectral  analysis  using  a  series  of  I -I)  optical  Radon 
transform  processors.  In  tins  paper,  coherent  optical 
implementations  of  the  FRT  are  discussed.  The  pur 
pose  ol  the  optical  processors  is  to  he  used  as  a  band¬ 
width  compressor  of  optical  shape  functions  as  well  as 
in  t he  use  o!  11)  coherent  opt icul  signal  processing  of 
multidimensional  signals. 

II.  Radon  Transform 

The  FRT  can  he  defined  in  a  number  of  spaces  and 
in  an  arbitrary  number  of  dimensions. -v  lf/u,\  )  is  the 
2*1)  input  distribution,  the  FRT'  of  this  f  unction  is  de¬ 
fined  as,  using  vector  notation, 

=  j  /lx )Aip  -  x  ■  { idx,  lit 

where 

x  =  *  a ,  +  v  a  > .  i  Ji 

£  =  cusfhi,  +  sin^a,  -  ^a,  +  £  a,  i  ll 

The  transform  parameters  (/>,/ M  can  he  considered  as 
defined  on  a  rectangular  grid  in  the  (pjf)  plane.  A  point 
(fh)M o)  in  the  transform  plane  corresponds  to  a  line  in 
the  x -y  plane  where  />,»  is  the  shortest  distance  between 
a  point  on  the  line  and  the  origin,  and  ft,,  is  the  angle 
lormed  by  this  shortest  chord  to  the  line  measured  from 
the  positive  x  axis  (Fig.  I).  Conversely,  lines  passing 
through  the  point  (.v,i,y0)  in  the  .v-y  plane  correspond 
to  the  sinusoidal  curve 

/»  =  i  ii  i  \  „  Mfi^  ( ^ ) 

in  the  (f>J0  plane.  'These  properties  are  the  basis  of  the 
line  detection  capability  of  the  FRT.  While  the  present 
discussion  deals  with  line  detection  in  noiseless  images, 
an  analysis  of  curve  detection  in  noisy  images  is  avail- 
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able.*1 Jl  Some  of  the  elementary  properties  of  the 
FKT  are  now  described. 

Translation  of  the  input  transparency  by  a  distance 
Xo.  i.e.. 


/Iwfx)  *  fix  +•  +  Vu>.  (M 

yields  a  FRT 

fijp.O)  -  hr  +  £  ■  xo.tM  uii 

Taking  the  directional  derivative  of  the  input, 

t;\\)  ~  A  '  (7) 

where  a  is  a  constant  vector,  yields  as  its  FKT 

f;ip.O)  =  (a  •  (M/Vap.  (8) 

while  the  linear  multiplication  of  the  input 

f\ix)  ■  <a  *  x)/(x)  no 

vields  a  directional  derivative  in  the  angle  FKT 
plane: 

*  “U  •  Vp/  (10) 


y 


Fig  1.  Cartesian  input  ami  radon  transform  plane  coordinate 
geometries.  The  shortest  distance  to  the  straight  hne  is  p,  while  6 
is  the  angle  formed  by  the  shortest  distance  cord  and  the  i  axis. 
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Fig.  2.  Coherent  Fourier  optical  lens  configuration  used  to  record 
the  FKT  of  a  2  I)  input  The  input  is  located  in  plane  /*,.  A  vertical 
slit  is  placed  in  the  Fourier  transform  plane  i*i,  while  the  output  FRT 
is  recorded  in  the  inverse  Fourier  transform  plane  P\.  For  a  fixed 
angle  of  the  input  in  the  plane  P i.  a  Cartesian  coordinate  slice,  with 
FKT  coordinates  as  a  Cartesian  grid,  is  generated  at  the  output  plane 
P i-  Rotating  the  input  in  plane  P\  by  an  incremental  angle  while 
simultaneously  linearly  translating  the  output  storage  material  in  the 
output  plane  P.\,  a  new  slice  of  the  FKT  picture  is  generated. 


Finally  the  FKT  of  t he  2-1)  convolution  of  two  inputs 
/i< x )  and  f2(x)  in  the  space  domain 

fix)  «  J  /i(y)/2(x -*  y)dy  (111 

is 

A p,fn  =  j'*  fii(i.nfAO.p  -  ndt,  «12» 

a  1-1)  convolution  of  the  FKT  of  the  inputs  in  the  Ibid  on 
transform  (KT)  domain,  This  property  of  the  KT, 
namely,  that  it  maps  2*1)  into  l -l)  convolution,  can  be 
extended  into  higher  dimensions.  Also  this  property 
allows  reduction  of  multidimensional  correlation,  con* 
volution,  and  spectral  analysis  to  be  performed  using 
l-l)  acoustooptical  signal  processing  architectures. 

III.  Optical  Methods  of  Generating  the  Hough 
Transform 

T  he  starting  point  for  the  coherent  optical  generation 
of  the  FRT  is  the  direc*  implementation  of  the  inte¬ 
gral 

hpM)  ~  XL  f{x,\  -  x  cost!  -  y  sintfldxtfv,  (II 

where  /U.y)  is  the  input  and/(p,0)  is  the  output  with 
a  delta  function  as  the  space-variant  impulse  response 
kernel  function.  A  particular  way  to  implement  tins  ns 
well  as  other  similar  kernel  functions  is  to  combine 
suitable  linear  and/or  circular  motions  that  will  map  the 
space-variant  kernel  into  a  space-invariant  kernel.  If 
the  coordinate  axes  U.v)  are  rotated  by  an  angle  0 ,  the 
new  coordinate  axes  (ii,r)  are 

u  =  i  eosfl  +  \  sin 0  e  =  —  x  sintf  -f  y  casfl.  (13) 

In  the  rotated  coordinates  the  input  transparency  is 
f~iu,i').  Thus  by  rotating  the  input,  the  output  FKT 
can  he  expressed  in  terms  of  a  space-invariant  kernel 
operation 

/(p.tf)  =  f~(u.v)b(p-u)dudc.  (14) 

This  integration  can  be  performed  either  in  the  space 
or  in  the  spatial  frequency  domain.  In  the  space  do¬ 
main,  for  example,  one  could  electronically  integrate  the 
rotated  function  along  the  i>  axis  and  then  use  a  2-D 
coherent  optical  processor  to  implement  the  second 
integration.  This  integration  can  readily  be  performed 
using  an  astigmatic  processor.25  Or  one  can  use  a  2-D 
multiplexed  holographic  processor  to  implement  di¬ 
rectly  the  space-variant  impulse  response  function  of 
Eq.  (1).  Alternatively,  since  the  kernel  of  Eq.  (14)  is 
space-invariant,  the  filtering  operation  can  he  per¬ 
formed  in  the  spatial  frequency  plane.  Letting  the 
Fourier  transform  of  the  rotated  input  transparency  be 
F~(zfw),  the  FRT  is 

/<p.0)  =  jy  F'(z.u')H~U,u;p,fh(l:dW'  (IS) 

where  the  filter  function  H~~  is 

-  ( l f'lx )6(i<> ) exp( (H?) 

This  filter  can  be  implemented  using  a  slit  and  linear 
phase-shifting  wedge.  An  experimental  implementa- 
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Fiji  3  Experimental  results  generated  using  the  optical  setup  de- 
sir»M  in  Fig  2.  In  in)  the  input  while  in  th>  the  FRT  slices  of  (a) 
are  shown.  The  vertical  axis  in  lb)  is/),  while  the  horizontal  is  the  0 
axis.  Here  the  angle  is  measured  from  the  ('artesian  >  axis. 


b 

Fig  5.  Experimental  results  generated  using  the  optica)  arrangement 
shown  in  Fig.  2.  In  (a)  the  straight  lines  input,  while  in  (b>  the  FRT 
output  slices  of  the  input  of  (a)  are  shown. 
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Fig.  4.  Computer  generated  FRT  slices  for  the  input  of  Fig.  3(a). 

tion  of  this  scheme  using  a  coherent  Fourier  optical 
processor  is  now  presented. 

In  Fig.  2  a  coherent  Fourier  optical  lens  configuration 
is  used  to  record  the  FRT  slices.  The  input  transpar¬ 
ency  is  located  in  the  plane  /V  A  narrow  vertical  slit 
prescribed  in  Eq.  (16)  is  placed  in  the  transform  plane 


P2.  The  1-D  linear  phase  shift  is  implemented  using 
a  second  Fourier  transform  lens  together  with  a  1-D  slit 
in  the  output  plane  P3.  In  the  film  plane  Pa,  for  a  given 
angle  0 ,  a  slice  of  the  FRT  plane  is  recorded.  By  ro¬ 
tating  the  input  in  plane  Plf  a  new  angle  is  selected. 
The  new  angle  is  recorded  by  linearly  translating  the 
film  to  a  new  Cartesian  grid  position  in  the  FRT  plane. 
In  this  fashion,  the  whole  FRT  output  plane  is  filled 
sequentially.  In  Fig.  3(a),  a  particular  input  containing 
both  horizontal  and  angled  lines  is  shown.  The  re¬ 
sulting  FRT  slices  are  depicted  in  Fig.  3(b).  The  ver¬ 
tical  is  p,  while  the  horizontal  is  the  axis.  The  FRT 
slices  are  spaced  45°  apart  with  zero  at  the  left  and  360° 
at  the  right.  In  Fig.  4,  results  of  the  computer  generated 
slices  for  the  input  of  Fig.  3  are  shown.  Note;  to  be 
consistent  with  the  experimental  arrangement,  the 
angle  0  is  measured  from  the  vertical  axis.  While  the 
results  are  sensitive  to  the  choice  of  the  origin  (see  Eq. 
(6)],  these  results  do  match  the  experimentally  obtained 
results  of  Fig.  3(b).  In  Fig.  5,  a  second  input  with  the 
corresponding  output  FRT  slices  are  shown.  Here,  as 
in  the  previous  figure,  the  straight  lines  are  thickened 
for  the  FRT  analog  computer.  While  the  ideal  HT 
descriptors  deal  with  idealized  line  segments,  practical 
boundaries  are  always  filled  for  better  detection.  Thick 
boundaries  are  also  necessary  for  the  optical  analog 
computer  to  improve  the  detection  SNR.  In  Fig.  6, 
computer  generated  FRT  slices  of  the  input  of  Fig.  5(a) 
are  shown.  With  the  identical  caveat  mentioned  pre¬ 
viously,  these  results  match  the  ex  peri  men  tal  results  of 
Fig.  5(a). 

Various  other  expressions  for  generating  represen¬ 
tations  of  the  FRT  can  he  derived  by  using  different 
representations  of  the  Dirac  delta  function.  These 
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representations  can  lead  to  new  optical  implementation 
methods  for  the  analog  evaluation  of  the  FRT  operation. 
The  kernel  of  the  FKT,  a  space  variant  impulse  response 
function,  is  a  function  of  four  variables.  One  may  iso¬ 
late  any  particular  variable  and  derive  new  represen¬ 
tations  for  the  FRT.  Consider  a  change  from  Cartesian 
U,y)  to  cylindrical  coordinates  where 

Jr  *  r  cos0  v  =  r  sin,;.  U7) 

In  cylindrical  coordinates  the  FKT  kernel  function  is 

=  6 Ip  -  r  cos (ff  -  (18) 

Using  elementary  properties  of  the  delta  function,  Kq. 
08)  can  be  expressed  in  terms  of  the  radial  variable 
as 

hir,*,p.O)  *  6|r  -  pi costW  -  <,M|/|cos(//  -  «/»l|  (  19) 

or  in  terms  of  the  angular  variable  as 

h{r,*:p,tf)  =  -  il  +  cos~*(p/r)|/r  sqr|l  -  (p/r)-|-  (JO) 

These  representations  of  the  delta  functions  are  now 
used  to  derive  new  expressions  for  the  FRT. 

Let  the  input  function  be  expressed  as  an  an¬ 
gular  Fourier  series 

fir,*  1=  V  /„(#*>  expt-puM.  (21) 

where 

fnir)  *  */jJr  fir,*)  ex\t(jn*)d*.  (22) 

Substituting  Eqs.  (20)  and  (21)  into  Eq.  (1)  we  have 


where 


lipJh  ~  expi~jn(hl„ip). 


mi 


/„  =  /„lr)|l  -  (p/r)-)_,/i  expl-;n  C!ts~‘(p/r)|(/r  (241 


The  FKT  is  expressed  as  periodic  Fourier  series  in  the 
angular  variable  and  a  forward  circular  transform  in  the 
radial  direction.  The  inverse  circular  transform  has 
been  discussed  in  the  context  of  computed  tomography 
and  in  terms  of  coherent  optical  implementation  by 
Hansen  and  Goodman.16  The  forward  and  the  inverse 
circular  transform  kernels  are  quite  similar,  and, 
therefore,  both  of  these  transforms  can  be  performed 
with  the  same  optical  configuration. 

We  now  use  the  second  radial  expression  of  the  kernel 
function  of  Kq.  (19)  to  derive  a  different  representation 
of  the  FRT.  Let  the  input  transparency  be  represented 
as 


where 


j[r,*\  -  \  ,nj 


J(» («,<,*>)  explur  )du  . 

hr 


(2M 


*  J  fir,*)  exp(— u  r)dr  (2fi) 

is  the  Laplace  transform  of  the  input  transparency  in 
the  radial  direction  and  Br  stands  for  the  Bromwich 
path  necessary  for  evaluation  of  the  inverse  Laplace 
transform.  Substituting  Kqs.  (25)  and  (19)  into  Eq.  (1) 
we  have  a  new  representation: 


where 


hp.0\  *  Hz*) 


j: 


(i  ~ (u  ,p,0)du , 


(27) 


(i~uv\pjn\  -/)  j 

1 11 


(1  —  u*)r tos_i(l/u  I  + 


X  exp (-puu  )du.  (28) 

Here  again  the  2-D  integration  has  been  decomposed 
into  two  1-D  integrations.  The  optical  evaluation  of  the 
Laplace  transform  operation  has  been  discussed,25  The 
second  integration  indicated  by  Eq.  (28)  can  be  per¬ 
formed  by  an  astigmatic  processor.26 

Another  method  of  optically  performing  the  FRT  is 
based  on  a  coordinate  distortion  technique.  Coordinate 
distortion  in  optics  has  been  discussed,27,28  and  in  par¬ 
ticular  use  of  coherent  optical  1RT  has  been  described 
by  Hofer.17  The  method  is  based  on  a  2-D  Cartesian 
form  of  the  kernel  function  ( 1 ).  Consider  the  following 
expansion  of  the  delta  function  kernel: 


b(p  -  x  rosfl  -  v  sin 0)  «=  *4? w  ^  exp(-;u  (p  -  x  costf  -  v 

X  sin 0)|d«\  (29) 

Using  this  representation  of  the  delta  function  we  can 
represent  the  FRT  as 


/W) 


Fiw  cos  Art'  smtf)  exp{  -jup)du\ 


where 


(30) 
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/ 


=  Jj  /I  a  )  e\p|/n  U  costf  +  \  sintf)|i/ju/\  (Hi 


Here  the  2-1)  Fourier  transformed  input  is  used  as  a 
building  block  to  perform  the  FKT.  The  2-1)  Fourier 
transform  is  performed  by  a  conventional  coherent 
Fourier  optic  processor.  By  performing  a  I  T)  Fourier 
transform  in  the  radial  and  imaging  in  the  II  direction  on 
the  2-D  Fourier  transformed  input,  the  FRT  can  he 
formed.  The  coordinate  distortion  using  an  astigmatic 
optical  processor  must  he  used  to  perform  the  radial  to 
Cartesian  coordinate  conversion  followed  by  a  standard 
astigmatic  optical  processor  that  images  in  one  and 
takes  the  Fourier  transform  in  the  other  direction. 

Finally  another  delta  function  representation  can  be 
useful  in  the  case  where  both  the  input  and  its  FKT  are 
limited  to  a  finite  region  in  both  the  input  and  FKT 
plane.  In  this  case,  both  the  input  and  its  FKT  can  he 
considered  a  periodic  function.  Consider  the  following 
resolution  of  the  delta-'*  function: 


oti  -  i'|s  4„  cosri ;rt /V.  cos/i  m 1.72) 

where 

a,,  -  1/7,  n  =  0  am  I  2  /.  n  s  0.  (7:0 

Substituting  Ki(.  (.72)  into  Ktj.  (1)  we  have  a  periodic 
function  description  in  the  p  direction  of  the  FRT: 

/tp,fb  *  H  An  l‘l»sn  7r/>/Cg„(0),  cm 

n  *u 

where 

gntfM  =  f  ( v  ,v )  cosn  n /l.[x  cos 0  +  y  sin(V|c/jrdv  (7.7) 

is  the  2-1)  cosine  transform  of  the  input.  The  2-1)  co¬ 
sine  transform  can  he  generated  using  the  conventional 
2T)  coherent  optical  processor.10 


IV.  Summary  and  Conclusions 

Analog  coherent  optical  generation  of  the  HT  shape 
descriptors  has  been  discussed.  HT  shape  descriptors 
are  one  of  the  most  efficient  boundary  descriptors  in 
pattern  recognition.  A  number  of  possible  analog  im¬ 
plementations  have  been  described.  The.optical  im¬ 
plementations  are  based  on  the  various  representations 
of  the  Dirac  delta  (identity)  function,  which  is  the 
space- variant  impulse  response  of  the  HT  kernel.  Since 
the  HTand  the  FRT  are  identical  for  binary  images,  the 
above  implementations  are  identical  to  the  optical 
generation  of  the  FRT.  Because  the  FRT  of  multidi¬ 
mensional  convolution  in  the  space  domain  is  a  1-1) 
convolution  in  the  Radon  transform  domain,  multidi¬ 
mensional  coherent  FKT  processors  can  act  as  optical 
preprocessors  for  coherent  1-1)  acoustooptic  signal 
processing  architectures.  Experimental  results  using 
2-1)  coherent  Fourier  optical  transform  lens  configu¬ 
ration  together  with  rotational  and  linear  motion  have 
been  presented.  A  good  match  between  computer 
generated  and  experimental  results  is  shown.  A  num¬ 
ber  of  other  coherent  optical  FRT  architectures  are 
proposed.  Some  of  these  architectures  are  similar  to 


those  proposed  for  evaluation  of  the  inverse  Radon 
transform  known  in  computed  tomography. 
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Estimation  of  Two  Closely  Spaced  Frequencies  Buried 
in  White  Noise  Using  Linear  Programming 

Jaroslav  Keybl  and  George  Fichmnnn 
Department  of  Electrical  Engineering 
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New  York,  N.  Y.  10031 

ABSTRACT 


Linear  programming  is  used  to  estimate  the  spectrum  of 
two  sinusoids  signals  close! y-spaced  in  frequency  buried  in  deep 
white  gaussian  noise  by  employing  a-priori  knowledge  of  the 
spectrum.  The  method  will  be  illustrated  by  a  number  of 
examples. 

INTRODUCTION 


In  the  calculation  of  the  spectrum  of  a  discrete-time 
signal  consisting  of  two  sinusoids  with  c 1 ose 1 y-spaced 
frequencies  embedded  in  white  gaussian  noise,  problems  arise 
when  only  a  small  portion  of  the  signal  is  available.  This  is 
known  as  windowing  of  the  data  and  becomes  evident  as  "leakage" 
in  the  spectral  domain,  i.e.  energy  in  the  main  lobe  of  a 
spectral  response  "leaks"  into  the  sidelobes  obscuring  other 
present  spectral  responses.  There  are  many  methods  used  to 
estimate  the  spectrum.  The  periodogram  fl]  performance  is  poor 
for  short  data  lengths.  The  SI ackman-Tukey  method  is  also 

hampered  by  spectral  distortion.  The  Burg  method  f3],  a  high 
f requency-resol u t i on  technique,  even  in  the  absence  of  noise, 
yields  spectral  line  splitting.  The  recent  algorithm  of  Cadzow 
[4]  outperforms  the  Burg  method  in  a  low  noise  environment. 

In  this  paper,  we  estimate  the  spectrum  of  two 
sinusoidal  signals  closely-spaced  in  frequency  by  employing 
a-priori  knowledge  of  the  Fourier  spectrum  of  the  signal  in  the 
form  of  linear  inequalities.  The  advantage  of  imposing 
constraints  in  spectral  restoration  process  has  been  pointed  out 
r 5 ] -  In  the  linear  programming  formulation,  there  are  many 
solutions.  Here,  we  select  a  solution  which  minimizes  the  1. 
norm  of  the  Discrete  Fourier  transform (DPT)  of  the  measured1 
signal  and  its  estimate. 

DETERV I  NATION  OF  THE  SPECTRUM 

We  assume  that  the  signal  estimate  s(k)  can  be 
represented  by  a  weighted  sum  of  past  signals 

p 

s  (k)  =  -  za(i)s(k  -  i)  +  e(k)  k  -  1,2,..  ,m  (1) 

i=l 

where  a(i)  are  unknown  weighting  coefficients  and  e(k)  is  an 
error  term.  This  expression  is  a  linear  pred ictor [f 1 .  By  taking 
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the  DFT  of  Eq . ( 1 )  ,  we  have 

5(n)  =  -  ?  a  { i )  exp{-j?7r  in/N)  S  (n)  +  E(n)  i<n<_  N  (2) 

i=l 

where  S  (n)  and  E  (n)  are  the  DFT  of  s  (k)  and  o  n.)  ,  i  1  v<>  1  y  . 

Because  S(n)  and  E(n)  are  complex,  air1,  nsii.’i  i.ulrr  's  fnimu’a,  we 
rewrite  Eq. (2)  as 

P 

SR(n)  =  -  r^a(i) [SR(n)cos  v  +  Sj(n)sin  v]  +  ER(n)  (3a) 

Sr(n)  =  -  a  (i )  [ST  (n)  cos  v  -,SR(n)sin  v]  +  Ej(n)  (3b) 

where  v  =  2^in/N.  We  have  2N  equations  with  p  unknowns.  When  the 
signal  is  real,  we  can  reduce  it  to  N  equations  with  p  unknowns. 
Assuming  that  two  sinusoidal  signals  are  present,  it  can  be 
shown,  the  a(i)  coefficients  are  the  pole  coefficients  of  the 
z-transform  of  two  sinusoidal  signals.  Therefore,  p  =  4  and  t-  hr 
range  of  the  a(i)'s  are 

-4<a(l)<  4,  -2<  a  (2)<  6,  a(3)  =  a(l),  a(4)  =  )  (4) 

Eq.(4)  adds  6  more  equations,  to  the  N  equation  generated  from 
Eq.(3),  for  a  total  of  N  +  6  equations  used  in  the  linear 
programming  formulation.  From  Eq.(3)and  (4),  we  can  now  solve 
for  the  a ( i )  coefficients  by  minimizing  the  1  norm  of  the  error 
-  E(n)  .  1 

NUMERICAL  RESULTS 

To  test  this  method  we  have  used  the  time  series 

s(k)  =  A1cos(2uf  1k)  +  A2cos(2Ttf?k  +  )  +  w(k)  (5) 

with  l<ki  J  and  w(n)  as  white  Gaussian  noise  with  zero  mean  and 
variance  a  .  The  two  sinusoidal  frequencies  are  normalized  so 
that  f  =  0.5  corresponds  to  the  Nyquist  rate.  The  individual 
sinusoidal  signal-to-noise  ratio's  (SNR)  is  given  by  201og  (A./ 
✓r0  )  for  k  =  1,2.  For  all _of  our  examples,  we  chose  the  signal 

amplitudes  A^  =  /7~,  A  =  *£o  and  the  signal  frequencies  f  = 

0.2168  and  f_  =  0.224R.  In  two  examples  we  introduced  a  1 
forty-five  degree  phase  difference  between  the  two  sinusoids. 

In  Fig.  1.  the  variance  of  of  the  noise  02  =  0.94  and 
M=  192.  This  corresponds  to  0.54  db  SNR  on  the  weaker  signal 
and  an  time-bandwidth  product  (TBP)  of  1.50.  The  spectrum 
calculated  using  the  periodogram  shows  random  fluctuations  and 
it  resolves  both  frequencies.  Our  method  yields  peaks  at 
frequencies  f.  =  0.2168  and  f2  =  0.2265  which  shows  the  ability 
to  resolve  the  frequencies  with  a  small  error  in  such  a  low  SNR 
environment  without  fluctuation.  It  should  be  pointed  out  that 
the  reason  both  peaks  are  of  equal  amplitude  is  because  we  are 
calculating  the  poles  which  make  the  function  blow  up.  In  Fig. 

2,  we  introduce  a  forty-five  degree  phase  shift  between  the 
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sinusoids.  Here  the  SNR  is  -  0.26  db.  U’e  see  that  there  is  no 
line  splitting  as  is  often  the  case  with  other  algorithms.  The 
peaks  occur  at  f^  *  0.2207  and  f?  =  0.2265,  again  resolving  the 
frequencies  in  aiow  SNR  environment.  Fig.  3  shows  the  estimated 
results  when  the  SNR  is  -  14.2  db.  The  frequencies  f  =  0.218P 
and  =  0.2305  yield  good  estimates  for  a  very  low  SNR.  In 
Fig. 4,  we  added  a  forty-five  degree  phase  shift  with  a  SNR  -  15 
db.  V?e  see  that  the  periodogram  is  unable  to  resolve  the  two 
frequencies  but  our  method  peaked  at  f^  =  0.2188  and  f^  = 

0.2324,  Again  we  note  that  there  is  no1line  splitting.  In  the 
next  set  of  exper iments , the  TBP  is  reduced.  In  Fig. 5,  the  SNR  is 
0  db  and  the  TBP  is  1.00.  The  periodogram  shows  peaks  while  our 
method  gives  the  frequency  estimates  f =  0.2160  and  f  =  0.23C5 
showing  that  we  get  good  quality  estimates  even  when  the  number 
of  samples  is  reduced.  Finally,  in  Fig. 6,  we  plot  the  spectrum 
for  a  SNR  of  -  14.5  db.  The  estimated  frequencies  are  f  = 

0.2168  and  f^  =  0.2363.  The  larger  frequency  error  can  Be 
attributed  to  the  very  low  SNR  and  the  small  TBP  environment. 

SUMMARY 

We  have  shown  a  new  method  for  determining  the 
spectrum  of  two  sinusoidal  signals  closel y-spaced  in  frequency, 
where  the  weaker  signal  has  a  very  low  SNR.  We  show  that  the 
method  is  not  affected  by  spectral  line  splitting  when  the  two 
sinusoids  have  a  phase  difference  of  45  degres.  Computer 
generated  results  of  this  method  have  been  presented. 
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A  nn*t  In  »•!  nf  restoring  I  lit*  dm  rcte  K*  »n  tut  transform  *  I  >1* T>  spe«  (rum  of  a  tiilft.it  limi  limited  (Pl.t  image  I  mm 
,i  narrow  i»tisrr\ .it ion '•rjatu  nl  of  tin  PI.  image  is  ptesriiicd.  The  PL  sp(-*  i r.d  n  -toral i«>n  pfu»  i*ss  i- t  hr  dual  *>f 
tin*  more  common  PL  1111.4:1’  rest or.it n*n  prongs  w  it h  the  rules  of  the  treqm-mv  and  sf»a<  e  re\ er-i-d  Appltt  at ions 
ol  .1  spn  t  rum  riwior.it  ion  un  lade  im  reusing  t  hr  I  nil  t  >1  \  lew  ul  exist  mg  imaging  and  ext  ra»  t  in*:  prec im-  Ire 

tpirm  x  t’oinponrnt  x  1  >1  a  large  PL  image  lw  using  on  1\  a  '•mall  segment  ol  thr  rntirr  intake.  This  111H  hot  1 1  011M  aLo 
he  employed  for  im.ip'  data  i  t»mprr->it>n.  w  hit  h  is  ol  mtrrrsl  111  digital  \  idro  appl it  at  urns.  Several  diltrrrm  rs  he 
tween  the  implementations  of  the  image  and  tin  sped  rum  rest  oral  ion  proiesse-  are  drserihed  The  estimate  is 
constrained  to  tiave  an  upper  hound  on  the  numht  r  of  frequency  components  »  ont.uneil  in  (lie  Fourier  spn  tnmi 
The  hound  is  the  number  ol  samples  at  tpnred  at  the  Nvquist  rate  tor  the  length  ol  the  image  The  magnitude  <*t 
the  PKT  spectrum  is  .iLm  hittindetj  l’hes,«  lon-liaml-  th  line  a  large  number  o|  possible  sohillons  Thr  di  mo  , ! 
solution  is  then  sflfi  ted  sin  h  that  the  distant' e.  detmed  in  .1  lutu  lion  theoretic  sense,  between  the  mra-ured  .md 
the  estimated  image*  i-  an  optimum  A  number >1  -m  li  measures  are  investigated  Numerical  e\t»eriments  di.*w 
that  this  appro.it  h  v  telds  re>ults  that  art*  highly  immune  to  measurement  noise. 


INTRODUCTION 

The  finite  aperture  of  any  physical  imaging  system  eliminate  s 
the  high  spatial-frequency  components  of  tlu*  object  I  mm 
appearing  in  the  image.  The  lack  of  high  frequency  drLiil 
results  in  a  loss  of  resolution  m  the  observed  image.  It  lias 
been  shown  that,  for  an  object  of  finite  extent,  an  exact  res 
torat ion  of  the  object  from  the  diffraction  limited  (Did  image 
is  possible.  The  restoration  of  the  DI.  object  can  he  viewed 
as  a  coot  inuat  ion  of  the  spat  ial  Fourier  sped  rum  beyond  t  he 
spatial  cutoff  frequency  imposed  by  the  DI.  system.  No 
merical  met  hods  id  DI.  image  restorat  ion  are  highly  unstable 
in  t  he  presence  of  measurement  noise.  The  restoration  pro 
cess  can  be  stabilized  by  imposing  additional  constraint^. - 
Many  restoration  methods  are  available,  such  as  those  of 
Kriedenr’  Schell/  Hiraudd  -lansson,  *  and  Hurg.*’  Recently, 
Howard1  demonstrated  a  computal iouallv  inexpensive 
method  using  a  least  -squares  approach.  This  approach  w  as 
further  elaborated  on  by  Rushforth  cf  uls  A  new  method  of 
obtaining  an  increase  of  image  resolution  using  linear  pro 
gramming  teehniques  was  recently  given  by  Mammone  and 
Kichmann/ 

In  t  his  paper,  the  dual  of  the  I  )l .  image  restoration  problem 
is  considered.  Here  t  he  exl  rapolat  ion  of  a  I  mil  e  segment  of 
t  he  I  >1 .  ti  e.,  spatial  haudiimited  1  image  data  in  the  present  r 
of  measurement  noise  is  presented.  'The  restoration  of  the 
discrete  Fourier  transform  (DKT)  spectrum  may  he  inti  r 
prrted  as  an  extrapolation  of  the  truncated  spatial  image. 
Applications  of  sped  rum  restoration  include  increasing  the 
field  of  view  of  existing  systems  and  extracting  rxm  1 

frequency  components  of  a  large  DI.  image  when  onlv  a 
smaller  image  segment  of  the  entire  image  is  available.  I  hi- 
reduction  woitlrl  he  also  of  interest  m  image  data  compres  ami 
applications,  such  as  the  transmission  ,»t  digital  video  im.i:  <  ■ 
Although  spectral  restoration  yields  the  same  mathemat h  d 


formulation  as  that  obtained  for  DI.  image  restoration,  with 
the  role  of  space  and  spatial  frequency  reversed,  there  are 
significant  differences  as  well.  Here  the  unknowns  are  Fou¬ 
rier  coefficients,  which  are  complex  numbers  This  fact  ap 
pears  to  double  the  numerical  complexity  ot  the  problem. 
Also,  the  circulanl  convolutional  operator  that  results  trurn 
approximating  tin*  exact  Fourier  transform  with  a  DF  V 
spectrum  leads  to  additional  consideration-  Finally,  al¬ 
though  the  restored  DI.  object  is  tmnncgai i\e.  the  rcsloted 
DKT  spectrum  need  not  he  real  or  nonnegativc.  Thus  the 
powerful  nonnegai ivity  constraint  used  in  t  he  l  >1 .  image  res- 
torat  ion  cannot  lie  directly  used. 

The  motivation  of  this  research  is  to  provide  high  resolution 
frequency  est  imates  of  data  available  only  on  an  meomph  * . 
observation  length  by  vising  the  computationally  ctln  1.  m 
last  Kourier-t ransform  afgont Inns.  To  simplify  t  fn-  disv  u- 
sion,  one -dimensional  images  art*  considered.  Inr  .1  DI. 
I  band  limited)  image  truncated  to  lie  within  an  interval  N.  the 
spatial  frequency  resolution  is  limited  b\  the  univrtaintv 
principle/*’  The  lack  ol  spectra)  resohition  is  related  to  the 
infinite  limits  of  integration  tn  the  delmition  ot  the  Fourier 
I  r.  1 11st  or  m.  Since  tin*  image  is  available  over  a  tin  Me  ob-er 
vatioii  It'iigth,  the  calculated  Fourier  transform  is  disinrted. 
It  is  convolved  with  a  mih  turn  turn  with  mam  Inhe  width  /  L 
For  the  case  ot  two  dMIerent  sinusoidal  images,  the  vnui  r 
taint y  principle  states  that  the  dtlfcreme  111  the  frequency 
between  the  images  cannot  be  less  than  />  ;  otherwise  tin  re 
will  be  a  significant  overlap  of  the  transforms.  This  overlap 
would  not  permit  the  two  separate  responses  to  be  distm 
gtiishaidc  Thus  the  frequency  resolution  is  (muted  to 

Ih  >  |  ,'N.  U> 

In  main  practical  situations,  the  image  is  not  available  lor  an 
interval  * >|  sufln  n-nt  length  N  lor  the  desired  frequency  res 
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nl  lit  inn  /)< .  It  is  then  desirable  1 1 » pn»t  ess  tin*  t  rttrn  of  the 
observed  signal  stu  b  that  the  spectral  rrsnlutinii  ts  **r 
than  that  given  hy  f'ornuila  1 1 ).  It  du  re  is  no  source  <>f  hum 
surement  error,  it  is  well  known1  that  the  spectrum  can  he 
exactly  restored.  However,  there  is  alwavs  m »i^t-  and  or  error 
present  because  of  the  finite  numerical  precision  neiessjrv 
to  carry  out  t  In*  computation,  if  for  no  »<t  her  reason.  There 
fore  noise  or  measurement  error  must  he  considered  as  part 
ot  the  lormulalion  of  tin*  problem.  In  this  case,  the  nMora 
tion  process  becomes  highly  unstable. 

It  h  as  been  found'1  that  the  imposition  ol  a  priori  con 
Straints,  such  as  nonnegativity  ot  the  estimate,  will  siahili/e 
the  restoration  process.  Since  the  maximum  number  of  in 
dependent  equations  that  can  he  generated  is  equal  to  the 
number  of  spatial  samples  r  required  at  the \vquis!  rate,  we 
shall  restrict  the  spectrum-restoration  method  to  situations 
in  which  the  number  of  frequency  components  is  less  than  or 
equal  to  r.  This  limit  on  the  dimension  of  the  estimate  lias 
been  obtained  in  several  other  studies.111  Further,  the  si¬ 
nusoidal  frequency  components  will  he  assumed  to  In*  har 
monies  of  the  same  fundamental  frequency.  'Die  latter  as¬ 
sumption  is  not  necessary  hut  is  made  to  facilitate  the  use  of 
(lie  fast  Fourier  'transform  (FFTf  approximation  of  tin* 
Fourier  spectrum.  The  periodic  frequency  spectrum  is  ob¬ 
tained  by  zero-padding  the  truncated  space  sequence  such 
that  the  total  numher  of  elements  is  /..  an  integer  that  is  a 
power  of  2.  'This  step  allows  a  liner  frequency  scale  to  he 
generated. 

The  optimal  data  titling  employs  linear  programming  il.lM 
techniques.  The  LI*  method  provides  both  cost  and  time- 
e/fedive  wavs  of  selecting  the  optimal  r  tuple  estimate.  LP 
methods  have  been  used  previously  tor  the  general  image- 
restoration  problem. v 1 1  M  The  advantage  ot  using  the  Id' 
techniques  lor  the  unstable  PL  image  restoration  problem 
lias  also  been  addressed.1'  In  this  paper  we  demonstrate  the 
advantage  of  this  approach  for  the  rank  deficient  spectrum 
restoration  problem.  Since  the  rank,  the  number  of  inde¬ 
pendent  equations,  is  less  then  the  number  ot  unknowns, 
many  solul ions  exist.  The  deficiency  of  the  rank  occurs  be 
cause  of  the  I )KT  low  pass  tilter  effect  that  eliminates  all  but 
r  discrete  frequency  components.  This  effect  is  in  contrast 
to  the  linear  discrete-convolution  case.'*  which  yields  a  system 
of  equations  that  has  a  lull  rank.  However,  because  the  ad 
jacent  equations  are  almost  identical,  when  a  finite  arithmetic 
machine  rrpresentat  ion  is  used,  the  system  of  equations  is  ill 
posed,  i.e.,  nearly  rank  deficient. 

FORMULATION 

For  the  following  discussion,  it  will  he  helpful  to  define  neg¬ 
ative  spatial  and  frequency  samples.  The  lirst  </  *2  +  U  el¬ 
ements  of  llu*  sequence  [n  or  /*’/,  correspond  to  the  positive 
spatial  or  frequency  samples,  respectively,  in  ascending  order, 
am!  tlie  remaining  1 /./2  —  I)  elements  correspond  to  the 
negative  spatial  or  frequency  samples,  respectively,  in  de¬ 
scending  order.  This  is  consistent  with  the  idea  of  extending 
the  sequences  in  a  periodic  manner.  The  spatial  truncation 
operation,  the  model  for  a  short 'observation  segment,  is 
characterized  hy  multiplication  of  the  spatial  sequence  hy  the 
unit  rectangle  sequence  //,,  where  r  is  llu*  number  ot  adjacent 
unity  elements  centered  about  the  zero  element  and  with  all 
other  elements  are  equal  to  zero.  The  OFT  <»1  the  rectangle 
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sequence  i-*  real  and  even.  Miiir  l\,  is  real  and  even,  and  is 
mii  nl. if  !u  ,t  .oupli  d  ‘-mi  lorn  i  imii  Let  1  hr  column  vert  or  h 
be  drille  d  a  \i  rt»*r  wb"M-  rlrinrnts  arr 

h  -PFTJ/L!  (21 

We  shall  use  thr  t  in  ul.mt  n»nvolutM*n  theorem1  to  examine 
tin*  rriatinU'hip  In  luo  ii  the  tnimatrd  and  tin*  exait  se 
queners.  Thr  t  in  ul.ini  » i*nv oh  it  h  >n  mat  ri\  tor  the  disrrete 
spatial  truncation  operator  Ii  is  given  hv  thr  cm  ulant  ma¬ 
trix 

h,t  /i,h.  ...  /i, 

Ill  i  h,h\  *•*  hf  : 

II  =  ;  ;  '*  .  CP 

hi  h  j  fi  (  ...  hu 

and  therefore 

f;*IIF  +  .V,  (4) 

where//  and  Fare  the  DFT’s  of  the  measured  and  the  actual 
sequence  vectors,  respe<  lively,  ami  .V  is  the  complex  error 
(noise)  vector.  It  is  well  known1  that  the  eigenvalues  of  a 
cin  ulant  mat ri\  II  are  tlu*  coefficients  of  t  fie  I  >FT  of  tin*  lirst 
row  ot  the  circulant  matrix  and  l lie  eigenvectors  are  the  DFT 
basis  vectors.  Therefore  the  rank,  which  is  equivalent  tot  fie 
number  of  nonzero  eigenvalues,  ol  II  is  equal  to  the  number 
of  nonzero  elements  e  in  the  rectangle  sequence  /i\  Since  the 
vectors  (i.  i\  and  .V  are  complex,  Kq.  tp  represents  2/.  real 
scalar  equations.  Hut  all  the  spatial  sequences  are  real,  and 
therefore  the  DFT  sequences  must  have  symmetry.  The 
symmetry  in  turn  rntrodoces  redmidannVs  into  Kq.  i  ll.  'The 
elimination  of  t his  redundancy  reduces  the  2b  to  b  real 
equations.  In  order  to  obtain  these  l.  equations  we  must  ex¬ 
amine  some  properties  of  llu*  DF  T. 

Similar  to  I  he  properties  of  the  Fourier  transform,  the  OFT 
ot  a  real  sequence  possesses  an  even  real  and  an  odd  imaginary 
part.  Since  we  wish  to  use  the  computationally  efficient  radix 
tvvo-FFT  algorithm,  the  sequences  must  consist  ot  an  even 
numher  of  points.  The  exist  mg  symmetry  can  he  displayed 
in  the  following  way.  There  is  one  element  lor  the  dc.  the 
zeroth  component  and  1/  2  -  1)  negative  as  well  as  positive 
elements,  and  one  remaining,  the  b  2.  element.  This  can  he 
seen  by  noticing  that  the  /..  2th  row  of  the  PFT  matrix  is  a 
sequence  of  alternating  positive  and  negative  unities  Finis 
the  center  and  the  tic  frequ«»ncy  components  are  always  real 
I  he  remaining  lb  -  2)  elements  consist  of  complex  conjugate 
pairs  between  the  lirst  and  the  last  elements  and  tlu*  second 
and  the  second -to-lasl  elements,  and  so  on.  There  are  (/.  2 
F  Ddfsiim  J  rra)  and  I /..  2  -  1)  imaginary  element''.  We  form 
the  concatenated  sequence  of  the  distinct  elements.  Let  the 
real  and  the  imaginary  parts  of  the  tth  element  he  denoted  hy 
the  subscripts  l\i  and  h.  respectively; 
then 

/»  -  (*!< i.  •  •  P’/o/ .  j 4 i>.  /»/,/  j-n.  *  *  *,  H/d* , 

F  a  | F/,’u,  I'm,  •  •  •,  Fin i  ■  i>.  F/u  j  -  n.  ■  •  F/ \\l , 

;V  -  lA'/iu.  Ar/n.  •  •  •.  A//,/ /■_»*  1 1,  ,VM  m  -  n,  *  *  *.  A i  A1 .  Cr>) 

where  tlu*  superscript  7*  stands  for  the  transpose  operation. 
The  relation  Kq.  (  It  can  now  he  written  in  rediued  lorm 

(1  =  IlF  +  N\  |i>) 


}l.  hLiinniuiiv  {tin)  l*.  )\>i  huumi 
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where,  if  we  decompose  the  degradation  matrix,  then 


H 


a  ;  it 

i 

t 

“T-l-n-- 
I  i  I 

I  I  « 

«•  !  <’  \f>\  I) 


(7) 


where  A  is  a  i/-/2  +  D  X  (L  'J  >  !  >,  It  is  a  t/./2  +  I)  X  l/./’J  - 
1 t,  (A  and  I)  are  [  l.  ‘1  —  1 )  X  ( /  ‘  1  ~  1  >  stibmat  nee*  <>)  H  and 
t\  and  p  are  unused  column  vectors  It  we  denote  the  mill 
vector  and  matrix  hy  <1  and  0,  respectively.  and  delihe 


\\  =  lulUdii, 

(SI 

0  =  Hiu(*!, 

<‘d» 

1)  =  jlMlD'T 

<  !<d 

where  the  superscript  / (  denotes  reverse  order,  then 


!? 


A  +  H  |  o 


» 

i 


o  |  n  -  r 


Thus  Kq.  <f>)  represents  /,  real  initiations  in  the  'J/.  variables 
b  and  S  Here  !1  denotes  the  space  truncation  operator  in 
the  DET  domain.  This  underdetermined  system  of  equations 
has  many  solutions.  I''<»r  example,  the  pseudo  inverse solution 
might  he  used, 1,1  although  this  would  produce  an  uncon¬ 
strained  estimate.  In  oriler  that  we  imp  he  aide  to  impose 
cont raints  on  and  simultaneously  choose  an  optimal  solution, 
that  is,  one  that  is  close  in  some  sense  to  the  measured  intake, 
we  use  a  method  of  optimization  known  as  linear  pro¬ 
gramming. 


CONSTRAINED  OPTIMAL  SOLUTIONS 

I’lie  haste  linear-programming  (LIM  problem  is  formulated 
in  the  follow ing  way1  : 

Minimize  the  cost  function  c  f  x 

Subject  to  \x  <  b, 

x  i  0.  l  KM 

where  c  and  v  are  tlie  known  cost  and  unknown  Af  dimen¬ 
sional  variable  vectors,  respectively,  h  is  an  r  dimensional 
constraint  wet  or.  and  A  is  .in  r  x  ,\f  constraint  matrix.  The 
most  often  used  Id*  technique  is  (he  sin-;  method.1 1  The 
output  obtained  Irmn  a  simplex  l  J*  algoru  r,m  will  correspond 
to  one  ol  l  he  following  l  hree  situat  loiis:  There  is  no  feasible 
solti  I  ion,  that  is.  all  I  he  1 1  nisi  raints  »  an  in  >t  simultaneously  he 
Stit  isl led.  »*r  t  he  opi  unal  si ■luhon  is  I rasi hie  hul  is  unbounded; 
the  third  pos»ihihl\  »>  that  the  solution  is  feasible  and 
bounded,  and  then  an  optimal  solution  is  provided. 

The  lir-t  >t ep  in  the  simplex  method  is  to  euiiwrt  the  in 
equalitv  constraints  tnio  equalities.  This  step  is  accomplished 
bv  lilt  r«*du»  mg  ii<  (<  1 1 1  ioi  ia  I  variables.  vailed  slack  variables, 
Tin  result  hi;;  s\stem  «'|  » qua!  mils  w ill  have  more  unknowns 
I  ban  eipialinfis,  h-ading  lo  an  umlertletr  limned  system  ol 
equations  Kn»in  the  manv  pos»ible  solutions  i (u*  simplex 
method  seeks  a  solution  that  minimize*  (he  cost  (unction. 
The  I  »asit  premise  *1  Id’  is  that  the  solution  vector  \  will 
consist  ot  at  most  r  non/ero  elements,  where  r  equals  tlie 
number  of  independent  equations,  and  I  he  remaining  M  -  r 
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elements  w  ill  lie  equal  to  zero.  In  a  vector-spare  interpreta¬ 
tion  ol  LI*.  I  he  inequalities  define  intersecting  hyperplanes, 
which  lorin  a  region  ot  possible  solutions.  'This  region  is  an 
r  dimensional  polygon  or  simplex.  The  simplex  will  always 
lie  m  the  postive  octant  of  the  vector  space  The  cost  function 
also  defines  a  hyper  plane  lor  a  fixed  cost  value.  Kor  this  value 
to  lie  minimum,  the  cost  function  by  per  plane  must  lie  on  an 
edge  ut  the  simplex.  The  coordinates  ot  each  edge  contain  at 
most  r  numbers.  The  simplex  method  is  an  iterative  tech¬ 
nique,  whit  h  begins  with  a  feasible  solution  at  an  arbitrary 
edge  of  the  simplex  and  progresses  to  an  adiacent  edge  in  such 
a  way  as  to  yield  a  redm  lion  in  the  cost  value  until  im  further 
(eduction  is  possible.  The  coordinates  of  the  resulting  edge 
yields  the  optimal  solution. 

In  a  LI*  tormulaliou  «•!  tlie  spectral -estimation  problem,  we 
now  take  Kq.  <b>  lo  lie  part  of  the  constraint  matrix.  To 
guarantee  that  all  the  variables  are  nonnegative,  a  necessary 
requirement  in  LI*,  we  next  define  the  vectors  in  Kq.  ltd  in 
terms  ot  th»*  nonnegative  component  vectors  F* ,  F  ,  .V  * ,  and 
X  such  that 


r  =  (i:» 

and 

N  =  S'  -  S~.  (14) 

If  vve  examine  the  sum 

S'  A.V  -  =  !  <7  -  II  F\,  (15) 

we  observe  that  Kq.  (15i  is  the  l\  norm  of  the  error  between 
the  spectrum  ot  the  measured  signal  <7  and  the  spectrum  of 
the  spatially  truncated  image  estimate  l  l  his  i>  (he  error 
measure  to  be  minimized  Thus  the  LI*  formulation  ot  the 
/i  spectral  estimation  problem  is 

Minimize  X '  +  S  ‘ 

Subject  to  (!  =  1!  (F*  -  F~)  +  X*  —  .V~, 

1  —  *  in. n  » 

r  <  >  o.  no 

The  /,  norm  is  not  the  only  measure  of  error  that  can  be 
minimized  'The  /,llt  norm,  or  maximal  deviation,  can  also  be 
employed.  This  can  be  seen  by  delining  a  scalar  F  to  be  t he 
absolute  value  of  the  maximum  element  of  the  error  vector 
given  hy  Kq  f  I5i.  Thus  the  Id*  formulation  of  t lie  spec¬ 
tral-estimation  problem  is 

Minimize  E 

Subject  to  (i  <  II  l  /'♦  —  F  )  -f  /•>„, 

r,  z  H  (/  •  '*  -  F  )  -  lu ,, , 

.-b 

F  <  .  F'J-.E  >  0,  (IT) 

where  e„  is  an  L-tlimensional  unit  vector. 

t  he  /  *  norm,  or  t  he  sum  of  t he  squares  of  the  error,  can  also 
he  minimized  hy  using  a  formulation  similar  to  that  ol  t  hi'  / , 
norm.  I  li  re  quadratic  programming  methods  are  u-rd. 
ljuadral ic  programming  techniques  provide  an  est  imate  with 
at  most  r  positive  elements,  with  all  other  elements  as  zero. 
In  general,  quadratic  programming  techniques  use  conqni 
tationallv  ellicient  Id*  methods  but  on  larger  augmented 
matrices.  'These  method*  are  well  known,  and  manv  puli 
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lished  computer  routines  are  available.1*  There  is  usuallv 
m«»re  than  mu*  least -squares  estimate  to  an  undetermined 
system  of  equations  such  as  Kq.  M).  For  example,  the 
pseudo  inverse  solution1  ’  provides  an  une -oust rained  least 
squares  estimate  with  the  lowest  l->  norm  as  well  as  the  mini¬ 
mum  squares  error.  Howard'  lias  demonstrated  a  eon 
strained  least  squares  approach  lor  the  dual  problem  incor¬ 
porating  a  penalty  function  to  force  the  solution  to  beeome 
noiincgnt ive.  Kushfnrl It  rt  «/. M  developed  this  approach  lor 
an  ill-posed  DL  intake restoration  problem. 

The  estimates  obtained  by  optimizing  the  / j.  /■,  and 
norms  are  maximum  likelihood  estimates  when  tin*  noise  n 
modeled  as  a  random  variable  with  uniform,  (laussian.  or 
exponential  probability  density  functions  t pdf’s),  respec¬ 
tively.*"  It  was  found,  however,  that  tor  the  low  noise  situa¬ 
tion  of  interest,  the  pdf  of  the  noise  dit!  riot  appreciably  allcct 
the  estimates.  In  terms  of  increased  frequency  resolution,  the 
/i  norm  consistently  offered  the  best  estimate,  irrespective  of 
the  noise  pdf.  This  fact  can  be  attributed  to  the  relatively 
fewer  numerical  operations  needl'd  as  well  as  to  the  strongest 
bounds  on  t he  individual  sped ral  components  for  t he  l\  es 
timate.  The  constrained  minimal  /_*  norm  using  quadratic¬ 
programming  techniques  requires  many  more  numerical  op¬ 
erations.  Also,  optimizing  the  U  norm  is  equivalent  to  mini¬ 
mizing  the  total  noise  spectral  energy.  Whereas  t  his  is  a  de¬ 
sirable  feature  in  general,  for  the  purpose  of  increased  fre¬ 
quency  resolution  it  tends  to  be  counterproductive.  1‘he 
minimization  of  the  global  spectral-error  noise  energy  dis¬ 
tribution  tends  overly  to  smooth  the  estimate,  leading  to  a 
decrease  in  frequency  resolution. 

NUMERICAL  RESULTS 

The  observed  spectrum  is  a  smoothed  version  ot  the  desired 
spectrum.  Figure  1  illustrates  tin*  inadequacy  of  simply 
taking  the  FI  T  of  a  truncated  sequence.  I  lere.  as  well  as  in 
all  subsequent  figures,  the  magnitude  of  the  Fourier  spectrum 
is  plotted.  The  dotted  line  represents  a  .'12-point  FFT  of  2d 
samples  of  a  cost/c/ .r  !  sequence,  where  »■/.  is  the  fundamental 
frequency  2n7d2.  Here  only  seven  data  samples  have  been 
replaced  by  zeros.  The  solid  line  illustrates  the  exact  spec¬ 
trum  of  the  full  d2-samp)e  sequent  e.  The  lack  of  resolution 
conceals  the  fact  that  only  the  fundamental  frequency  is 
present.  The  size  of  the  FFT  as  well  as  t  he  number  of  samples 
acquired  is  known  a  priori.  This  information,  togt  t  her  wit  h 
the  bounds  on  the  amplitude  spectrum,  is  used  to  estimate  the 


Fig  I  The  aj  point  KKT  iii  sample  pomK  ot  t  /  r  Hdntled 
dashed  curve l  and  the  exact  DKT  spectrum  tsolnl  line) 


l{  M  unit . .  <  f\h  f.m.mn 


Fig  2  Tin*  r»"'t,,r.ili«in  "t  ih>  i  *1  I  »  i  u.-m  1 1  ♦  point 

FI  I'  «>l  7  vp.uui  dom.nn  '.«mpi«-  |><-niO  « d » *  -i  d.i-hrd  <  irw  i  and 
ihe  estimated  sju-i  train  is,,ii<(  Inn*' 


Fig  ■!  K\t r. ipol.it ii*n  <•!  » *s«;.  »  i  '  in  tin-  sp.m.d  d'  nu:n  nt  line 

curve)  from  7  sample  points  l solid  lint-i 

actual  sped  rum.  Figure  2  illud  ral  es  1  he  perlormam  t  o!  t  he 
/ j  norm  when  only  7  of  the  !VJ  samples  o|  tin*  hindament.il 
cosinusoid  are  available.  The  dotted  line  curse  represents 
t  lie  overly  smoot h  estimate  of  tin  direct  FFT  spectrum  I  he 
minimal  l\  estimate,  the  dotted  line,  is  close  to  the  <irigin.il 
spectrum.  The  doited  line  curve  lies  <  nmplelelv  on  the 
solid  line  curve  Thus  a  complete  restoration  is  possible  even 
when  only  7  of  T2  samples  are  available.  Figure  a  illiM  r  ‘s 
I  lie  ext  rapolated  space  domain  sequence  <  or  respi  Hiding  1  o  I*  ig. 
2.  'Flu*  acquired  samples  are  illustrated  hv  Hu1  solid  line 
segment  and  the  extrapolated  curve  hv  tin'  dashed  hue 
curve. 

The  method  stands  up  quite  robust  |\  to  measurement  noise. 
Figure  J  illustrates  the  restoration  of  tin*  fundamental  co¬ 
sinusoidal  waveform  I  mm  I  d  samples  <<i  a  .12  point  FI1’  P  w  it  li 
white  ( muss  inn  noise  added.  The  \ .main  e  nt  the  noise  is  one 
tenth  of  the  amplitude  nt  the  function.  'Thus  an  increase  of 
the  order  of  a  factor  of  2  in  resolution  is  obtained  even  in  t  lit* 
preseme  nt  signjluant  measurement  noise.  The  corre¬ 
sponding  spatial  extrapolation  is  depicted  m  Fig.  T>  1  lere  the 
incomplete  observation  is  shown  by  tin-  dotted  line  cur\e. 
The  ext  rapolated  waxelorm  tlhe  c  ham  dotted  curve!  follows 
the  exact  curve  (solid  line)  closely.  In  many  applications  it 
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l ,  iit-r  li/urt-  n  deni"ii-t  rates  tl.*-  pcrlnriuitMi  t*  *»l  the  l , 
hft  in  in  t  l.i'  i  a-r  Hi  r<«  tin-  1 1  r -!  and  tin-  111  1 1.  harne-nii  ^  art* 
sampled  with  I  1  "J  tin-  [m-'-tlM'  :tj  -ample-  '1  iu*  direr!  ap 
pin  i«  «|i  »it  l  lu-  1*  IT  |u>.\  nl<-  t  hr  Muffed  -pr<  If  mi  indicated 
l,v  the  I  ii.un  1 1. .11  n  1  i  ur  .  e  Thr  r-miiaird  ITurier  -pet  truth 
I-,  |t  pit  -t  iitrd  l'1'  1 1  ir  ’•■Td  hlir.  .u%t  i ! . .  lli*-  r-timate  thr 

I  \.,i  t  -pn  i ruin  -A  *  1* *-n>  th.ll  tin*  t  urxe-  are  mdi-tinnuish- 
.,1  »ir 
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hi^  K\t r.ipi Til 1 1 »n  ul  tm* 1  i.  ;  t  I  l«  li.nn  dulled  »  urv*  »  in  tin*  -pal  ial 
<|i.fiuin  <  i.rre-pundmi:  U»  lrnpu-mx  donum  j »1  *  t  ..I  1*  t  I  In 
di-tud  d.i'tird  iiirsc  indii.iU-  tin-  mioinplete  « ■  1 » - « - r *• . 1 1 1 * < n  ss  th 
inea-ure  nn-nt  error  1  hr  t*xt  r;ip*  ■  hit  *-«  i  »ursr  i»  !..nn  dotted  Inn* 
j i *1  It >ss s  u.t*  ,ui w.d  \s,i\rli*riu  i-olnl  lim*  t  nrsr i  i  m-ei\ 
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SUMMARY  AND  CONCLUSIONS 

A  i u as  in*  thod  "1  rxtr.ijnd.it  ;n;:  ii  DI .  uu;i^r  tr«*Mi  n  measured 
-.pal  lt,lj\  Hum  ,di  d  m,.i:*r  i-  pre-ei,Pd.  Tin*  method  dim  tly 
.nidi  r--r-  tin-  l  — m  ul  tin*  til.iir  di  UN  r  nl  lrrrd<*lli  ut  the  UL 
J  Mini!.**  1  tn  tea -line  linage  r-l  itli.il  e-  iifr  t  "11 -Hat fled  to  ha  Vi* 
,tu  upper  ti'.uinl  **i.  tin-  numh.r  "t  Iretpienty  components 
i  "iiipriMi..;  tin-  umc.r.  A  ueometru.  nppn*.uh  ha>  hrt  n  takni, 
that  i-.  x.irmu-  i-rr^r  inr.i>urr-  wrrr  "plmu/ed  -uhirrl  tu  I'iven 
i"U-tr.u!ii-  I  Tin  *  rr**r  ii"rm- ’Arp’ du*  u-n-d  t  fu*  (\,  /.•.  and 

;  :  m-rti.-  1  lu-  ii. r’. In *< i  permit-  the  ndditmn  nt  "thrr  linear 
i "U-t r mi d -  < >n  the  dt  -.reti  linage.  That  nu\  mequalitv 
a  Imenr  i  *  nihin.il j**?i  "J  Hit*  ujikimun  samples  on 
"Hr  -nle  "1  tin*  ineiju.il it x  .uni  a  -  i  .d.ir  \nlue  "ii  t hr  "thrr  side 
i  an  ai-*  hi  .jddrd  1  In*  Hu  hi-iull  "!  o-ll-t Taint-  ot  this,  type 
ha-  l»rrn  I  **u  j  ul  l"  -lahdi/r  thr  "thrrv\i-e  uii-t. ddr  prnhlem. 
Nun  a  in  a)  rt  -ul:  -  dejm  *u-t ;  ate  tin  nn  rea-e  in  resolution  as 
w,  li  a-  it-  strrie.dh  a. '..mot  mr.t-ui *  iia-nt  noise.  In  terms  ul 
re-"».iti"ii,  tin* ! ,  norm  ua-  loumi  t"  ^ive  eonsi-tently  the  best 
spr*  t  r.d  e-i  jmates. 
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Two-dimensional  optical  filtering  of  1-D  signals 
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I.  Introduction 

In  t lit*  area  of  optics  dealing  with  tin*  optical  pro 
cessing  of  1  I )  signals  such  as  speech,  radar,  sonar,  or 
general  communication  signals,  the  2-1 )  nature  of  the 
optical  system  has  been  viewed  primarily  as  the  means 
for  providing  the  potential  for  greatly  increased  infor¬ 
mation  throughput.1  T  homas'  has  shown,  for  ex¬ 
ample,  how  a  coherent  optical  system  with  a  large 
space  bandwidth  product  in  one  dimension  can  he  used 
to  perform  spectrum  analysis  of  a  I  D  signal  with  a 
time  bandwidth  product  of  the  order  of  a  million. 
More  recently,  a  large  number  of  space  variant  coherent 
optical  processing  operations,  performed  on  a  1  1)  sig¬ 
nal,  have  been  added  to  the  repertoire,  such  as  fre¬ 
quency-variant  spectral  filtering.1  variable  spatial 
magnifications,  and  geometrical  distortions, ,c;  as  well 
as  various  1  l)  signal  transforsm  such  as  Mellin,7-*  Abel,s 
and  Laplace.  An  excellent  tutorial  on  space-variant 
coherent  optical  processing  is  given  by  Walkup.-* 

Coherent  optical  processors  can  also  be  used  to 
evaluate  attributes  of  two  identical  nr  diflercnt  I  D 
signals  such  as  convolution  or  correlation  between  the 
two  signals.  A  particular  signal  representation,  tin* 
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so  called  radar  ambiguity  function  tAFl,10  1:*  uses  the 
whole  output  plane  of  a  inherent  Fourier  processing 
system.  Mere  one  axis  represents  a  time  delay,  the 
range  in  format  ion.  while  t  he  ot  her  axis  represents  t  he 
I  )oppler  frequency  shilt,  a  measure  ot  the  target  mot  ion. 
Allot  her  interpret  at  ion  of  l  he  inlormat  ion  displayed  on 
I  he  I  wo  axes  is  i  hat  one  axis  represent  s  I  he  cor  re  I  at  ion 
while  the  other  axis  represents  the  energy  or  power 
sped  rum  of  the  1  D  signal.  Recently,  another  signal 
representation  was  introduced  in  optics,  the  Wigner 
dist  ribut ion  turn  (ion  ( WDl. 1  1  1  ’  T  he  \V1 )  of  a  single 
1-1)  signal  represents  the  spatial  variable  on  one  axis 
and  the  energy  or  power  spectrum  on  the  other  axis. 
Another  interpretation  of  the  \YD,  which  is  similar  to 
the  AF,  is  that  the  displayed  information  represents 
convolution  on  one  axis  and  the  energy  or  power  spec¬ 
trum  on  the  other  axis.  The  twosignai  representations 
are  not  independent.  In  fact .  hot h  t he  A F  and  t  he  \\  I) 
are  special  cases  of  a  general  class  of  represent  at  inns  of 
a  I  D  signal. 1,1  In  both  eases  the  I  - 1 )  signal  is  repre¬ 
sented  in  a  2  1)  space.  The  two  independent  coordi¬ 
nates  represent  important  physical  information  about 
the  1-1)  signal. 

The  two  spatial  axes,  coherent  Fourier  optical  im- 
plementat ions  of  the  AF  display  and  their  modification, 
lead  easily  to  2-D  spatial  filtering  of  1  1)  signals.  T  he 
purpose  of  I  his  paper  is  t  »  "xplnre  various  2  D  opt  ical 
filtering  schemes  for  l-i.  signals.  Such  signals  have 
applications  in  system  modeling,  bandwidth  compres¬ 
sion  of  speech  or  video  signals  as  well  as  in  t  he  general 
understanding  of  the  fundamental  limitation  for 
lime-frequency  represent.*!! ions. 
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II.  Formulation 

Tin*  shirt  inf*  point  for  -1)  filtering  of  I  D  signals  is 
I  hi*  generalized  span*  frequency  (('.SI*')  reprcMiitat  ion 
of  1-D  signals  introduced  hy  Cohen. *“  This  signal  class 
is  defined  by  the  triple  integral 

( ’!  i  .n  .// 1  =  1  jjt  |  |  |  <*xp[ni  x  -  .*u  -  Mi t| 

x  iiiu ( I  * 

where 

tilt,;)  -  fix  4-  :  'J \f\  \  -  i I 

where  fix)  is  the  11)  signal,  the  *  represents  complex 
conjugate,  and  Ilir.z)  is  an  arbitrary  kernel  function. 
Particular  kernel  functions  lead  to  different  space 
frequency  representations  of  fix).  Two  of  tin*  most 
important  ( JSF  representations  of  the  I  D  signal  are  1  he 
AF  and  the  WD. 

The  AF  is  defined  as 

=  j  expi-/(  K  )c<  i .?  ></r  gO 

An  alternate  representation  of  the  A  Fean  be  expressed 
in  terms  of  t  he  Fourier  transform  of  fix  ): 

F\u  )  -  |  fXpi-Mi  »  \J\  1  )ih  14  > 

The  AF  can  now  be  expressed  in  terms  of  the  Fourier 
spectra  as 

A  U\:  )  =  ’  jJT  j  VX])ijU  J  l(  r  if  .n\<ltl  ,  (f‘< 

where 

(Hr.u  )  =  h'iw  +  i  :2)F*itr  -  r/2>.  MU 

The  WD,  another  (JSF  representation,  is  defined  as 

n  U'.m  t  -  |  «*xp<  —  ju •*  i  ,j i<^  (74 

Similar  to  the  AF,  the  WD  can  be  found  from  the  Fou¬ 
rier  spectra  as 

H'(jr,u')  =  export  W,*(r.rrl</e  (Si 

The  two  representations  are  not  independent.  In  fact 
the  AF  and  t  he  VVDare  related  by  a  scaled  2-1)  Fourier 
transform: 

/t(r,z1  -  D rr  j  |  t'xpj— jit'x  -  tt  z  )|  IV  (  t  ,n  \tiuiit>  (JU 

The  WD  of  a  I  D  signal  can  be  interpreted  as  the  dis¬ 
tribution  of  signal  energy  over  space  and  spatial  fre- 
quency,  d  ims  a  2-1)  display  of  the  WD  represents  the 
energy  content  of  the  waveform  localized  in  space  and 
spatial  frequency.  However,  t bis  interpretat ion  is  not 
always  true  since  there  are  eases  where  the  Wl)  is  not 
positive.  One  of  the  most  important  properties  of  the 
WD  is  that  a  shift  in  either  space  or  spat ial  frequency 
of  the  1-D  signal  leads  to  a  proportional  shift  in  the 
corresponding  WD.  This  is  not  true  lor  the  AF  where 
such  shifts  lead  to  an  additional  phase  factor.1 ;  In  fact. 


the  magnitude  of  the  AF  is  completely  insensitive  to 
shifts  in  space  or  spatial  frequency.  I  detailed  proper! 
of  both  the  AF  and  WD  cm  be  found  in  tin*  litera¬ 
ture.1  *-,H 

Hot h  t  be  AF  and  1  he  Wl )  art*  special  cases  of  the  ( 1SF 
distribution  ('.  In  fact,  I  be  ( JSF  fund  ion  f  *  ran  lie  ex¬ 
pressed  either  in  terms  of  the  AF  as 

Cl  \  Jt  Jh  ~  1  ...t  ||  »  xp1  Hi  *  -  (.4  n  tsV  a.v  MO) 

or  in  terms  ol  t  he  \\  I )  as 

Cl  t  .o  .//  I  -  1 .  .t  ||  hi  i  -  ;  .tt  -  i  )  IV  i  *  ,i  !</*</,  ,  (]|| 

where 

/tH.rrl  -  |  |  e\p{Mrt  -  f.  1 1 // » f  hit  *1:  M2t 

is  a  scaled  2  D  F< airier  transform  of  tlie  kernel  function 
//  'Hie  kerne)  function  h  i.v./c)  represents  a  2D  im¬ 
pulse  response  filter  acting  <m  the  space-frequency  de¬ 
pendent  WD.  Since  the  Wl )  and  the  AF  are  related  by 
2  1 )  Fourier  Iran*- format  ion,  the  kernel  function  //  acts 
as  a  2-1)  filter  on  t  he  A F. 

An  all-pass  filtering  function 

// 1  =  l  <  t:u 

t ransforms  t  In*  1 1SF  funct ion  ( 1  in 1 1 1  t be  WD.  An  im¬ 
pulsive  kernel  funct  ion 

Hi:  I  -  2**1  x  -  :  »«>(n  -cl  M  l) 

transforms  tin*  (JSF  function  C  into  the  AF.  Other 
tiller  functions  have  appeared  in  the  literature.  For 
example,  tin*  all  pass  filter 

Hi:  >  -  rxpld ’.•((  .*) |  (1M 

transforms  the  ( JSF  distribution  C  into  an  energy  d is- 
Iribution  function,  This  lifter  function  was  first  pro¬ 
posed  hy  Ribac/ek.1'*  The  energy  distribution  function 
has  properties  similar  to  tin*  Wl).  Another  filter 
function  is  a  (Jaussian  kernel 

Hu  )  =  i*xj»|  -n  n  *  f  *  h-\)  (1  f>) 

I di i s  function  lias  been  suggested  by  deBrui jiv1' as  an 
aid  in  the  interpretat  ion  of  the  WD.  a  non  negative 
energy  distribution.  It  can  be  shown  that,  by  choosing 
l hi*  parameters  u  and  b  in  a  particular  way,  the  modified 
WD  will  always  be  real  and  non-negative.  Further,  this 
type  of  (JSF  distribution  is  then  consistent  with 
Heisenberg’s  uncertainty  principle. 

lit.  Two-Dimensional  Filtering  of  1-D  Signals 

The  (JSF  distribution  can  be  interpreted  therefore 
either  as  a  generalized  Wl)  (( JWD)  or  as  a  generalized 
AF  t(JAF).  While  the  present  discussion  will  deal  with 
t  he  ( JWD.  wit  hunt  any  loss  of  generality  it  can  be  car¬ 
ried  over  to  a  discussion  of  the  («AF.  Particular  fil 
tering  of  the  ID  signal  fix)  leads  to  a  new  signal  hi.  v  ). 
The  Wl)  of  the  new  signal  is  related  to  the  Wl )  of  fix) 
through  a  linear  2-1)  filter.  This  (JSF  filter  is  the  WD 
of  tin*  1-Dl  ilter  funct  ion. 
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For  example,  the  1-1)  linear  spare  invariant  filtering 
of /(  (  I  is  represented  as 

tul  =  /ul4/lt)=  |  fw\l it  (171 

where  / 1 1 )  is  the  spatial  impulse  response.  The  \VI)  ol 
I  In*  tillered  signal  k  t  \  )  is 

t *  =  H\ 1 1 .11  »  -  j  H’a \  ,ii  )U  a  i  -  \  .u  )</\ .  i  IS) 

where  U  ^  is  t  iu*  Wl )  of  the  tiller  lunct  ion.  ( 'omparing 
Ki|.  ill)  with  I  he  present  equat  ion  we  note  that  I  he  now 
Wl )  is  t  he  ( ISF.  The  new  \\  I )  may  also  be  eoiisidered 
as  a  generalization  of  the  Wl)  of /(.v).  Filtering  the 
signal  in  i  he  spare  domain  also  leads  to  filtering  of  t  he 
W! )  in  t he  spare  domain. 

Multiplication  of  the  signal  fix)  in  the  spare  do¬ 
main 

A(x )  -  f[x  Pm u  i  (ip) 

leads  to  1 1  Wl): 


•  ! .  t  T  f<  ~  ii~  -1 

Fig.  1  liask- 1  «ifu‘rt‘itl  Mptit  al  .irr.iii^cuu  nl  to  ^(‘Merrill*  tlu*  AKand 

1  iu*  in*  «t  it  hit  \\  I  *  In  ftl.'itir  /  ’|,  t  to*  J  t )  si^n.il  it  rt*  pl.u  **»t  m*xt 

t*i  f.»*  li « *(  1  it  r  vulli  l  hr  ni.i^k^  turiu*<i  pu°  ;ip,irt  Tlu*  ni.isk  jissrinhly 
o  r**l.ilnl  pc  nl.ilnr  i«»  (lu*  Iriirt  symmetry  axis.  In  plant*  !*•  the 
ina;:niluil»-  « *1  llu  i*  «l  Wl)  is  generated,  vslulf  in  plant*  /*  i  l lie 

majlinttale  t»l  tlu*  At*'  is  feller. »led. 


Therefore,  the  spectrogram  is  also  a  member  of  the  (»SF. 
In  this  rase,  a  fully  2- 1)  fiber  .smooths  the  Wl)  of  the 
signal  fix). 


'This  equation  represents  1  -I )  filtering  of  the  Wl)  in  the 
spatial  frequency  domain.  Here  again  IT,,,  is  the  Wl ) 
of  the  mult iplieation  function  m(r).  This  is  another 
( i \\  1 )  of  t  he  signal  fix  ). 

Another  type  ol  tillering  ran  he  obtained  hy  weighing 
the  signal  with  a  sliding  window  function 

Alt )  =  /{  .%  Pm  I  t  -  y )  121) 

The  TiWl)  of  the  signal  k ( v  )  is 

U\u.u  )  =  1  ;n  j  VVfU./HV„,u  -  y.u  -  tuit  (22) 

If  we  consider  values  of  the  («WI )  only  on  the  .v  =  y  axis, 
we  obtain  a  new  version  <JW|)  called  a  pseudo- W I ). 
"There  are  other  t  iller  opera! ions  t  hat  can  he  performed 
on  either  fix)  or  on  F[w)  that  will  lead  to  (JWI).  Bv 
defining  the  filtering  operations  on  the  2  1)  Fourier 
transformed  Wl ),  a  ( *AF  ran  also  he  lound.  I  )epending 
on  tlu*  type  of  information  needed,  spare  spatial  fre¬ 
quency  or  space  spatial  frequency  delays,  either  CAYl) 
or  (iAF  ran  he  generated. 

Another  2-1)  display  of  a  1  -  I)  signal  is  the  short -spare 
Fourier  transform  (SFT).-1  'This  transform  is  given 
as 

(/U.u  )  8  |  t*xp(-jn  \  )A(i  ,\  U/v,  l2:t> 


IV.  Optical  Realizations  of  2-D  Filtering  of  1-D 
Signals 

A  starting  point  for  the  optical  realizations  of  the  2-D 
filtering  of  1  !)  signals  is  the  large  number  of  optical 
realizations  for  lhr  AF.  'There  art*  coherent  space - 
integrating,1’1*11  coherent  time-integrating,—  incoherent 
lime-integrating/-’ 1  coherent  hybrid  (time  and  space 
integrating),-1  coherent  joint  Fourier  transform  ho¬ 
lography,^’  as  well  as  other  methods'-’*’  -*1  available  to 
optically  display  the  AF.  While  all  the  above  imple¬ 
mentations  ran  be  modified  to  be  used  as  a  starting 
point  tor  the  2-1)  filtering  of  l  -l)  signals,  here  we  discuss 
some  experimental  space- integrating  coherent  optical 
implementations  of  these  filtering  operations. 

The  basic  coherent  optical  realization  of  the  Wl)  is 
shown  in  Fig.  1.  'The  field  in  plane  I\  is 

u u  ,\ )  -  / 1 t  / \  2u  +  y  )\f  *  1 J  \  2u  -  y )) -  (2t>) 

In  plane  T  we  have 

o<  i  ,\ )  -  j  j  g U ' ,v ' )  expL/fc /2/U X ’  )| 

x  v\\)\-jh/'lfiy  -  x’)‘\dx’dy (27) 

where  k  is  tin1  wave  number  and  f  is  the  focal  length  of 
l he  cylindrical  lens.  I  )ue  to  the  cylindrical  lens,  we  have 
a  Fourier  transform  in  the  v  direction  and  a  Fresnel 
kernel  transform  in  tin*  other  direction.  'The  above 
integral  ran  he  expressed  as 


where  k(x,y)  is  the  expression  given  hy  Kq.  121).  The 
spectrogram  .V(r ,/r )  is  obtained  by  taking  the  square  of 
the  magnitude  of  the  SFT  for  all  possible  window  po¬ 
sitions: 


(»U.\)=  |  |  llu  uM’Xi  -  i  ’}  t*X|>|  - /A  ■  2/1  v  -  \  \'\ilx  </\  ’ 

(28) 


Spc.u  )  *  |</U ,n  )|-  (21) 

1  lowever,  the  spectrogram  ran  also  be  expressed  as 

S  l  x  ,i  r )  =  V.’ir  (  (  Wf  t  V  .2  \  tv.„  ( v  (I  ul\<h  l2'») 


In  plane  a  filtered  version  of  the  Wl)  is  generated. 
Using  a  spherical  lens,  the  convolution  operation  can 
be  mapped  into  a  product  ot  Fourier  t  rauslorms  Tims 
in  plane  I*  \  we  have  the  AF  multiplied  by  a  quadratic 
phase  factor  'The  magnitude  of  the  product  j.s  the 
magnitude  ot  the  AF.  The  phase  factor  ran  he  elimi- 
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taining  five  different  spatial  f« 


41) 


input  is  a  phnso  mask.  The  transilluminaUui  mask  is 
in  turn  imaged  using  twu  spherical  lenses  to  plane  /' 

[n  plane  /'.i tin*  magnitude  oft  In*  energy  distribution  is 
recorded.  Figure  *1  shows  the  energy  distribution  of  a 
signal  with  live  different  spatial  frequencies. 

'There  are  a  number  of  coherent  opt  icnl  filtering  op 
orations  that  can  he  performed  on  the  1  1)  signal.  For 
example,  using  the  same  arrangement  as  in  Fig.  l.with 
t he  exception  of  placing  a  vertical  slit  on  the  two  masks 
in  the  input  plane,  another  type  of  filtering  operation 
can  he  realized.  Figure  5  shows  the  results  of  (his  fil¬ 
tering  operation  on  two  and  live  spatial  frequencies, 
with  and  without  tic  filtering  in  the  Wl)  plane.  'The 
results  depict  the  magnitude  of  the  modified  Wl).  The 
mudifieat ion  is  the  additional  Fresnel  tillering  due  to 
the  cylindrical  lens.  Note  the  much  reduced  effect  of 
the  sidelohes. 

Figure  ti  depicts  the  coherent  optical  implementation 
of  the  SFT.  'There  are  a  number  of  window  functions 
discussed  in  the  literature.’10  In  our  experiment  we  used 
a  uniform  gate,  approximated  hy  a  45°  slit,  as  the  win¬ 
dow  function.  'The  running  filter  can  also  he  placed  in 
the  Fourier  plane  to  affect  the  SFT.  Placing  (he  run¬ 
ning  filter  in  the  Fourier  optical  plane  of  the  1  -1)  signal 
may  he  advantageous  if  we  wish  to  realize  a  real-time 
acoustoopt  ic  (AO)  implementation  of  the  SFT.  Since 
there  are  physical  limitations  on  the  size  of  the  active 
area  that  the  present  AO  transducers  will  support,  the 
•If)0  slit  will  degrade  the  performance  of  the  SFT. 
Placing  the  slit  in  the  Fourier  plane  removes  some  of  the 
burden  on  the  AO  modulator.  'The  magnitude  square 
of  the  SF  T  is  the  spectrogram.  Hy  placing  a  square  law 
device  in  the  SF  T  plane,  sik  h  as  photographic*  film  or 
a  gamma-corrected  video  camera,  the  spectrogram  of 
the  1-0  signal  is  generated.  Figure  7(a)  shows  the 
spectrogram  of  a  double-sided  frequency  ramp.  The 
spectrogram  with  an  impulsive  running  filter  function 
has  been  termed  as  the'  local  spectra1  ’  of  the  l-l )  signal. 
It  is  possible  to  generate  space  delay  and  spatial  fre¬ 
quency  delay  spectrograms  hy  magnitude  squaring  the 
*2-1)  Fourier  transform  of  the  SF  T.  Figure  7(h)  shows 
the  delay  spectrogram  of  the  double-sided  frequency 
ranq>.  For  an  impulsive  running  filter  function  this 
spectra  may  he  termed  as  the  local  Doppler  spectra. 
The  local  Doppler  spectra  does  not  have  the  many 
sidelohes  characteristic  of  the  global  Doppler  spectrum. 
Figure  8  shows  the  local  frequencvnnd  the  local  Doppler 
frequency  spectra  of  five  different  spatial  frequency 
quasi -periodic  signals.  Note  the  excellent  frequency 
resolution  of  both  local  spectra. 

V.  Summary  and  Conclusions 

A  new  signal  representation  lias  been  introduced  that 
allows  Mu*  simultaneous  space  and  spatial  frequency 
filtering  of  1-1)  sij  aals.  'The  new  representation,  de¬ 
pending  on  the  particular  filter  function,  can  represent 
the  Wl),  the  AF,  various  filtered  Wl)  and  quasi-WD 
functions  as  well  as  different  spectral  representations 
such  as  the  spectrogram,  the  local  spectra,  and  the  local 
Doppler  spectra.  Using  the  2-1)  nature  of  optical  signal 
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processors  all  the  above  (iSF  can  he  displayed  as  well 
as  manipulated.  Kxperimontal  results,  using  coherent 
space-integrating  optical  processors,  have  been  pre¬ 
sented.  Other  types  ot  optical  processor  can  also  In* 
modified  to  display  I  he  l  ISF.  These  processors  \\  ill  find 
applications  where  the  real-time  display  ot  time-fre¬ 
quency  information  is  of  interest. 

'This  work  was  supported  in  part  under  a  grant  from 
the  Air  Force  Office  of  Scientific  Research  AFOSR 
SI  -Oliil)  and  a  contract  from  the  Rome  Air  Development 
(  enter  FllH  128-80  (MlOKfi. 
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AB6TRAT 

In  this  pojer,  omilyzo  a  l  an,  Liu  TW:  A  access 
protocol  foi  sat  cl  1  i to  channels.  Tin*  analysis  is 
based  on  obtaining  the  l  ilt*  and  busy  jn  nd  dis- 
tribal  Km.  Tin.*  average  j«K’k'’t  gueuurg  d*.  lay  in 
the  user's  but  for  is  obtained.  A  rtdiluxl  iano«m 
TWA  is  .suggested  to  inif  tove  the  (oi.  im.-hv.  p:- 
plof.ojitat  ion  of  the  schem*  in  r.mlt  i-sp  »t  W.msnto- 
lilies  is  pios«*nt(\l  and  the  delay  on  Lon  id  the 
satellite  is  obtained. 

1.  mm  mi  it  i  itj 

In  [ticket  switched  satellite  cxjunumcat  ions 
systems,  various  multiaccess  r.chn**s  hive  b»vn 
proposed  in  order  to  accomplish  flexible  cxxuiv*cti- 
vi ty  and  efficient  sliate  of  satellite  channels. 
Those  channel  utilization  proUx\:>ls  nay  L«*  | ** i l t  i  — 
tioned  into  Uitec  mini  categories.  The  fust  is 
tlie  conventional  protocol y,  such  its  tuqueney 
division  milt iplexing  access  VWA,  and  time  divi¬ 
sion  Multiple  access,  TWA.  These  ate  effective 
for  heavy  traffic,  but  with  a  bursty  nature,  such 
a  fix<d  allocation  of  cliannol  capicity  is  extreiii*- 
ly  wasteful.  lor  such  a  situation,  larnkm  access 
schoik’S,  like  the  AIOIIA  protocol ,  tin*  Uni  scheme, 
an<1  the  Collision  Ri'soluti on  Algol  ith«;s  ITiee  Al¬ 
gorithm)  ,  are  mt  induced .  The  third  category  is 
the  dynamic  assignment  protocols.  These  include 
the  polling  and  reservation  schtines. 

In  Lius  pip.T  wv  analyze  a  landun  TWA  nul- 
tiplc  access  protocol.  The  analysis  is  Uisid  on 
obtaining  the  busy  and  idle  p'rinl  distributions 
similar  to  the  approach  use!  in  (l AM  711)  in  the 
ana  lysis  of  TWA.  The  ntdifnd  land  in  TWA  is 
tlien  projofied  as  an  liuptovemrnt  to  the*  land- «n 
TWA.  Other  randomized  TWA  piotctols  can  l*e 
found,  for  cx*uiplc  in  (KU;i  78 J  and  IKPUH  80|  . 
Implementation  of  the  scheme  for  multi-spit  Lx\im 
satellites  13  considered. 

A  priiuuy  concern  m  satellite  cMiiitnnic.it  u nr. 
la  to  rover  a  wide  area  with  u  high  gun  sp>t  hum 
antcnrui.  A  multiple  spat  L-e.un  antenna  his  a  sig¬ 
nificant  cl  li*ct  in  nuking  the  use  of  avail.Ude 
frequency  hinds  and  satellite  i  cannier*.  nno 
efficiently.  If  we  consider  a  sysl**;?  ufo.v’  an- 
tonivi  is  arrangid  to  fonn  N  non -over  lapping  sp»t 
beams,  each  serving  a  sej urate  zone.  IVich  h  am 
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nvvrisi's  an  uplink  ami  a  d*  will  ink  at  different 
carrier  1 1 o juencies,  with  evun-x't  ions  rad**  h’tweest 
up  ami  divnlwiks  thiough  a  : w itch  nutnxWn  Uie 
satellite,  slxv.ii  onceptuully  in  Pig.  1.'  The 
•witch  Matrix  int  eiconnects  leceiwis  and  trans¬ 
mit  tern.  acootding  to  a  \i  v-\  i  erenh  d  algoi  1 1 1  cn . 

I'oi  sinplicity,  c»rpb*te  ttaltic  spnetiy  is  assum- 
t\i:  each  z» -r.e  contains  the  sire  nunh*i  of  stations 

b,  wtuch  is  j.!kvii  by  aii;  i  =  l, . ,?1  anti  p-1, 

. ,h.  huh  stalit'n  tians.uts  pn>.ets  to  tile 

sia'  n.nb't  i  i  st.it  i-  ns  Nl.  and  all  interst at  K  n 
pi  'M.'t  lates  aie  i^jual.  Th**  ariival  p'i^o.*ss  is 
husson  aixl  all  n«  :a»iges  consist  of  a  suvjle  pic¬ 
ket  .  All  bulh*r  capicitios  are  infuute. 

In  (ClAIiK  801,  tlie  authoi  s  ar.ipn o  different 
piot^cols  ftn  aciX'ssing  a  Mult  i-lo.un  .satellite. 

I  .'r  exar.pl e,  U.  t  hi*  I'lXixl  A:  t,i<jnju*nt  TWA;  th.e 

I I  am*  sliui’luie  is  sh  v.n  in  ITg.  I'.irh  *-aitli 
station  d  1 1  t*ct  s  its  tint  tic  to  the  desiitvl  des- 
tir-it  ion  zoi.e  l  y  1 1  an*<MUt  t  ing  into  the  peassiuiuvl 
rb  t  dvuir.c  wtitch  the  switch  is  set  Id  tii.it  d»-sti- 
nat  ion.  The  iwitcinng  smjui  nci*  in  an  vi.-Uuid 
Udlth  luid  ix  bdliw:.  a  J  I  <  ai  :  u; ,  *d  sc) ^ \\ \ b»,  p*r.- 
rddy  mxhlul  at  mi; oqnent  inteiwils  ly  cciitiml 

50  dials  turn  tl.e  gtouivl.  !\Uo  slottivl  AIOHA/TIM 
sy*’.t,ii*t;  with  rndtiple  uj'liiiks  h*iv«*  h*cn  investiga¬ 
te!  toi  broadcast  mg  satelliti*s  in  ISIV>A  80], 

U'AVi:  79]  and  (bh^i  70]  . 

In  sect  ic»n  2.  wv  Jej.ci  ibethe  i  and, in  Tlf-A,  and 
th*.  n  aiuilyz»*  it  in  mvtirn  In  m  ct  n*n  \  we  dis- 
cii: «s  the  inph Mentation  of  tht.*  r.clurt*  in  nulti-spT 
1*mmi  rat  el  1  it  es. 

The  hukliin  TWA  lTot»vol 

H> *i e r  we  will  dc:.ci  ih*  the  hmdiin  TW*tA  pioto- 
i\d  let  a  gloUil  Iv.un  covet  ing  the  enliie  area  to 
h*  retvicixl. 

W*»  aii.*  conssii  i  in  i  M  kmI  us«*is.  Tin*  M 

i*  **is  aie  il l v nl \!  into  N  giimps.  |.»>Ti  uieup  i\*n- 
srl’i  of  P  us*  i  (i.e.,  M  ti|.)  .  luch  i  e  '  i  ■,  hr;  an 
intitule  bill  lei.  Messages  ai  rival:;  at  t  a,  h  u-i-i 
f’l  Itw  the  statistics  ol  .1  l\t  is:. on  pinv:  !'.  wit  It 
i.de  t  hich  n«'s:.ag*.*  is  a*.si:::»^!  to  O'nsisl  o|  a 
•;  l  ng  le"  picket  .  Chinn»*l  tir«*  is  slotted  an-1  all 
U'  in;  aie  synclii  win  /•*«)  to  1)n*  iiunin'1  t  Hi*.  The 
slot  length  is  ojual  to  tin*  duiatu  n  ol  ,  me  pic- 
ki  t.  Chinnel  slots  me  gn  up>l  in  ti.un  s  eacli  (d 
w*ii  h  consists  of  N  consi  out  ive  slots. 

hull  gi* 'Up  is  nssign«d  a  lixid  lavoioi  slot  in 
tin*  linin',  as  m  the  case  of  dedicated  slots  of 
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ordinary  'll MA ,  and  no  two  groups  Uiv»*  the  same 
favored  slot.  Mach  user  (or  earth  sl.it  ion) 
within  a  group  has  the  same  favored  slot. 

When  a  user  lias  a  jacket  ready  tor  trans- 
mission  (i.t.,  the  iiiiunt  a  jacket  arrives  at  the 
Hbnd-of-tho-I.ino  (liOI.)  position  in  th.e  burim),  ho 
looks  at  tin:  next  N  slots  and  schedules  tr  an.  .mission 
in  t he  f awrcJ  slot  of  hi s  g rouj »  with  pi ob. ih  1 1  r t y 
a  and  in  any  other  funicular  slut  witnin  the  next 
N  slots  with  prolxibi Uty  b  Thus; 

a  +  b  (N-l)  -  1  (1) 

lx?t  us  explain  tins  process  using  a  roulette, 
which  has  the  numbers  1  to  N.  One  of  these  numbers 
aji«ars  with  probability  a,  the  others  apj«vir  with 
probability  b.  IXiring  busy  jeiiod  when  there  are 
sene  jxickets  in  tiro  buffer,  as  sevn  as  a  jacket  at 
MOL  is  transmitted,  the  loulette  is  tuimd  anti  hence 
the  MOL  jvicket  transiiussi on  is  sebduled  according 
to  the  result  of  the  loulette.  During  idle  jariods 
when  there  is  no  jackets  in  the  butler,  as  s.*an  as 
a  jocko  t  arrives  at  the  buffer,  the  roulette  sche¬ 
dules  the  jacket  transmission ,  Fig.  3. 

If  a  jacket  encounters  a  collision,  which  is 
acknc*/ledg«tl  after  a  round-trip  tin*'  delay,  it 
joins  the  buffer  at  tiie  Knd-of- the- lane  (iX)h)  josi- 
tion  in  the  sano  way  as  a  jacket  rwvly  genoratml. 

The  queueing  ni  xie  I  with  feedback  is  shev/n  in  Fig. 4. 
The  message  arrival  rate  including  collided  pac¬ 
kets  is  denoted  by  X 

G 

3 .  Analyst  s 

let  Di  bo  tire  exjocted  jacket  delay  due  to 
the  user’s  buffer,  The  exjocted  packet  delay  in 
the  syr, t an  is  then  given  by 

D  =  2P  +  Dj  +  E  (D i  +  2P)  (2) 

wtiere 

E  -  tlie  average  number  of  retransnissions 

2P  =  round-trip  delay. 


F  “n  in  a  pu**f  similar  to  Uim's  jutxif,  Dj  is 

11  Vi  *11  by 


l,T  +  I‘(ir*(Fj)  ,6) 


After  scire  nkuiijAilation  sbxvn  in  Aji  endue  A  we 
obtain 

I»‘  |s)  =  (a-b) 2F!'T+|2ab*  (N-2)b2)  t  e'tl  Si> 

l=i 


(7) 


2  /  o"s(y*n1  dF(y)  + 


1!  is)  =  (a-b) 


N  -  1  ST  ) 

•  |2ab+(N-2)bJl  1  e  N  _G_ 
i=l  >,s 


sT 

~N 


v;r 


l-Q 


(8) 
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F(y): 


-1] 


1-c 


-XT — L  4  u(y"IHN) 


-  ll-if  N  "  yl)  0  T 


19) 


y  *T 


whom 


U  <x)  =  ^  1  x  :  .) 


0  x  •  0 


(10) 


Fum  eqs.  (6),  (7),  (B)  wo  c«m  easily  blaim  the 
f  1 1  st  ami  seivml  in.iit.nts  of  jacket  service  tines, 
b,.  b2,  6r  b,  . 

3.2  Tin*  Average  f'n’  U-r  of  pj't  i  aiv-i  ii^sioiLS  F 


3 •  1  Tire  Ave rage  Packo t  Waiti ng  T iito 

To  find  the  value  of  P ,  we  use  an  aj preach 
similar  to  the  one  used  by  Lam  (IAM  77]  in  his 
analysis  of  "DMA.  We  consider  two  kirvls  of  ser¬ 
vice  time  distribution  function  B(x)  and  B(x)  in 
a  single-server  queue  with  Poisson  arrivals  at 
jackets  jx'r  second  B(x)  is  the  distribution  of 
a’juckot  which  initiates  a  busy  period  witJi  first 
anti  second  nunents  bj  and  b^.  All  subsequent  jac¬ 
kets  in  the  same  busy  period  Imve  service  times 
drawn  independently  frein  the  distribution  B(x) 
with  first  and  sccoivi  nx  incuts  bj  ami  b~.  let  N 
b'  the  nunber  of  jackets  in  the  system  (both  inL 
cpieue  and  in  service)  at  tine  t.  We  now  define 
tlie  transforms: 

P(z)  =  i;  Zn  P  :  P  =  1  nn  Prob IN  -ii)  (3) 
n=0  n  n  f  -  L 


n  (o) 


/  e  0><  dB(x) 
o 


(4) 


CO 

D  (s)  -  y 
O 


e  d  6  (x) 


(S) 


Tlw*  average  nred-u:  <u  tet  i amavassions  icquir- 
(\1  j  jacket  ec'ner i  >  lelatovl  muier  ojuili- 
Luiur  conditions  to  the  system  t brought Mt  and 
G,  tlie  clumnel  traffic: 


1  +  E  =  |  HD 

where',  s  -\T  aul  G  -  V  .T.  Mere  a.,  the  ptoli- 
ainlity  tiuf  a  jacket  t  ran;j;ti  t  lid  in  a  favored 
slot  and  b  ,  tJn.*  j»tolability  tikit  a  jacket  is 
trani3mtt«d  in  a  m »n-t iiwie^l  slot,  are  given  by 
a  -  aV.T  and  b.  -  bv  .T.  l\!finin«j  Fr  to  1«'  tlie 
j>i olvibiVity  tlkd  a  jacket  is  ti.miaiuVuxl  succens- 
iully  and  P  to  b*  tit--  jiobibility  that  a  jacket 
o>llid*s  with  :  •«  nu  *  i  *thft  jackets,  and  assuMi.g  in- 
d«  -|  x  ‘iii  J<  nee  ,  average  nuiib-'i  s  of  i  i't  I  ansrd  si*  ion 
.it tempts  is  given  by 


E  =  1  P.1 

i-O  11  1 


r  iv 


•'-■V2 


(12) 


W’i  ic  re 


V 

s 


r  i  ob(  A  suecessd  ul 
Slot  l 

Pi ob l A  success Lul 
ml  slots;) 


a  ho  *  b  p 


20  ‘ 


tronsru  an  ion 


t  ransjiu  ss  ion 


in  a  lavoicxl 
in  i on- favor- 
(13) 
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rf  =  I’nb  (col  1  isk'h:;  in  a  faunal  '.loti 
+  rub  [col  lisions  m  non-favouxl  slots) 

Ml  M-l 

*  a  ;■  v  +  iw-i)  b  :  r  (U) 

i  1  1  i-1 

I>u  =  I’tob  (i  pockets  m  (L-l)  favor  ixi  and  (M-I.) 

non-tawuxl  pickets  aie  Uan:juitU\l  simul¬ 
taneously)  . 


]=0  i.-i  i  i  1  i  M-i  ‘-1 

=  au  a  (i,  I  i-l)  (lrd 

=  rtobi  i  pickets  ill  L  tavouxl  anil  (M-h-1) 

1  non- 1 avos t\i  packets  aie  tiansnutttxl  s iruul- 
turuvusly  1 


v  e>  a  3  M  _  ^  *  h- 1 -  ] 

j-oL  3  1  1  1 


M-L-lCi-jbl 


M-L- 1-1+] 


u-b,  r  *■ 

nan  ti,L) 


4.  Miilt  i -JVv\;u  Satellite 

lion-,  we  t\  nsi.ioi  the  iiq  li in  station  cl  i  an- 
ikiTi  TDiiAn  r.rult  i-*>p>t  In  Min  satellite.  The  system 
nixie i  is  the  s.un'  as  m  Knj.  1. 

podr  user  ii.ii  N  butti’ts  each  of  which  a  ties* 
(Olds  to  the  ihlb’ient  best  n;. it  ion.  huh  but  ter 
cloys  the  aO  ve  t.»  i it  n  i »  '  1  taiulrt  aece.  s  r.elo.  r. 

The  mi  in  ptt_v«dare  peidiar  to  the  rul  t  sys- 

tun  in  the  assigns*  sit  el  the  1  oi . y  s  let  .  Tin* 
favor  ixl  slot  is  assnjn«d  into  the  h  “  sh  t  m  a 
{ivuiv*  for  the  pickets  which  ate  1 1  at  asm  1 1  id  tiom 
uplink  A  to  ek vail  ink  A  ,  where 


h  =  f  in-1 


( m-i 

\  l+ri-N 


l  ♦  )  *N  ♦  1 


When  a  jacket  is  ackm  xvhd  n  d  to  cikx'imtei 
a  collision  after  a  certain  belay,  it  is  letians- 
mittixl  at  the  PUL  in  the  user's  butter  vis.  bes- 
ciuWxl  in  Six*.  3.  tmly  nen-tvl  1  i.lid  pickets  are 
accept  ixl  into  t)io  puv»"srni  satellite  with  butter 
ini]  c.i|  uln  1  i  t  y  anil  1 1  anum  t  t  *xl  to  the  appropriate 
downlink  on  the  Ursis  of  ’ll jM. 

In  tlus  case  the  expect  id  (ticket  belay  is 
given  by 

n  -  2r  +  d.  +  n.  +  i:  (n.  +  m  ( IB) 

M  J  *.  J 

wtiete  L)^  “  tji«*  exp 'Cl  oi  picket  belay  at  the 
satellite.  Tin'  il  l  f  1 1  •  i  *  ’lice  hatt  h|.  (d  is  l  h».* 
teim  P,  due  to  the  quell  I  n> }  bebiy  111  the  ».  ilollltO. 
Hit*  exfocl oil  nunt)>er  of  pickets  in  a  satellite 
buffer  fot  a  specific  destination  is  given  by, 

(see  Appendix  B) . 


a,MH-l)  h, 

N  a  -jr- . fc.-T  1  - 

Il-dj-fN-Db,] 


1,-  a\t.TL,  ti,  v,  b \  yn . 

e’.  aj plication  ot  hit  tie’s  toiriulu  gives  the  ex- 
:•  ‘  <d  picket  b<  lay  P,  *  n.m* 

1  erne  S  shws  the  analytical  and  sj::uiutn>n 

i  esili  s  h*i  the  average  picket  queueing  delay  in 

i  h* '  eaith  station  v»u  sus  '  ,  th**  cli.irj-.el  traffic 
rate.  Both  result:',  ivur.e  j!t  •  witJi  each  ether.  it 
u  also  uot»d  thit  tr.e  ett»v*t  c-i  a  t,  is  i-uill  at 
1  v.'  a  rat*  an!  quite  l.uqe  at  hiuliei  cl.uJU’.el 
trait  ic  rates  line  to  increase!  o*llisi<"*. 

injure  (•  s.h  xvs  th*.*  cuuilytical  and  sum  bit  n  n 
r  >  ‘-ills  lor  tire  total  av*»r  a  ;e  picket  belay  P  as 
upon  vu  e|.  (1),  versus  throughput.  l*ae  tluit 

t  :  hiabcL  p  (i.o.  close  to  IlflA)  wv  obtain  U  tter 
j «  a  U'i  rrvmce . 

Tl.e  ler.ults  can  U»  lrqiovi'xl  at  bwr  Uin.uqh- 
p.'  by  rrokiag  U'.i*  of  th**  t  act  tiot  the  nlle  p>- 
ir  »isf  (wlreu  ne  pi  kets  are  in  tire  user's  larllei), 
are  i  :■  i*i ;  i  i.  it  it  .it  bver  thaoucigut  aivi  1  rrsy  jern»is 
an*  itiranar.t  at  hi  jher  t  ]a  on  rhpirt  .  Theju  d  it  n  d 
ran  fir  ’ll t'J\  is  b«  ser  lU.xi  as  loihvs; 

Purina  nile  pa  ioi.,  tb.e  arrivixi  pickets  .ire 
a  .nnroi  U>  t ;  an:u”a  t  t  oi  t»*ien'  t }..  >  r:,.|  *,  ♦  the 
1:  u't  ,  wbach  bis  the  lavetoi  sb't  ^  the  last  slot. 
Mi  l  ly,  whin  a  picket  arrives,  at  l  J*  si  t  m  a 
t:  ir  i,\>ns  ist  ir.a  ef  n-sh'ts,  it  is  1 1  anntat  t»xl 
Ubc'en  ( i  ♦  1 1 1  * 1  a:.d  ntl1,  when*  t).e  h*atni(  oi  the 
hr,  lid  slet  is  t):e  jjth  v%  1 1  i  j  t\jaal  piob- 

abriity.  The  •.  lrulat  joi .  r«.ults  t-r  t:  •  aveiaae 
p‘-’t  quiuomq  bebry  are  :  !.  wi  in  pep  !>.  it 
'  b  ••c-  ,:q  r  ••veil  :;t  s  at  bv.ej  t  la  oi.  rJ-.j  ut  inpaoi  to 
’lWU  an:  rarulci  lirarc  t'  :-}j  v.;  1 1 .« •  sir. rula- 

tn  n  l  e  s  ■» ;  1 1  s  to:  the  belay  tia  euahj  \it  x  h-u  acte:  i  s- 
tic-.,  which  ah.  iLm;  th.i*  ad.vautuqe  of  tiio  nixh- 
f  :-d  r and  n  T.KA  at  h  vci  tiuouqhput. 

i\  n  ihix  A 

In  this  a:  p  niix  we  be:  ive  b.p;,  (7)  to  (g)  . 
hr:  *  ,  :  all  we  'j,:  ib<a  !  ic  relation  U'tW’een  two 
P  .  .  cut  rve  pi.  k.  t  t  :  .u. ecu  ssn  ns.  let  i  (1  b  N) 

a:.  ;  i  (b.ph)  n  rie.,.  ut  tin*  slot  addxes'.es  of 
the..*  o  ns< s art  ive  picket:..  Kh.en  w>'  asiirva*  a  lav- 
cr . d  slot  uddi»s>  with  #1,  iciHiiwinl  slot  abdres- 
.  with,  #.?  to  arq  wt*  obtain  the  piclvibility 

1  lot)  li,  j]  f  a-  ,  i  :  1 

^  .ib  /  i  1 ,  )/ 1  or  l/l ,  ,)=d  l All 

,  b'  ,  i/l ,  i /  l  . 

Par  lira  bu*o,  pried:.,  the  piclvibility  Ui.it  the 
rvuv  tine  is  kT,  N  is  given  by 


hb)  (i,  j  w  ( 1 1  k )  }  ,  i  ik  .M 


|h*(Nd)l',  /*(N  l)t>;  1 

l  ^  I )!>-,? J 


w  fx)  =  |  x 4N  ;  xlO 

<  x  ;  0‘X^N 

lx  fi  ;  N*  x 

Th.eietoi  e,  we  ,  *1.1  a  in 
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iUMx)  _ 


k-l  i-l 


Viob  li,  j-v  (i«k)IMx-  -‘tJ  )  (A-W 


'll--  Z-l  * as-  !  •  >t  n  ol  bp  (1*1)  is 


(x)  ia  a  <5  -timet ion. 
filiations  (4),  (Al)  and  (A4)  qive  (b)  . 

rkv  wv  pioceix.1  to  deal  with  the  onto  tlut  aie 
idle  poi  nxl  l*;  included  lvtw  i  n  lw>  eur.s«\i;i  ive 
jurket  t lansmissions.  In  this  case,  i  asi  )  ie- 
presont  the  slot  addiesses  toi  the  last  nansub- 
sion  of  tho  piecedinq  busy  (ornxi  ,unl  the  tnsi 
1 1 amsiussion  of  the  itl i*'  |vih»l.  H  ue  wo  in- 
t  reduce  a  varial  y,  which  is  Hoi  mol  to  U*  1 1  a  * 
period  between  the  first  annul  and  the  cotics- 
jondinq  locution  in  the  niiruil  I  x\t«’*a  of  the  last 
dejxiiture,  as  slkvn  in  Fiq.  4.  Since  the  variable 
lus  tho  same  nxunvnu  as  that  defined  in  (I AM  781, 
tho  distribution  t  unci  ion  is  cjivon  by  Dj.  pi). 
Here,  when  y  is  in  the  ivqi on 


-l-  ■'  T  i  v  < 

N  * 


1  Shi  N 


the  pioWitulity  that  the  service  t  ine  is  y  +  ~r 
is  qiven  by  ' 


Prob  ( i ,  j  =  w  Utk-1 )  } 


-  <a2  Z»  1-a^)  {i 


i'  .1  -a  iti.t  I ;  («<:>  11.!.  t  la  ‘  h  :  la  k-r?.:s  hi:,.- 

tias.-t  :  :•*  bp  i:iCir.»'  7<],  we  luve  ti**  ru  i  .  .t -qcnei  - 
at  i:i>]  1 1*!.‘ -t  i<  n 


OC.-I  :  (I-,  ) 


I  *  i  )  V  t  .* ) 


.1/  '  7=1  r  ‘S  1 


"hi  •:  eh  e,  t  ne  cxpvtid  niri  ■  *i  ci  :nssa;>*s  m  a 
s«tt  <  ‘1 1  it  •'  bv.t  *  or  h  u  a  :  ,*vi !  K‘  >i«  :i  uat  r  'll  is 
ijiven  by  N  -  o.u)  winch  ijivi-:.  ioju.  (I1)). 

"d;;” 


1  i  h*k-l  i  N 


Therefore  wo  obtain 

-j~  =  >•  I  l-ivibli.j  vi(i*k-i.)l 

k=2-h  x=l  d'  y^x-fc* 


r;~“  ~T«  x  c  T  ( A‘3 ) 

Equations  (5),{Al)  and  (AS)  yive  eqn.(’) 


In  this  ap(*’ndix  wo  derive  lip  (2J)  .  Here  we 
think  alowt  a  concent! ator  with  N  teiTninuls.  The 
concentrator  out  jut  tran:suis:»i  on  line  rate  is  one 
d*U  a  unit  j«rr  ^  sixt.  Each  tenranal  is  synchroni¬ 
zed  to  this  output  late.  One  ot  these  terminals 
feids  in  either  one  data  vuut  with  pu  liability 
a  -  i 'very  ^  sec,  t  he  others  with  pi  cl  lability  b> 
oveiy  A  sec.  In  this  case  the  probability  Hint  k 
{tickets  arrive  in  the  concentrator  small  tuneously 
is  ijiven  by 


( 1 -a  ^ )  (l-b2)‘ 


tl'^,H-lCkb2k‘1-b2'N'l'k 


*2  N-lCk-IbZ 


k-l  m  _  v 

b2  ,Hl2'  |.k.N 
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input  is  a  phase  mask,  The  transillummated  mask  is 
m  turn  imaged  iimiij:  h\u  spherical  trusts  to  plant*  /' .. 
In  plant-  /\  t lu-  magnitude  of  tht*  energy  distribution  is 
recorded.  Figure  l  shows  t lu*  energy  distribution  ut  a 
signal  with  li\t-  different  spatial  Iroqueinies. 

II  a- it-  a  if  a  number  of  coherent  opt  ical  t  die  ring  op 
orations  that  cun  be  performed  on  the  1  I )  signal.  For 
example,  using  t  In*  sinu*  arrungetne nt  a>  in  Fig.  1 ,  wit h 
llu*  exception  of  placing  a  \ut ion l  slit  on  tin*  two  masks 
in  tin*  input  plain*,  another  tvpf  »  a  tillering  operation 
tan  lu*  realized  Figure  a  shows  tin*  result".  ot  this  til 
tering  operation  on  two  ami  live  spatial  frequencies, 
with  anti  without  tit  liltermg  in  tin*  W  D  plant*.  I  hi* 
rt-sults  i  If  put  t  lit*  m.ipiit  Utlt*  ot  t  In*  moil  it  iftl  W  1  ).  Hu* 
moditit  at ioti  is  tin*  ailtlilional  Fre-aiel  tillering  thu*  to 
thf  cylindrical  h*ns.  Not t*  tin*  much  reduced  oiled  ot 
t  In*  sitlt* lobes. 

Figure  n  depicts  tin*  coherent  op'.ieul  implo mentation 
til  t  lu*  SFT.  Flu* rt*  art*  a  numhe r  of  w  imlow  functions 
discussed  in  llu*  litt-raluif . In  t>ur  experiment  wo  ust*tl 
a  uniform  gate,  approximated  hy  a  F~>°  slit ,  as  t  ho  w  in 
i low  t inn  t  ion.  Tin*  running  t iltor  oan  also  lu*  platYti  in 
tin*  Fotirior  piano  to  at  feet  t  ho  SF  T  Flaring  tin*  run 
mug  filt  or  in  t  ho  I  on  nor  opt  teal  piano  oil  In*  1  1 )  signal 
may  ho  advantageous  if  wo  wish  to  roah/o  a  real  time 
acoustoopt  ic  t  AO  *  in i pit* inont  at  ion  of  t  ho  SFT.  Sinoo 
thoro  art*  physical  limitations  on  tho  size  of  tho  aotivo 
aroa  that  tho  prosont  At )  t ransducers  will  support,  t ho 
lo°  slit  will  degrade  tin*  porforntanoo  of  tho  SFT. 
Placing  tht*  slit  in  t  ho  Fourier  plain*  removes  somo  of  t ho 
I*  ui It* n  on  t  ho  A (>  nuuiulator.  Tho  ma;:nitmlo  square 
ot  tin*  SF  T  is  tho  spectrogram.  Hy  placing  a  stpiaro  law 
device  in  tin*  SF  T  piano,  such  as  photographic  lilm  or 
a  gamma  o *rreded  v  ideo  camera.  tin*  sptvfn^ram  <>J 
tin*  1  D  signal  is  generated.  Figure  Tint  shows  tin* 
sped  regrum  of  a  tlouhlo  sided  frequence  ramp.  1  ho 
spool rug  rum  with  an  impulsive  running  fillor  function 
has  boon  tor  mod  as  t  ho  local  spoct  ral  ’  of  tho  1  l  >  signal. 
It  is  possihlo  to  generate  spaco  tlolay  and  spatial  fre¬ 
quency  delay  spectrograms  hy  magnitude  squaring  tin* 
2  H  Fourior  t  ransform  of  t  ho  SF  T.  Figure  Tthl  shows 
tho  dolay  spoct rogram  of  tin*  double  siilctl  frequency 
ramp.  For  an  impulsive  running  filter  function  this 
spectra  may  ho  termed  as  tlu*  local  Doppler  spectra. 
Tho  local  Doppler  spectra  does  not  have  tlu*  many 
sidelohes  characteristic  of  tho  global  Doppler  spectrum. 
Figure  S shows  the  local  frequency  and  the  local  Doppler 
frequency  spectra  of  five  different  spatial  frequency 
quasi  periodic  signals.  Note  tho  excellent  frequency 
resolution  of  both  local  spectra. 

V.  Summary  and  Conclusions 

A  new  signal  representation  has  boon  introduced  that 
allows  tlu*  simultaneous  spaco  and  spatial  frequency 
filtering  of  11)  signals.  'Flu*  now  representation,  de¬ 
pending  on  the  particular  filter  fund  ion.  can  represent 
tin*  \VD,  tho  AF,  various  filtered  \VD  and  (piasi-WD 
fund  ions  as  well  as  dil  to  rout  spool  ral  represent  a  l  ions 
such  as  tho  spectrogram.  tho  local  spectra,  and  tho  local 
Doppler  spectra.  I  Ting  the  2-1 >  nature  4 »t  optical  signal 


processors  all  tho  above  (IS1*  can  hi*  displayed  as  well 
as  manipulated.  Kxpomuont.il  losults,  lining  coherent 
"pact*  integrating  optical  protestors,  haw  lurn  pn* 
Sfiitod.  (  U  In*  r  typos  of  optn.d  processor  can  .d>o  la* 
iiu  »d it  n  il  it*  display  t  lu-  ( iSF.  Thest*  pri  <t  esst  w  ill  t  iial 
applications  whore  the  real-time  display  ol  tune-tre* 
quoin  v  inform. itioii  is  of  interest. 

'FI i is  work  w.s  supported  in  part  under  a  gr.int  from 
tho  Air  F* >rct *  Ollier  nt  Sen  nlilif  lu -e.Mi  h  Al  tiSH 
S1  1 1 1 1 P.  t  and  .i  t  to  it  rat  t  Irmn  tho  Koine  Air  Development 
t  Tutor  FlPtijs  sw  C  tuipf,. 
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r .  a  ..i  ;  iw  a  -  r'.  n  2  3  a  -  u  J 
J.-ioa 

ABSTlvVT 

In  this  coje:,  v.e  -uoly/e  a  lan.kiv.  'll/*. A  access 
protocol  lot  sat  fl 1  1 t e  channel';.  Tae  an  i  lysis  ;  s 
based  on  obtaining  t J> ■  1  lie  and  l;;.,y  jeied  dis¬ 
tribution.  The  aveia  ;e  ; - 1  * K . ’ c  .,.-ieu.  *i:.g  d*.  lay  m 
tile  user's  buffer  is  nU.uj,nl.  Arui.tmi  iu:«u»i:i 
TOIA  is  sug-jesud  to  imjievo  tin*  In- 

plerentatlon  of  t  h*  -■  schciin  r :  L  t  i  - 1  ; l  u-.m  sate¬ 
llites  is  [  resunttd  and  the  delay  1:1  Loan!  the 
satellite  is  obtained. 

1  *  Intwjoi  on  t« 

In  picket  switched  satellite  cxmnunicat iuis 
systems,  various  mult  1  access  schir.it s  hive  Loon 
proposed  in  order  to  ucccuiplish  flexible  conn* nati¬ 
vity  and  efficient  shiie  of  sat  el  lit*?  chunn«ds. 
These  channel  utilization  protiouls  r.uy  L-*  j.irti- 
tioned  into  tluee  rum  cativjei  ms.  The  tirst  j  s 
the  convent i or, a  1  protocols,  such  ns  iiiAjnmcy 
division  multiplexing  access  hK\,  and  t  nu*  divi- 
sion  multiple  access,  11 fA .  lhi*se  aie  cMivtive 
for  heavy  traffic,  but  with  a  bursty  nature,  such 
a  fixed  allocation  of  ehmnel  eanicily  is  cxitriti'- 
ly  wasteful.  Tor  such  a  situation,  lauLin  access 
schemes,  like  the  AIOilA  protocol,  the  Urn  schist.*, 
and  the  Collision  Hi ‘dilution  Algorithm  (Tree  Al¬ 
gorithm),  are  introduced.  The  third  category  is 
tlie  dynamo  assiguient.  protocols.  These  include 
the  [oiling  and  reservation  sehirv's. 

In  tins  [wgor  wv  analyze  a  laukr;  TO*1A  mul- 
tiple  access  protocol.  'Ihe  analysis  is  based  on 
obtaining  the  bu.sy  aid  idle  [<riid  distributions 
similar  to  the  approach  us«.d  in  [I AM  7b)  in  the 
analysis  of  ibMA.  The  mxlifitd  iand-»n  rlUA  is 
then  pmjosud  as  an  n::prover«  nt  to  the  rand*  in 
ItYlA.  Other  random  zed  IMA  [  tolocols  can  to 
found ,  f  o  r  exurg » 1  e  in  [  K I .F  I  7b]  and  [  hi 'I  111  GO)  . 
ImplfTi/ntat ion  of  tin?  scheme  for  mul t l-spat  beam 
satellites  is  consider'd. 

A  prinuiy  concern  in  satellite  o  maim  cat  i»  ns 
is  to  cuver  a  wide  area  with  a  high-gain  ‘.pot.  ban 
antenna.  A  multiple  sjot  team  autumn  Jus  a  sig¬ 
nificant  edict  in  rvtkirrj  the  use  of  available 
froguency  hands  arid  satellite  res'iiices  mon? 
efficiently.  if  w<'  consider  a  s,  stun  wlvv.e  an¬ 
tenna  is  arrang'd  to  fnnn  N  imi-» '.vrl,i|  pinq  :g « »t 
beiuns ,  each  serving  a  sep.irate  zone.  Fnch  bo.iiu 
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ci*  rises  an  u;  link  anti  a  d  link  at  a:  Merer. t 
cat  i  ier  1  tor...  os  ,  with  o  :.!  •  *ct  u  ,s  :vi  ;•*  L*  tvoen 
a;»t  ii  viil  ink:,  tin  'uuk,  a  nwitt.'h  ivitrix'in  the 

1. C'  llite,  «;!kwi  on/«,;tuiily  m  Fiu.  1.*  'Il'.e 
:.v;tc;.  r«it  rix  ll.t  e;  vi’iiin  •**!  r,  nvsVfrs  aixl  traus- 
"dtus  aco.'idnv;  to  a  [  u  -j  u  ■  cnbrl  algonUs-'.. 

F  u  si.rjlicity,  txrjba*.-  tialtic  :  yi  it  iy  is  assum¬ 
ed:  «  ach  z*  c,  stainav  th*’  s  l-:o  n\s:oi  ot  stations 

I.,  which  is  :  !k.*v;',  iy  ail;  i  =  l, . ,N  a:ul  ]*-l, 

. ,h.  Lack  statn-n  tn.rsj  its  |ocki-l  s  to  the 

n«:-o  issuer  uf  Mat  i  ns  Nh  and  all  ihU  :  st  at  u.n 
l«?  ket  ratiio  aie  oiual.  li.a  artival  p.x'ess  is 
Is;  iss»ih  ai«.i  all  nme-a.;es  of  a  suvgle  pic¬ 

ket..  All  buffer  cajmcities  aio  liifinite. 

In  [t'A-F.  Hf1]  ,  th.e  autlv  >i  s  oi  <*iie  different 
l  i  otocols  for  acv<rsini  a  sea  1 1  i 
1  sj  tx.ii'i'le,  in  the  Fixid  Amenrent  T.  *"w\;  th.e 
fi  *»•  litmctuie  is  si..  Wi  m  1  i<j.  7,  eack.  i>aith 
slat  ion  duio'ts  Jts  ti.dfic  to  th<*  desi  i  vi  des- 
tnuition  zone  by  1 1  an:cu  1 1  i  no  into  the  ;  reassigned 
si.t  during  which  the  switch  is  s*.>t  ter  tint  desti¬ 
nation.  'Jne  ir.sitchi::g  St^;a,  n*<*  m  ait  n-h^ud 
switch  n«it  i  ix  bdhvs  a  |  i  o.u  :  an;«  >1  sc!,  dale,  ;  <  5- 
ril  ly  ndifid  at  lnfM^rauit.  inteivals.  i  y  ci  ntn  1 
sj  :nals  frail  tl.e  at  mind.  Also  slot  Wd  Ali'iilA/'ilf-l 
!:y;hro  with  rn  1 1 1  j  j  1  *  *  uplinks  f«ive  h1*  :i  invest  ma¬ 
t'd  for  bn.odcas!  inj  satellites  in  |SlTVi  GO}, 

( VSvJl'.  79]  and  fl;id^)  70]. 

In  section  2.  we  d--  sci  ihoffie  i 'and.  in  'IMA,  as.,1 
th<  n  aitollyse  it  ;n  s-act  u  a  3.  In  srctlon  ■!  wv  d.is- 
cuss  t )  10  ir.pl . : »  nUiticn  of  the  sciari'*  n,  nulti-spot. 

h  Mil.  r.atel  1  ites . 

2.  'Ihe  Pcmd>  in  ’IMA  Frotcx^ol 

Ii*  a  o,  we  will  lirscmU'  the  l-h;uki  .  'iM-A  j  nio- 
ex>l  tot  a  gioUil  Umih  covot mg  the  entiut  area  t.? 
h‘  SeiVicixl. 

We  aie  c.’nsidt  r  m  :  M  id*  nt  ical  um'i  in,1  M 
us*’!!*,  aie  divnlixl  irur  N  ai.Mps.  I'.uh  »;i..uj  a 
s  i  ’  f  s  of  I,  us,,  is  { I  .e  .  ,  M  hi.)  .  I  a. ’ll  n:  >'i  h.  i  ;  an 
ltd  m  it..*  Ini!  I «  i  .  M.  a  j-  awivals  at  u.l:  u:  <  a 
f'd!^vr  t)je  slat  istic:;  of  a  l\>i*;son  pj  :  us';, 
late  Inch  i:»  a  ;•?  is  to  O'lv.ial  -  t  a 

s  mail1’ picket  .  rhann<‘l  tun  is  sK  tin!  nvl  all 
u  'is  aie  :>yi;chio::i/i>]  to  tfu:  ehmnel  t  uiu  T  l  ;<  ■ 
sh  d  b  iiijth  is  mpi.i1  tr>  th«*  ilmat  u  n  ol  n<’  ;ac- 
k<  t  .  C'hanr.e]  slats  ai*’  ijiims^xI  in  tram  s  eacli  rf 
w^iuh  consists  of  f»  (ni.srnit  ivs  slots. 

b»eh  gnnip  is  . i : • :% i « in.  x  1  a  tixt\i  ta\'>  ii.i  sh  l  in 
tii»*  (i.imy  as  m  tlie  rajn*  of  i  i  ;<\d  slots  ol 


7F.3.1 
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unliiuuy  *i UA,  .ed  u-  two  ,  invi*  t  1 ..I 

f  avoi  oil  1)  I  t  .  1‘ach  Us i  i iu  r  u  (  s  i  1 1  i  c »x i ) 


within  a  group;  ha:j  t, 


f  a vo i  ed 


u  us «*r  Uis  a  jock**'  ready  hr  Uans-* 
-lioaioil  LllO  ,1  jacket  Ur  M  VS  Ut  the 

IVud-of-the-'  If:#'  (ilOh)  j*u;i  ♦  i' n  jn  t ),•  ).uh<!  ,) 

**  *>■'■  "V  >1  -ilo-.  :M-d  :)<*..  »x  ... 

IU  th»*  hiV  ll*d  <  1  hiv,  v^itn  Jiti.S  .)A*y 

a  ai  id  j  n  any  <  t  i .« ■  r  jam  cn  lar  si ...  t.  wuJa  t  h. .  t 
N  slots  with  j  i ^  i.xiij 1 1 1 ty  b  "has; 

a  a  b  (h-l)  { ,  j 

Lo  t  us  c.vt  1  a  in  th  l  s  pi  cmc  s s  as  i r..  j  a  i  or  1 » T  •  e 
vKich  Kvs  t.Se  r.rivis  1  to  N.  Oo  at  t.V  o  r  j  tv  :s 
^>v.a; »  n;h  Aw>t  x%  th.-  .  ..v 

(>'»•*  i  r;n  I  I  ■  ,  ;;  ,.  .... 

!.  i  f  l  i  •  ''••  I:.  ‘  ■  t  it  ».  .  , . . 

llh  is  t  r  naoiit  Sd,  tin*  i.'uhttr  is  *  i.r  n- **  1  a:.,  i  I. 
th»G  Hub  {ticket  t.r  UnSiUl  S»S  li.U  IS  SC.1  i1  >  1  i  ]  • i  ]  ilOV 
to  the  result  o t  the  u»nli't  us  L"  in;. ;  idle  't  *'i i  »ds 
wtion  tin'll*  is  no  jackets  m  the  os  ;  *  a:; 

a  packet  arrives  at  the  butter,  the  roulette  scho~ 
dules  the  packet  tran.^mssion,  Y:  i.  3. 

If  a  ja.v.et  encvunt  oi  s  a  collision,  wluuh  i$ 
ackncvr  lo  i<  j*  %  1  a  1 1  e r  a  r  i  >v«  d - 1 1  ;  p  t  i  :* » ■  •  i  ■ .  *•  t  n 

]oms  tho  hut  for  at  the  Ilnd-of -the- Line  cm)  ^ra¬ 
tion  in  the  s. i;  «'  way  as  a  j^iskot  r>  *wly  nvr.>  a  «it « •*). 
The  qumioinj  :*. » h  1  with  tvutMfk  is  :.:,  vn  in  fn.k 
V>.c  nx’ssa-je  anival  rate  mcluimg  coliia*xi  jac¬ 
kets  is  ii“iv.u.4jii  by  >  . 

3.  Aju\ly_s  i  s 

U*t  1*}  ti.c  j^i.rkrt  do  lay  due  to 

tiiG  user's  beffoi.  ’i'i.e  o>;-oeted  joekot  delay  in 
the  systMTi  is  then  given  b> 

D  a  2P  ♦  Di  +  i:  (D;  +  21)  <2) 


t:  =  ti.e  averare  iu£il«»r  cf  retransmit ssio;is 
2P  =  round- trip  tie* lay. 

3 . 1  Ihe  Averag*. ‘  Pae kiyt  V\’a  1 1 1 T  i r * » 

'ro  fiul  the  value  of  D.  ,  wo  use  an  appaouch 
sunil  at  to  th*  one  used  try  Li*:i  1 1  AM  77  J  in  I,  is 
analysis  of  'II/IA.  We  ct»/i.snier  tw  knvls  of  sor- 
vie*1  tire  distrilut  ion  function  li(x)  and  U(x)  in 
a  sinqle- server  cju*ue  with  i'Oir.son  ai  rivals  at 
|k*ckets  (»*r  s<*conl  U(x)  is  the  distribution  of 
a  jocko t  which  initiates  a  busy  period  wil.li  first 
and  second  mu  mints  aid  b7.  Ml  subsequent  j'ac- 

kefs  ii  Lh*  sum  *  i  jc."  i ». .  : '  **!  hiv*‘  s*’;vi-'#  t  Hi* *• 


-1  ;;  u‘  x  pi^d  nm.ilar  to  I.ri'ii  \  i  *A,  U  is 
'Jivnby  1 


1  .U  -b.) 


1  2  t  1 


Mho  voe  ranijuiaticx'.  sijjwri  ui  A  we 

v<a  a  m 

*  '>-'•**  'y  *« 

a  *,4t;2^MN-2)^]  ;;  u~uhi 


5?  ishU-ir  ;  dep.  1 


N  *d 


*  \2id n  (h-2)b; 


-‘.d'i-pl 


+  u(y-T*N) 


^  ll-rWr'  N  “  yi]  0  :y.  T 


c  t  x  i  —  ,  ]  x  -  .) 


Kia:i  ops.  (n)  ,  ( 7)  ,  dd  wv  e*in  easily  Llaim  the 
t;i  st  aj.d  ui  nn.iits  of  jacket  service  Ua’fi 

br  h'  h-  h  • 

3,2  Tli1'  Av«*i\ss*  !VrU  r  of  pj*t  r an*. « .issi:jv:5  T 

ilv*  avcii.Ji’  r*i  i •  1 1 tin:-*  i ss;<  ns  torn 

<>1  j  *  i  |«a«.  ket  f,(  rm.i.;-'!  ;s  r<da*.oi  uvl'*r  iquili- 
hi  un  c\  n  lit  ions  t  >  hr  th«*  system  U\rciMli|’Ut  and 
G ,  tJie  cli« time l  trait jc; 

1  +  E  - 1  11 

wh'-re,  s  *'  T  aj»i3  G  -  '  .T.  lieie  a.,  tdie  pacd>- 
ahility  t  ii,i  t  a  j-acki-t  i?i  t  r.intsiii  t  toil  in  a  tavoie 
slot  and  t U’.o  probability  tlvtl  a  jack-.-t  is 


tiUt  f  if  >  ■>'  t  A  a  <  ■  t  ■  ,,  <  f 

*•  th/*  f.'S'J^r  of  k*'l  s  in  th"  :;y  <  t  \  [\t  ••  h  in 

picue  an*)  in  .e'lVie*)  at  t  inr*  t.  We  r rw  d'dinu 

►  f »  r ;i.‘  - f ; 


fid  I  /  u.d  P  t  or.-!  i..  j  r  }*.!■;  •  i  /  *  *  1  :  1  '  1 

collides  witli  •/ -Hi*  mi  l»  i  pi'h  ts,  a;  f  1  a. .sen.  i  in* 
a  ,,  .  i,  t  ,  .  -.ip  rm-N-r*:  d  i  **l  1  ansi’i  ‘-'a  >n 


B  (o)  «  /  e 


*  t  il  i  j )  -* 


P  =  ProbiA  successful  tranira ssjon  in  a  favored 


r  =  l  ul)  JodllS! ill  a  t-‘-  1 

+  Pieb  l  col  1  is  lens  ill  M- 1  ivei  1*1  i  :•  let :« \ 

M-l  M-l 

—  a  !  r.  *  in-ii  b  .  i\  (.4) 

.11 

i  i  i  i 

V  -  bob  li  |*u  v’l'i  in  lh-1)  1  r,'  i-'l  an  1  (•'  I  ) 

'l  1  ir  : i - 1  uwtoi  pickets  an.?  1 1  an/cu  1 t  ed  s  i ; *. :: * l - 
tanev  sis  1  y  ]  . 


_  y  ,  ]  . 

L-1C]  ‘*1  *  _‘V  1-1 

i  -  »uin  1 1 ,  1,-1 )  U 

P  -  piot-ii  pickets  in  I.  hcao-d  an.!  iM-I.-i) 

i.  n-  :  .iwi  >il  pickets  <uo  — 

Umtvuslyi 

-  V-C'V.-f,- V’  ' 
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Abstract 


A  reservations  scheme  for  multiple  access  suited  for  local 
networks  of  broadcast  channels  is  considered  and  analyzed.  The 
analysis  is  based  on  two  imbedded  Markov  chain  models,  one  for 
the  individual  user  and  one  for  the  entire  set  of  users.  The 
average  queue  size  as  well  as  the  average  delay  time  are 
obtained.  The  proposed  scheme  adapts  automatically  to  varying 
traffic  conditions. 
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The  problem  of  efficient  transmission  of  voice  and 
computer  data  through  shared  channels  remains  a  challenging 
and,  in  many  ways,  open  problem  in  the  communications  area. 

In  particular  the  problem  of  optimum  or  “good"  design  of 
multiple  access  protocols  for  radio  channels  continues  to 
receive  wide  attention.  There  is  a  special  need  for  protocols 
that  can  be  analyzed  so  that  their  performance  can  be  predicted 
and  evaluated  before  implementation. 


As  is  widely  known  by  now,  multiple  access  protocol,  s  may  be 
partitioned  into  three  main  categories.  The  first  is  Lhe  class 
of  fixed  allocation  protocols,  like  frequency  division  multiple 
access  (FDMA)  or  time  division  multiple  access  (TDMA).  These 
are  known  to  be  effective  only  for  heavy, non*bursty  traffic. 
Otherwise  they  are  known  to  be  extremely  wasteful  of  channel 
capacity  since  dedicated  allocation  of  resources  is  ill-suited 
to  low  duty -cycle  users.  The  second  class  of  protocols  consists 
of  contention -based  random  access  schemes,  such  as  the  ALOHA 
protocol  and  its  many  variants.  Their  main  disadvantages  are 
known  to  be  an  inherent  instability  and  a  low  channel 
utilization  rate.  The  third  broad  category  of  protocols 
includes  dynamic  assignment  schemes,  and,  in  particular, 
schemes  based  on  reservations.  Few  of  the  proposed  reservation 
schemes  have  been  analytically  studied.  Although  comprehensive 
reviews  of  multiple  access  protocols  abound  in  the  literature 
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(1-3]  we  briefly  recall  some  specific  protocols  selectively 
to  create  the  context  for  our  study.  In  [4]  a  protocol  is 
described  where  the  channel  is  divided  into  two  subchannels; 


one  operated  in  slotted  ALOHA  mode  for  reservation  requests, 
and  the  other  operated  in  a  dedicated  mode  for  data  packets. 
Reservation  schemes  such  as  this  one  are  particularly 
attractive  when  a  significant  part  of  channel  traffic  consists 
of  multi-packet  messages,  became  the  successful  transmission 
of  a  single  reservation  mini- packet  is  sufficient  to  reserve 
as  many  s 1 o t  s  as  are  necessary  for  the  entire  message. 

In  [5]  Crowther  et  a 1  .  introduced  the  Reservation -ALOHA 
idea.  In  this  scheme,  time  is  slotted  and  the  slots  are 
grouped  together  in  fixed  length  frames.  Slots  are  periodically 
assigned  to  users  as  in  T  D  M  A  ,  but  unlike  T  L)  M  A  t  h  e  s  e  are 
temporary  assignments.  When  any  particular  slot  is  idle  in 
one  frame,  the  slot  is  opened  up  for  contention  transmission 
in  the  following  frame.  A  user  who  successfully  transmits  in 
such  a  contention  slot  will  then  be  granted  ownership  of  the 
slot  as  long  as  he  has  packets  to  transmit.  This  sell  erne  was 
later  analyzed  by  Lam  [6], 

In  [ 7 ]  the  Round  Robin  reservation  scheme  was  introduced 
in  which  a  slotted  channel  with  fixed  length  frame  structure 
is  assumed.  Each  user  is  permanently  assigned  one  slot  p e  r 
frame,  and  unowned  slots  (which  are  present  when  the  number 
of  users  is  smaller  than  the  number  of  slots  per  frame)  art? 
made  available  as  needed  to  users  requesting  them.  Users  are 
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permitted  to  transmit  in  slots  owned  by  terminals  which  currently 
have  no  packets  to  transmit;  such  inactive  terminals  can  easily 
regain  ownership  of  their  slots  when  needed  (following  a  one 
f  r ame  delay)  . 

In  [8]  a  reservation  scheme  is  proposed  in  which  channel 
time  is  divided  again  into  frames,  each  of  which  consists  of 
reservation  time  slots,  preassigned  (PA)  time  slots  and 
reservation  access  (RA)  time  slots.  In  the  reservation  time 
slots  each  user  sends  a  reservation  message  (consisting  of  the  number  of 
packets  in  his  buffer  minus  one).  Each  user  is  assigned  one 
PA  slot  permanently  and  a  number  of  R A  slots  matching  the  amount 
requested  in  the  reservation  message. 

Wieselthier  and  Ephremidos  [2]  introduced  the  Interleaved 
Frame  Flush-Out  (1FF0)  protocols,  in  which  a  long  propagation 
delay  is  assumed  (satellite  channels).  Slots  are  reserved 
by  the  different  users,  but,  because  of  the  d e 1  a y  ,  there  may 
be  slots  that  will  go  unused  unless  they  are  opened  to 
contention.  The  performance  of  the  protocols  was  evaluated  by 
computation  and  simulation. 

In  this  paper,  we  consider  and  analyze  a  simple  reser¬ 
vation  scheme  essentially  identical  to  that  proposed  by  Rot  ha user 
and  Wild  [16]  and  studied  by  Bux  [17]  and  similar  to  the  ones 
proposed  by  S p a n i o 1  [18]  and  Mark  [19].  This  protocol  has  high 

adaptivity  properties  and  is  eminently  suited  to  low  propagation 
delay  environments,  namely,  local  network  applications.  The  usual 
assumptions  are  made  about  user  statistics,  channel  properties, 
etc.  as  will  be  described  in  the  sequel.  Our  analysis  is  new  and 
is  based  on  a  modeling  idea  introduced  in  [9],  that  consists  of 

two  coupled  Markov  chains,  one  that  describes  the  status  of  the  bullet*  contents 
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of  a  typical  user  (the  User  M . «  r kov  Chain  )  and  one  that  describes 
the  status  of  all  the  users  of  the  channel  (the  System  Markov 
Chain) .  Such  an  approach  p remits  us  to  account  for  the  interaction 
between  the  users,  quite  similarly  to  the  case  studied  in  [9J. 

In  section  II,  we  describe  the  channel  frame  structure  and 
the  user  model.  In  section  III  we  describe  the  system  Markov 
chain,  and  in  section  IV  the  user  Markov  chain.  Delay  is 
calculated  in  s e c  t  i o n  V  .  Finally,  numerical  solutions  and 
simulation  results  are  presented  in  section  VI.  Two  Appendices 
provide  the  analysis  details. 

II .  The  Model 

In  this  section,  a  description  of  the  reservation  scheme 
is  presented.  This  is  done  by  first  describing  the  channel 
time  frame  structure  and  then  i  he  user  model. 

a )  T  h  e  F  r  a  m  e  structure 

The  frame  structure  is  shown  in  Fig.  (l.a),  where  the 
first  slot  in  the  frame  is  the  reservation  slot  which  is 
subdivided  into  M  mini-slots,  ( M  is  the  number  of  users  in  the 
system).  Each  user  is  permanently  assigned  one  mini-slot  in 
the  reservation  slot,  in  which  he  sends  a  request  whenever 
he  has  a  packet  ready  to  transmit.  In  particular  user  i  is 
assigned  the  ith  mini-slot  in  every  reservation  slot.  To 
equalize  delays  among  the  different  users  this  assignment  may 
be  cyclically  rotated  so  that  no  individual  user  has  a 
permanent  advantage  over  the  others. 


Mini-slots  are  very  s  h  o  r  t  since  they  need  only  a  c  a  u  mm  od  a  i  e 
the  transmission  of  a  single  tune  indicating  the  need  to 
reserve  a  slot  in  the  next  franc.  The  length  of  the  main  frame 
is  a  random  variable  and  depends  on  the  number  of  users 
requesting  transmission  at  tin*  beginning  of  the  frame.  For 
example,  as  depicted  in  Fig.  i ,  users  number  2  and  number  4  are 
requesting  transmission.  Since  user  1  has  nothing  to  transmit, 
user  2  knows  he  will  be  the  first  to  transmit  after  the 
reservation  slot  is  over.  Also,  since  user  3  did  not  request 
transmission,  user  4  knows  he  will  transmit  after  user  2  ,  and 
then  the  new  frame  will  start.  Thus,  the  frame  length  equals 
to  three  slots.  Note  that  the  frame  length  cannot  be  less  than 
one  slot  and  cannot  be  m ore  than  (M+l)  slots  (corresponding  to  the 
case  that  nil  M  users  are  requesting  transmiss  Lo  n )  .  W  e  a  s  s  u  m  e  that 
messages  consist  of  single  packets  each  of  which  fits  exactly  in 
one  slot.  Also  we  assume  that  t  h  e  channel  is  m  o  n i t  o r  e  d  by  all  users 
without  delay  or  noise.  Note  also,  that  a  user  transmits  only  one 

Dacket  ner  frame  regardless  of  the  number  of  packets  in  his  buffer. 

b )  The  User  Model 

The  system  consists  of  M  terminals,  e a c  1 1  of  which  has  an 
infinite  size  buffer.  The  arrival  process  at  each  ol  the  M 
terminals  is  assumed  to  be  a  li e r n o u  1  i  i  process  with  rate  c  . 

The  user  may  be  in  one  of  two  states,  an  idle  state,  if  his 
buffer  is  empty  or  an  active  state,  if  h  i  s  buffer  is  non-empty. 

As  shown  in  Fib.  l.b,  whenever  the  user  has  a  packet  in 
his  buffer,  he  sends  a  reservation  request  at  the  beginning  of 


the  frame,  and  at  the  same  time  the  packet  at  the  head  of 


line  (HOL)  is  moved 

t  o 

the  transmitter  and  wa 

i  t 

s  t  h 

e  r  e  until 

it  is  transmitted  in 

i  t 

s  a ss  if'iied  slot  . 

We  need  the  following  du  t  i  n  i  t  i  oils  : 

Let 

TJ 

A 

steady  state  j>  r  o  b  a  b  i  1  i 

t  y 

o  1 

h  a  v  i  n  g  i 

i 

= 

packets  in  the  buffer 

a  t 

t  lie 

beginning 

of  a  fra  m  e . 

Hence 

no 

= 

p  r  o  b  a  b  i  1  it  y  o  1  a  b  u  f  f  e 

r 

h  e  i  n 

p>  t  y 

= 

probability  of  a  us  e r 

be 

i  n  g 

idle. 

1  ~,:o 

= 

probability  of  a  user 

l)  e 

i  ng 

active  . 

Such  a  model  produces  1 1 

non-standard  q  u o u o  i  n  g 

syst  e 

m  o  f  in  t e  r 

acting  queues .  The 

service  time  of  a  p a c  k  e  t 

i  n 

o  a  e 

q  u  e  u  e 

depends  on  the  status  of  the  other  queues.  If  we  let  n  represent 

K 

the  contents  of  the  input  buffer  of  the  user,  then  a 

complete  description  of  all  M  users  requires  the  specification 
of  the  joint  probability  distribution  P  ( n  ,  n  ,  .  .  .  ,  n  )  .  The 
determination  of  these  probabilities  is  very  complex,  if  not 
impossible,  and  demands  the  solution  of  a  1  a  r  g  o  n  u  in  her  of  sets 
of  equations.  Fayolle  and  Iasnogorodski  [10]  showed  the 
limitation  of  this  direct  approach.  They  considered  the  simple 
example  of  two  exponential  s  or v  e  r s  t  h  e  s  e  r  v  i  c  e  rate  of  e  a  eh  of 
which  depended  on  the  status  of  the  other  server's  queue. 

Using  the  theory  of  analytic  continuation  they  reduced  the 
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determination  of  the  joint  probability  distribution  of  the  two 
queues  to  a  Riemann-Hi 1  be r t  p  rob  ’em  and  were  able  to  obtain  a 
very  complex  closed  form  for  the  solution.  Furthermore  extension 
of  their  analysis  to  the  case  of  more  than  two  users  is  not  feasible. 
Leibowitz  [11]  presented  an  approximate  method  of  treating  muitiqueue 
systems.  He  studied  the  case  of  N  queues  with  a  single  server,  in 
which  the  queues  are  served  in  cyclic  order.  He  made  the  assumption 
that  the  distribution  of  the  queue  size  of  one  queue  is  invariant  if 
it  is  the  same  for  all  queues  during  one  cycle. 

11  a  s h  i d a  [Id,  13]  used  L e  i  b  o w  i  t  z  1  s  approximation  in  the  analysis 
of  muitiqueue  systems.  Also  in  |  S ]  ,  the  authors  used  Lei  bow  i  t  z  '  s 

approximation  a  Ion g  with  an  i  n d  e p end e  n  c  e  as  s  u m p t  i  o n . 

We  consider  that  the  main  contribui ion  of  our  paper  is  the  proposed  mathematical 
model  for  t  h e  analysis  of  interacting  buffered  terminals.  The  model, 
as  we  mentioned  earlier,  is  basically  a  combination  of  two  Markov 
chains  imbedded  at  the  beginning  of  each  frame;  one  Markov  chain 
that  describes  the  state  of  the  user,  r  e  f  e  r  e  d  to  as  the  user  Markov 
chain  jand  another  Markov  chain  that  describes  the  state  of  all  the 
users  in  the  system,  ref ered  to  as  the  system  Markov  chain. 
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111  The  System  Markov  Chain 

The  state  of  the  system  is  .escribed  by  the  number  o i  users 
requesting  transmission  at  the  be  g  i  nni  n  g  of  the  f  r  a  m  e  .  first,  we 
need  the  following  definitions:  bet 

j  -  n  u  u  her  of  users  rei’uost  i  n  g  t  runs  m  ission  at  the 
bei;iun  i  n g  o  i  a  f  r  a r. <. ■  . 

t.  «=  t  rame  length  when  j  users  request  t  v  a  u  s  i;i  i  s  s  i  on  at  the 
J 

b  e  g  i  n n  i  n  g  of  t  b e  f r  a  m e 
=  1  +  J 

p  j  A  s  t  e  a  d  y  state  p  r  o  b  a  b  i  1  i  t  y  o  f  j  u  s  e  r  s  r  e  quo  s  t  i  n  g. 

transmission  at  the  be.*,  inn  iiq;  of  the  l  r a e 

o  .  A  V  r  o  b  a  b  i  1  it  v  a  n  idle  u  s  e  r  v:  e  n  e  r  a  L  e  a  t  1  e  a  s  t  o  n  e 
J  - 

packet  during  a  f  r  a  m  e  ul  leu  g  t  b  t  1  t  j  )  slot  s 
=  1  -  (l  -  o)i+1 

F  =  p  r  v)ba  b  i  I  i  L  y  o  i  one  ;m  c.  be  t  on  l  y  in  the  butler  a  t 
the  beginning  of  a  t  r  a  m  e  g  i  v  e  n  user  is  a  c  t  i  v  e 


The  system  M arkov  chain  is  ill  u  s  t  r  a  t  e  d  in  Fig.  2 . a .  w h i 1 e 
its  transition  probabilities  are  obtained  in  a p p e mi  i x  A . 

IV  T  h  e  1)  s  e  r  M  a  r  k  o  v  C  h  a  i  n 

The  state  of  the  user  is  described  by  the  number  ol  packets 
in  L  h  c  buffer  at  the  b  e  g  i  n  n  i  n  v;  o  t  t  h  e  f  r  a  m  e  .  hot 

n  =  number  of  packets  in  b  u  f  f  e  r  at  t  h  e  b  e  g  i  n  n  i  n  g  o  L 
a  f  rame 


tt  =  steady  state  probability  of  n  packets  in  the  buffer 
n 

at  the  b  e  g  i  n  n  1  n  g  of  a  f  r  a  m  c 


Let  ;r  ( z )  be  the  generating  function  of 


1  .  e  , 


TT  (  Z  )  =  E  71  ZU 

n  =  0  U 

In  Fig,  2  .b  ,  we  show  the  user  Markov  chain.  Note  that  only 
one  packet  can  leave  the  buffer  per  frame.  We  show  the  derivation 
of  the  transition  probabilities  in  Appendix  b. 
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V  Delay  Analysis 

The  average  total  delay  per  packet,  D,  Is  the  same  of  the 
average  waiting  time  in  the  buffer,  W,  plus  the  average  time  from 
the  beginning  of  the  frame  during  which  the  packet  is  scheduled  for 
transmission  until  the  packet  is  transmitted  In  its  assigned  slot. 

We  refer  to  the  latter  as  the  service  time,  S  . 

Hence , 

D  =  W  +  S 

Using  Little's  result,  we  can  write  W  as; 

h  -  a 

o 

where,  Q  is  the  average  queue  size. 

Assuming  the  cyclic  assignment  discussed  earlier  in  section  II, 
the  average  service  time  equals  to  one  plus  half  the  average  number 
of  requesting  users  at  the  beginning  of  the  frame,  i.e, 

s  =  1  +  j  J 

V I  Numerical  Solutions 

Equations  (A-l)  through  (A- 5)  provide  the  values  of  the  quantities 
that  were  needed  in  order  to  solve  the  state  transition  equations 
for  the  system  probabilities  p ^ .  However,  some  of  these  quantities 
are  not  solely  expressed  in  terms  of  the  parameters  M  and  o,  but 
also  in  terms  of  F,  i.e.,  the  probability  of  one  packet  in  the  buffer 
at  the  beginning  of  a  frame  given  user  is  active,  which  is  a  function 
of  it  ^  and  tt  ^  .  But,  tt  ^  and  tt  ^  both  depend  in  turn  on  {  p  ^  }  .  Thus  ,  in 
total,  we  have  a  set  of  simultaneous  coupled  non-linear  equations  in 
tt  ,  tt  and  {p  },  In  Fig.  3,  we  show  the  flow  chart  for  the  numerical 
procedure  used  to  solve  these  equations.  We  start  with  an  initial 
value  for  F,  then  solve  (M+l)  linear  simultaneous  equations  in  { p  ) . 

We  then  determine  and  and  use  Iiq  .  (B.l)  and  Wegestein’s  iteration 


J  7 


i 


scheme  [14]  to  obtain  the  new  value  of  F.  Two  to  five  iterations 

were  required  for  convergence  to  a  solution  for  F  within  a  10  ^ 

tolerance  in  the  values  of  F  and  In  }. 

n 

For  M  =  4  ,  Fig.  4  shows  the  probabilities  7i  ^  (the  prob  .  of  an 

empty  buffer),  p  ^  (the  prob.  of  zero  users  requesting  transmission) 

and  p  ^  (the  prob.  of  M  users  requesting  transmission)  versus  channel 

throughput.  Notice  that  as  throughput  increases  both  ^  and 

decrease  while  increases.  Figs.  5  and  6  show  the  average  waiting, 

service,  and  total  times  versus  throughput  for  M=4  and  8  respectively 

Throughput  is  defined  as  A  =  Mo  .  Such  a  definition  presupposes  steady- 

state  stability  so  that  the  average  input  rate  is  equal  to  the  averag 

transmission  rate.  Although  detailed  stability  analysis  is  not 

considered  here,  it  is  fairly  obvious  that  since  there  are  no  wasted 

slots  (due  to  contention  or  idleness)  except  for  one  reservation  slot 

in  each  frame  and  since  all  users  are  assumed  similar,  any  value  of 
M 

A  less  than  produces  a  stable  equilibrium. 
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V  I  Cone  Iasi  on 

In  tli  is  paper,  we  have*  intr  need  a  :i  d  analyzed  a  res 
scheme  ,  in  which  reservation  pac  ts  are  bro accost  in  TbM 
over  a  traction  of  the  nain  channel.  The  average  wait  in.;, 
and  total  times  h  a  v  e  bee  n  o  b  tai;u-  I  .  ike  n  a  1  y  s  i  s  in  b  a  s  e 
combination  of  two  imbedded  Mark  v  chains,  in at  partially 


for  the  inter  p  e  n  d  o  u  c  e  between  t : .  e  u  s  era.  b  a o  h  n  a  p  p  r  e  a  c 
i  tv.  p  1  e  ::t  e  n  t  e  d  by  the  a  u  t  h  o  r  s  in  t  h  anal  y  s  i  s  of  si  o  t  t  e  d  -  A 1,  C 

finite  n  u  :u  b  e  r  of  b  altered  a  s  e  r  s  [  V  ]  .  A  1  s  e  in  'lb;,  we  y  e 

the  a  n  a  1  y  s  i  s  o  l  the  reservation  n e  h  e e  p  r  e  s  e  n  ted  in  i  h  i  a 
a  c  Ivins  o  1  Col i  is  i,  n  R  e  s  o  1  u  t  i  o  n  A  1  o  r  i  L  hr.  s  o  i  the  C  a  p  e  t n 

a  s  s  a i  n  >’  a  f  i  n  i  t  e  u  ;.tb  e  r  of  b  u  f  f  e  r  e  d  u  s  e  r  n  . 


(j 


A  P  P  l  N I  ■  l  X  _A 

System  Markov  Chain  for  the  Reservation  S  c  1 1  e  m  u 

Let 

K  =  number  of  active  users  out  of  the  j  requesting  users 
The  distinction  between  the  number  of  active  users  and  the  number  of 
users  requesting  transmission  is  necessitated  by  the  assumption  that 
if  there  is  only  one  packet  in  the  buffer,  it  is  automatically  moved 
to  the  transmitter  box  at  the  beginning  of  the  frame  and  hence  the 
corresponding  user  is  considered  idle  at  the  beginning  of  that 
frame.  Thus 

K-  K.  =  number  of  idle  users  at  the  beginning  of  a  frame 
Recall  that  F  is  defined  as; 

F  =  Probability  of  one  packet  only  in  buffer  at  the 
beginning  of  a  frame  given  user  is  active 


.»  1 


Now  we  can  write  h ( i , K | j ) ,  the 


joint  probability  distribution 


of  i  and  K  conditioned  with  respMl  i o  j,  as  follows: 

h(i,K|j)APr[i^out  of  M-K  id  I  <.■  asi;rs  generate  packets  during  a 
frame  and  K  out  of  the  j  requesting  users  are 


active  given  j  users  are  requesting] 


(•'  K)o![  1-0  .  ]M~K_i(l!)  (1-F)K(F)  j-K 
1  J  J  K 


*  q(i»K,j)C^(j) 

Referring  to  Fig.  (2. a)  we  can  write  A^(n),  the  transition 


(A.  4) 


probability  front  state  j  to  state  n,  as; 
in  i  n  ( j  ,  n  ) 

A  .  ( n  )  =  I'  P  r  (  ( n  -  K  )  arrivals  and  K  users  active  given  j 

3  K=0 

users  are  requesting] 


m  i  n  (  j  ,  n  ) 

h  h  (n-K  ,  K  ,  j  ) ,  1  <j  <M  ,  0<_n<_M 

K  =  0 

(A.  5) 

z  M v  n  . .  x  M-n 

(n)o  (l-o)  ,  j  =  0,  0  <^n  <_M 

Of  course  then  we  have 
M 

p  =  I  A  .  ( n  )  p  .  ,  0_<  n  <  M 

n  j=0  J  J 
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A  V  P  KMJl  .*  ji 

I  h  e  User  Mar  k u  v  d ha  i  n  fur  r  he  K  e  s  e  r  v  a  i  i  o  n  Scheme 
In  this  Appendix,  we  determine  the  transition  p r o b a b i 1  i  t i e s 
lor  the  user  Markov  chain  shown  in  F  i  g .  2 . b .  Also,  expressions 

tor  :i  ^  ,  the  probability  ot  an  e n. \  t  y  b  utter  and  ,  the  butter  size 
are  obtained.  First,  we  need  the  following  definitions;  Let 


b  ,  1  ( n )  =  Y.  P  r  o b  .  [one  packet  leaves  and  no  packets  arrive 

n  + 1 

J 

given  j  requesting  users  at  the  beginning 

o  f  a  fra  m  e  ] 

M 

=  l  '•?-  g(0,j+l)p.]/(l-p0)  n  =  0 ,  l  ,  .  .  .  (B.3) 

j  =  l  J 

Notice  that  the  initial  value  of  the  summation  is  j  =  I,  since  at  the 
beginning  of  the  frame  we  have  at  least  one  user  requesting  transmission. 
For  the  same  reason,  we  normalize  by  dividing  by  ( 1 “ P  q  )  • 

Similarly, 

f ' n  C n )  =  ^  ^r*  [one  packet  leaves,  and  one  arrives  during  a 

frame  given  j  users  requesting  at  the  beginning 
of  the  frame] 


?  1 


M 

::  g(i,j  +  i)p.]/(i-pn) 

j«i  J 

M-l 

E  g(0, j  +  l)p  ] /  (  1-p  ) 

j  =0  J 


n  r'O 


n  -  0 


(B.4) 


Note  that  for  n*0,  we  can  have  a  ;r.  a  x  i  m  u  m  of  ( M  - 1  )  requesting  users; 
hence,  we  need  to  normalize  by  dividing  by  the  probability  of  having 
at  most  (M-l)  requesting  users. 

Similarly  , 

B  ^  ^  (  n  )  =  £  Pr.  [one  packet  leaves,  i  + 1  arrive  given  j 

j 

requesting  users] 


M 

=  [ B(i+1  ,j  +  l)p.l/(l-Py) 

j  =  i  J 


n  =  2,3,... 

1  ,  m  i  n  ( n  -  1  ,  M  ) 


and 

M-l 

BQ(n)  =  [  g  (n  ,  j  +1 )  p  )  /  U  -p  )  I)  =  ],...,  M 

j=n-l  J 


Hence , 


n  -  I 

B  .-.(h)77  ,-,+B  ( n  )  7T  +1  B  .  ( n  )  tt  ,+Bn(n)r  2<n<M 

n  +  1  n  +  1  n  n  .  .  n-i  n-i  0  0  —  — 

i=  1 


M 

lB  ( n )  77  -Hi  (n)7i  +  £  B  ,  (n)<r  .  ,  n>M 

n+1  n  +  i  n  n  i - 1  n-in-i 


(B.5) 


(B.6) 


( B  .  7  a  ) 


ttq  =  B1(0)7t1  +  B0(0)tt0  (B.7b) 

t  ^  (  1 )  1T2  +  Iil  ( 1 )  7T1  +  B0  ( 1  >  *0  (13.7c) 

To  obtain  Tr(z),  we  need  to  define  the  following,;  Let 
M 

G ( 0 ) A  £  g(0,j+l)p  ( B . 8 ) 

.1=0  J 

M 

G (1 ) AE  g(l,j+l)p. 

~  1=0  J 


T  / 


(  B  .  9  ) 
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Multiplying  E  q  .  (B.7)  by  z  and  summing  over  n,  we  get; 


tt  (  z  )  =  E  tt  z 

n  n 
n  =  0 


=  1  B  *  ( n )  it  z  + 1  U  ( n  )  tt  z  1;. 
n  n+1  n  +  1  n  n  n 

n  =  0  n  =  0  n  =  0 l = 1 


min  (n-1 , M ) 


B  .  ( n  )  t:  .  z  u  (  n  - 

n  ~  i  n  -  i 


+  ^  B  0  (  n  )  11  0  Vm  u  ^ n  ~ 1  ^ 

n  =  0 


(B. 10) 


Using  Eq.  (B.3),  we  get 


1  B  .  ( n )  tt  z 

n  n+1  n+1 

n  =  0 


yi— [G  (0  ) -g  (0 , 1  )  p0  1  z“  1  { Tt  (  z  ) -t:  0  ] 


(  B  .  1 1  ) 


Using  Eq .  (B.4),  we  get; 


n  1 


Z  B  ( n )  7r  z  - 
A  n  n 
n  =  0 


From  (B.6)  we  get; 


y—  )(G(l)-g(l,l)p0n^(2)-:r0j  + 


[  G  (  0  )  -  g  (  0  ,  M  + 1  )  p  ^  ] 


(  B  .  1  2  ) 


M  j 


£  B  ^  ( n  )  ii  0  z  n  =  K  }'.  g(n,j  +  l)2np.un/(l-pM) 


n  =  1 


j=l  n=l 


j  0 '  '  M  ’ 


0 _  [H(z)-G(0)-H(z,M+l)pM+g(0,M+l)pM] 


(1_PM) 


(B.  13) 


Finally,  we  have 


00  min(n-l,M) 

Z  Z 

n-2  1  =  1 


M  w 

B  .  (  n  )  tt  ,  z n  =  1  Y.  B  .(11)77  .  z  n 

n-i  n-I  ,  ,  . , i  n - l  n-i 

1  =  1  n  = i + 1 


M  M 


*  T~ -  E  £  g(i  +  l,j+l)p.zi[7i(z)~7r  l 

A'p0  i=l  j=i  J  u 


i-p 


M  j  j. 

-  y  I  g(i  +  l,j+l)p.z  l u ( z )  —  n 
0  j=l  i=l  J 


0 


1  - 1 

- — —  X  z  lH(z,j  +  l)p  -g(],j  +  l)p.z-g(Otj  +  l)p  )(:i(z)-7i0] 
0  j=l  J  J  J 


~(1  _  p  ~ )'  [  (  z  )  -n  q  ]  (H(z)-zG(D-G(O)  ] 


(B. 14) 


From  (fi.ll)  -  (15.14)  and  (15.10)  wo  obtain 


,  \  X(z) 

1,(50  =  yTzT^o 


(15.  15) 


whe  re , 


1“P0  ' 

XU)  =  t 7 Tn  -  -  z-l]H(z)+p  (g(0,l)+zg(l,l)>-  -7-n: - rll(z,M  +  l)p 


IB. 16) 


Y  ( z )  =  (l-p0)z  -  ll(z)  +  p0(g(0,l)  +  zs(l.l)) 


(15.17) 


To  obtain  n^,  we  take  the  limit  in  E  q  .  (15.  15)  as  z  •  1  ,  hence, 


tTq  =  Y  (  1  ) 


(B. 18) 


whe  re 


Ul)  ,  -  i  -  p0d-«)  -  inn 


(B. 19) 


X^>  ■  Q-; p")  Up,,-  ^)H(l)  +  l-p0  +  op0(l-pH)-(l-p0)p^(M+ 1)0  +  1)  ) 


(B.20) 


To  obtain  Q,  we  differentiate  (B.15)  and  take  the  limit  as  z  -►  1  ; 


Thus 


X  (  1  )  rr  Q  -  Y  (  1  ) 
2Y  (1) 


(B.21) 


whe  re 


3  z  M 


[ H ( 1 , M+l ) t-2  H  (  l , M+ 1 ) ] J 


(B. 22) 


Y(l)  =  H ( 1 ) 


(B. 23) 
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AN  ALGORITHM  FOR  SCENE  MATCHING 


By 


Tarek  N.  Saadawi  and  George  Eichmann 


Abstract 

A  fast  algorithm  tor  scene  matching  is  studied  and  compared 
with  other  algorithms .  The  measure  of  similarity  is  the  sum  of  the 
absolute  values  of  the  differences  of  image  intensity  between  the 
template  and  the  scene.  Conparison  with  other  measures  of  s  ini  lari 
is  presented. 
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(I)  Introduction 


The  need  for  analyzing  the  images  of  the  earth  taken  by  various 
sensors  on  high  and  low  altitude  aircraft  is  increasing  at  a  rapid  pace. 

These  images  have  been  used  for  several  purposes f  including  navigational 
guidance,  geographic  and  topograpni c  map  matching,  natural  resource  analysis, 
weather  prediction,  and  environment  studies. 


A  common  problem  that  arises  in  these  applications  is  scene  matching. 
Given  a  pictorial  description  of  a  region  of  scene,  determine  which  region 
in  another  image  is  similar.  Tho  method  used  to  solve  this  problem  in  its 
simplest  form  is  called  template  matching.  An  image  of  template  is  given, 
an^  it  is  desired  to  determine  all  locations  of  the  template  in  another 
image.  If  a  sufficient  number  of  corresponding  regions  could  be  determined 
then  the  image  could  be  accurately  registered,  facilitating  location  of 


corresponding  objects. 

Including  the  basic  template  matching  method,  several  approaches  and 
technique 7  have  been  used,  such  as  sequential  similarity  comparisons, 
matching  with  invariant  moments,  and  matching  with  edge  features .  An  imaae 
edge,  as  defined  here,  is  a  local  chance  or  discontinuity  in  imaae  luminance. 
Generally,  imago  matching  can  be  shown  in  the  fol Jawing  functional  diagram. 


Sensor. 


Image 

Input 


Preprocos 


Edoe  Match i no 

Extraction. - Algorithm 


Reference 

Image 


Preprocessing  is  an  important  procedure  in  imaae  matching.  Its  basic 
objective  is  correcting  two  degraded  images  (such  as  two  images  from  two 
difference  sensors)  to  get  new  images  to  meet  cei.ca.in  requirements  of  accuracy 
and  projection.  Generally,  it  includes  geometrical  correction  and  intensify 
correction.  Geometrical  correction  is  a  transform  to  change  each  pixel  of 
distorted  imago  from  original  coordinates  (:<,  \0  to  metrical  ly  more 
corrected  coordinates  (u,  v) .  For  intensity  correction,  a  rat;  chin:  v  runs  ion.- 


BC/) 

at  ion  nt1  the \1  cal  led  Karhonon-I  jocwc  transform  was  widely  used. 

Edue  information  is  important  for  both  human  and  raehine  perception 
and  a  considerable  amount  of  research  has  been  done  in  attempts  to  extract 
edge  information  from  pictures. 

In  this  paper  wo  study  an  algorithm  for  image  matching,  Tno  algorithm 
requi res  less  computation  time  and  achieves  a  high  registration  performance. 

In  the  next  section  we  discuss  the  algorithm  and  in  Section  III  we  analyze 
the  results.  Section  IV  compares  the  performances  of  the  different  measures 
of  similarity  for  image  registration.  The  conclusion  is  given  in  Section  V. 
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i  1  An  A !  v  :o  1 '  i  tin  u  ft  i  r  Scene  *  Match  i  1 1>  • 

n- 1 1  1  i  wm  iuitl  M u~t  in ,  in  (1),  nro:  * \sed  a  rou:*sc  si  urdi  corr «. -1  at  if >n 
t< ‘dm iquL*  for  image  i\x.ii  strut  ion,  tine  siAiuential  similarity  detection 
algorithm  (SSDA)  .  They  used  the  normalized  correlation  coefficient  as 
tiie  measure?  of  match.  The  normalized  correlation  coefficient  is  dr-fined 


C(i,j)  = 


M-'e..-  4  e 


l  (M  a,  -  cj  (M* 


o  •'  i ,  j^L  -  M  -I  1 


where 


K  M 


q=i 


\Sij (p,q) f 


M  M  ii 

0,  =  C  c  S  J (p,q) *T(p,q) 

p=-  q=i 


Where  T(p,q)  denotes  the  template  with  dimension  MxM  pixels,  S1^  (p,q)  is 
MxM  sub  image  of  the  scene  located  at  (i,j)  coordinates  with  dimension 
LxL.  Hence  the  search  area  would  be  (L  -  M  +  1)  x  (L  -  M  +  1)  .  The 
formula  in  equation  (1)  requires  a  large  computation  time. 

A  critical  issue  is  hew  to  choose  a  thresholding  technique  for  terminating 
the  calculation  of  C(i,j)  early.  If  the  thresholds  are  too  high,  we  may  elimi¬ 
nate  the  best  registration  because  of  the  early  behavior  of  the  calculation  of 
C(i,j).  On  the  other  hand,  if  the  thresholds  are  too  low,  many  poor  registra¬ 
tions  will  avoid  early  termination  of  eqn.  (1).  The  authors  in  (1)  proposed 
an  empirically  derived  function  as  the  threshold  sequence  dependent  on  the 
subset  number  b. 


li-uv,  Wt 1  i  :i\\  st  i .  re .  ■  ur.  in:  t  i  u  *  i/.im  o!  th<  a}  aw.]  ut 


ala?it\  :-i  -usure::i  *ni  i  1  # 


i  s  c  k  ■ :  i  ru  \ :  as 


valns  c  :  th 


!  ^  \  mc  '•  ;  x. -t  v.\  '•!  i  two  ironies  .  i : '  t  ’  k  •  ::  *  \  isu  ■  a:  si  mi  ]  :m  it.'/.  Tji 


A(i,  1) 


;  Sl  j(p,q)  -  T  (p,q) 


The  calculation  procedures  arc  similar  with  the  SSDA  ah.:or.i  thm.  First, 
wo  use  a  subset  of  tlio  two  data  fields  (S1^  and  T)  at  a  1  ow  resolution  tc ■ 
calculate  an  estimate  of  A^i,  j)  .  For  example,  7q(i,j^  would  bo  calculated 
using  one-sixteenth  of  the  data.  This  A  (i,  j)  could  then  be*  compared  to  a 
threshold.  If  it  passes  the  threshold,  the  next  estimate  A,  { i ,  1 )  at  the 
next  high  level  resolution  is  calculated  using  an  additional  one-sixteenth 
of  the  data.  This  procedure  is  continued  until  all  of  the  data  points  are 
used,  giving  the  best  estimate  A^(i,j)  or  until  a  threshold  is  not  exceeded, 
which  aborts  the  evaluation  process. 

The  sampling  rule  of  the  data  set  is  shewn  in  Fig.  1.  The  coarse 
spacing  is  taken  in  a  pseudo-random  manner.  Each  pass  requires  that  one 
additional  pixel  be  chosen  from  each  block.  The  choice  of  the  sampling 
method  is  important  and  affects  the  matching  procedure.  This  will  Lx: 
s t ud i ed  i n  the  f o 1 1 owing  section . 

The  threshold  sequence  is  taken  as  a  constant  value.  It  is  shewn  that 
the  threshold  constant  is  between  301  to  501  of  the  average  value  of  the 
intensity  per  pixel  in  the  template  image. 
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intensity  distribution  and  background  no  is*.  . 

Tho  U mpl ate  consists  of  character  'T"  an.:  its  dimension  is 
its  intensity  distribution  has  boon  quantized  into-  12S  lewis.  Within  "I"" , 
the  intensity  has  Guassion  distribution  plus  Gaussian  noise-  with  : xto  near: 
and  ~  4 3 ,  tiio  background  (outside  ”t-,M )  has  randan;  noise  with  intensity 
quantized  levels  between  0  and  30. 

Within  "1  "  the  intensity  distribution  can  be  expressed  as 

Vx'y)  =  128  xky  \t!X;)["(x~y-k'r  ''  “T  ]  .e>X)[-(y-y._k)  7ak  ] 

•  rect  j  (x-xp/.p.  j  .  rect  f  (y-yk)/.‘yk  j  (3a) 

'••hero  "F"  is  broken  into  K  blocks  (K  =  3) 

rect  (X/.ix)  function  is  defined  as  the  "slit"  function; 

1  x  •_  ;.x/2 

roct(X  7.x)  -  .  . 

0  |>:  -  ^ 

and  tiie  Gaussian  noise  has  tho  form 

i,v,  (?.  )  =  KXP  [-  iq  -  \43x43)  ]  (3b) 

U.Lsido  "P,"  the  background  noise'  distribution  is 

'nil  =  lNT{  150  x  (l  +  2  x  |RAN(i)  -  0.5 

New,  we  express  the  template  as 


t 


Xnb 


1  +  0.5 


(4) 


The  template  and  scene  images  are  shown  in  fi cures  2a  me:  2b.  The 
total  and  average  intensity  of  both  the  template  and  scene  arc  giver,  in 
Table  1.  The  ratio  of  average  noise  intensity  to  average  signal  intensity 
is  very  high. 

B)  Similarity  Measure 

In  image  registeration,  the  important  criteria  is  the  decree  of 
similarity  between  the  two  images.  The  measure  c:  similarity  should;  be 
able  to  provide  a  lurce  separation  between  the  match  point  and  all  other 
test  locations.  Different  values  for  the  measure  of  similarity  are  shoo:, 
in  Fig.  3.  We  notice  that  the  matching  points  (denoted  by  +  )  have  lower 
distinct  ( i ,  j )  values  from  those  of  tiie  non -matched  points. 

C)  The  Sampling  Rule 

The  rule  to  pick  up  the  samples  can  affect  tip  mat  china  procedures. 

In  Fig.  4,  four  different  coarse  spacings ,  randomly  chosen ,  have  ixvn 
considered.  Figures  f.a,  5.b,  5.c  show  the  dependence  of  the  r  ensure  of 
simi  lurity  7\^(i,j)  for  different  (i,j)'son  the  sanpling  points.  The 
figures  demonstrate  that  block  1  converges  faster.  In  block  1,  each 
additional  pixel  is  taken  as  far  as  possible  from  the  vicinal  pixel. 


D)  Thresholding  Sequence 

From  Fig.  3,  we  notice  the  clear  distinction  between  the  b^st  match 


L''AD-A13V  923  ACQUISITION  IMAGE  AND  DATA  COMPRESSION! U)  CITY  COLL  NEW 
YORK  D  L  SCHILLING  ET  AL .  30  APR  83  AFOSR- TR -83- 0505 
AFOSR-8 1  -  0 169 


UNCLASSIFIED 


F/G  17/2 


IV  Performance  compar  1  son  o  1  the  different  criterion  function 


We  have  compared  the  performance  of  different  measures  of  similarity. 
(These  results  have  been  shewn  in  Table  II.  It  turns  out  that  the  absolute 
difference  measure  is  faster  than  other  measures.) 

These  measures  of  similarity  are  as  follows: 

1.  Average  Absolute  Difference  Measure  (ABS) : 

M  M  | 

ABS(i,j)  =  —  ::  Sx](p,q)  -  T(p,q) 

M'  p=i  q=i  ' 

2.  Absolute  Difference  Measure  (AB)  : 

M  M 

AB(i,j)  =  ^  Slj(p,q)  -  T(p,q) 


Ob) 


I V  Per tonnance  caparison  of  the_ dif f crei it  cri  ter i on_  funct i on 

Wo  have  compared  the  performance  of  different  measures  oj  similarity. 
(These  results  have  been  shewn  in  Table  11,  It  turns  out  that  the  absolute 
difference  measure  is  faster  than  other  measures.) 

These  measures  of  similarity  are  as  follows: 

1.  Average  Absolute  Difference  Measure  (ABS) : 


ABS (i , j)  =  — 


M  M 

v  y 

p=i  q=i 


Slj(p,q)  -  T(p,q) 


2.  Absolute  Difference  Measure  (AE)  : 


AB(i, j) 


M  M 

JL  1  E 

4M  p=i  q=i 


Slj  (p,q)  -  T(p,q) 


3.  Correlation  Coefficient  Measure  (COE)  : 

M2o2  -  o4o5 

COE(i,j)  =  _  _ ;  - 

|  M^-  om21/M2o3-  o52 


where  0] ,  o  ,  ojf  o4  and  o s  as  given  in  equation 


4.  Normalized  Correlation  Function  (NCF)  : 

M  M  ±. 

EES  J(p,q)T(p,q) 

NCFU.j)  =  - 

MM-. 

E  E  rS1D(p,q)f 

p=i  q=i|- 

5.  Direct  Correlation  Function  (DCF) : 

M  M 

1  E  E 

p=i  q=i 


DCF(i, j) 


Slj(p,q)  T(p,q) 


The  results  of  the  conparison  between  the  five  measures  of  similarity  are 
shown  in  Table  II.  It  turns  out  that  the  ABS  measure  gives  the  best  match 
between  the  template  and  the  scene.  The  COE  is  time  consuming.  The 
thresholding  sequence  f  (b)  has  been  used  as  in  reference  (1),  where  f(b) 
is  an  errpirically  derived  function. 

f  (b)  =  (b:’  -  0.25b  -  0.25)/  b(b  4-  1) 

The  NCF  measure  is  slightly  slower  than  the  COE  measure.  DCF  function  is 
the  worst. 
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V  Conclusion 

In  this  paper  a  fast  algorithm  for  image  matching  has  been  studied. 
The  image  similarity  measure  used  is  the  sum  of  the  absolute  intensity 
difference  between  the  template  and  the  scene.  The  thresholding  sequence 
has  been  found  to  be  a  constant  and  equals  approximately  50%  of  the 
average  intensity  per  pixel. 

The  performance  of  the  algorithm  has  been  tested  by  computer 
generated  images  (English  letters)  and  Gaussian  noise. 

Comparison  of  the  algorithm  with  other  measures  of  similarity  proves 
that  the  average  absolute  intensity  difference  is  faster  than  others.  We 
intend  to  investigate  the  algorithm  performance  with  real  images  and  at 
real  time. 
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Block  1 


Fig.  1  Sampling  manner  of  the  data  set 
(each  nass  requires  that  one  additional  pixel  he  chosen  from  each  block  ) 
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